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Scope of JET PROPULSION 


Jet Propu.sion, the Journal of the American Rocket Society, is 
devoted to the advancement of the field of jet propulsion through the 
publication of original papers disclosing new knowledge and new de- 
yelopments. The term ‘‘jet propulsion’’ as used herein is understood 
toembrace all engines that develop thrust by rearward discharge of a 
jet through a nozzle or duct; and thus it includes systems utilizing 
atmospheric air and underwater systems, as well as rocket engines. 
Jer PROPULSION is open to contributions, either fundamental or ap- 
plied, dealing with specialized aspects of jet and rocket propulsion, 
such as fuels and propellants, combustion, heat transfer, high tem- 
perature materials, mechanical design analyses, flight mechanics of 
jet-propelled vehicles, astronautics, and so forth. JeT PROPULSION 
endeavors, also, to keep its subscribers informed of the affairs of the 
Society and of outstanding events in the rocket and jet propulsion 


feld. 
Limitation of Responsibility 


Statements and opinions expressed in JET PROPULSION are to be 
understood as the individual expressions of the authors and do not 
necessarily reflect the views of the Editors or the Society. 


Subscription Rates 


One year (twelve monthly issues)..................00005 $12.50 

Foreign countries, additional postage............... add .50 


Change of Address 


Notices of change of address should be sent to the Secretary of the 
Society at least 30 days prior to the date of publication. 


Information for Authors 


Preparation of Manuscripts 


Manuscripts must be double spaced on one side of paper only with 
wide margins to allow for instructions to printer. Submit two copies: 
original and first carbon. Include a 100-200 word abstract of paper. 
The title of the paper should be brief to simplify indexing. The 
author’s name should be given without title, degree, or honor. A 
footnote on the first page should indicate the author’s position and 
affiliation. Include only essential illustrations, tables, and mathe- 
matics. References should be grouped at the end of the manuscript; 
footnotes are reserved for comments on the text. Use American 
Standard symbols and abbreviations published by the American 
Standards Association. Greek letters should be identified clearly for 
the printer. References should be given as follows: For Journal 
Articles: Authors, Title, Journal, Volume, Year, Page Numbers. 
For Books: Author, Title, Publisher, City, Edition, Year, Page 
Numbers. Line drawings must be made with India ink on white 
paper or tracing cloth. Lettering on drawings should be large 
enough to permit reduction to standard one-column width, except 
for unusually complex drawings where such reduction would be pro- 
hibitive. Photographs should be clear, glossy prints. Legends must 
accompany each illustration submitted and should be listed in order 
on a separate sheet of paper. 


Security Clearance 


Manuscripts must be accompanied by written assurance as to 
security clearance in the event the subject matter of the manuscript 
is considered to lie in a classified area. Alternatively, written assur- 
ance that clearance is unnecessary should be submitted. Full respon- 
sibility for obtaining authoritative clearance rests with the author. 


Submission of Manuscripts 


Manuscripts should be submitted in duplicate to the Editor-in- 
Chief, Martin Summerfield, Professor of Aeronautical Engineering, 
Princeton University, Princeton, N. J. 5 


Manuscripts Presented at ARS Meetings 


A manuscript submitted to the ARS Program Chairman and 
accepted for presentation at a national meeting will automatically 
be referred to the Editors for consideration for publication in JET 
Proputsion, unless a contrary request is made by the author. 
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For further information, please write to: 
Mr. Mark Peters 

TECHNICAL RECRUITING, BLDG. 100 
CINCINNATI 15, OHIO 


Mr. William E. Potter 
ROCKET SECTION, MALTA TEST STATION 
BALLSTON SPA, N. Y. 
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The performance records of rockets as power 
plants for supersonic missiles and aircraft have 
led to accelerated design and development 
activities in this field. Our huge satellite program 
is now under way. Engineers who wish to 
contribute to advances in propulsion will be 
interested in the following openings created 

by expansion at GE. 


ROCKETS 


Design and develop rocket engine valves, seals and 
piping. Establish and maintain standards and 
instructions for their manufacture, test, 
installation and operation. 


e Design, develop, construct and test rocket engine 
seals and piping. Prepare standards and 
instructions for their manufacture, test and use. 


e Assemble and test rocket engines from their 
functional components. Design and 
develop supporting members needed to 
mount engine components. 


e Design, develop, construct and test 
propellant pumps. 


© Design, develop, construct and test turbines. 
Conduct performance analysis. 


Create advanced inlet designs for high 
performance rockets ; provide data on 
aerodynamic design and performance of such 
inlets, and conduct theoretical and experimental 
investigations of designs created. 


You'll be working with top men in the field, 

in small groups which stimulate exchange of ideas. 
Here, your talents can be easily recognized and 
rewarded. Educational assistance programs and 
specialized technical courses are also 

available to you. 


Openings in both 
MALTA, NEW YORK 


and 
CINCINNATI, OHIO 
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Are You One of the 728 Engineers Who Needs This 
HYDRAULIC PUMP 
nd 
This New Series 
HYDRAULIC PUMPS 
Answers Many of Your Demands for... 
\ 
\ The quiet efficiency of this new 65W Series of STRATO- LIGHTER WEIGHT 
' POWER Variable Delivery Pumps excites the interest of 
design and project engineers because this efficiency spells . 
SMALLER ENVELOPE 
Again, STRATOPOWER has come forward with a signifi- 
; cant development in hydraulic equipment, geared to your # 
advanced thinking. Compacted into a smaller envelope, 
and with a remarkably low weight/horsepower ratio, HIGHER SPEEDS 
STRATOPOWER 65W Pumps operate at system pressures 2s 
to 3000 psi and the smaller sizes at 10,000 i continuous 
speed! HIGHER ALTITUDES 
Whether your project concerns jet propelled aircraft, & 
rockets or guided missiles, consider STRATOPOWER 
Hydraulic Pumps as the heart of your Hydraulic system. HIGHER PRESSURES 
WATERTOWN owision 
THE NEW YORK AIR BRAKE cours’ (fy) HIGHER TEMPERATURES 
STARBUCK AVENUE ° WATERTOWN N.Y. 
INTERNATIONAL SALES OFFICE, 90 WEST ST., NEW YORK 6, N. Y. 
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0C-7C 
SEVEN SEAS 


ANNOUNCING 


A CRUISE CONTROL SYSTEM 
WITH NO MOVING PARTS... 


The Giannini-Douglas Differential 
Pressure Flight Angle Computer 


Differential pressure probes with no 
moving parts sense cruise control 
air data on Douglas Aircraft’s new- 
est overseas transport, the DC-7C. 
The probes are an integral part of 
a completely new angle-of-attack 
measuring system, the Giannini- 
Douglas Differential Pressure 
Flight Angle Computer, which was 
designed to have the greatest relia- 
bility and longest operating life pos- 
sible in a cruise control sensing 
instrument. 

In the Giannini-Douglas developed 
unit, small impact probes are accu- 
rately positioned on the head of a 
short stub boom mounted on the 
side of the fuselage; or for flight test, 


G.M. GIANNINI & CO.,INC. 
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on a free airstream boom. The 
probes are connected to sensitive 
Giannini pressure transducers 
which supply electrical signals pro- 
portional to air data to a passive 
network computer having no vacu- 
um tubes. The output of this integ- 
rating unit can be fed directly into 
an automatic flight control system 
or can be used to activate a panel 
indicator. 

System accuracy in the control 
range is +0.1° to +0.2° and angle 
of attack data can be sensed over a 
range of +20° from Mach 0.3 to 2.0. 
Less than 0.25 ampere at standard 
aircraft voltage is required for con- 
tinuous operation. 


918 EAST GREEN STREET 


Thoroughly proven in wind tunnel 
and flight test, the Giannini-Doug- 
las Differential Pressure Flight 
Angle Computer is one more out- 
standing example of recent ad- 
vances in aeronautical progress 
made possible by the ingenuity and 
skill of today’s research and design 
engineers in the field of avionics. 


AIRBORNE SYSTEMS DIVISION 
REGIONAL SALES OFFICES 

NEW YORK 1, N.Y., Empire State Bldg., CHickering 4-4700 

CHICAGO, III., 8 So. Michigan Ave., ANdover 3-5272 

PASADENA, Calif., 918 E. Green St., RYan 1-7152 
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THE SKY IS FULL 
OF VITAL 
CLARY CONTROLS 


.../n missiles, rockets 
and aircraft / 


LIKE SO MANY missile developers, the 
designers of the Firestone “Corpo- 
ral,” the North American “Navaho,” 
and four other classified missiles, 
specified Clary for vital automatic 
controls. Every control that bears 
the Clary name combines delicate 
precision (up to seven-one-millionths 
of an inch tolerance) with bulldog 
construction and supreme reliability. 
Our long experience in manufactur- 
ing to the stringent requirements of 
missiles also gives us an outstanding 
background in producing vital 
rocket and aircraft components... 
for air-to-air rockets, jet engine 
after-burner actuators, fuel system 
emergency valves and others that 
cannot be named. 
WE INVITE You to come to Clary, 
where our expanding engineering 
and manufacturing facilities provide 
you with— 

» DESIGNING AND TESTING to 

established specifications and 

envelope drawings. 

PRODUCTION-ENGINEERING 


of parts or complete components covered 
by your prototype sketches or drawings. 
MANUFACTURING of precision 
components to established drawings 
and specifications. 


TYPICAL CLARY AUTOMATIC CONTROLS 


Servo actuator Gate-type propellant 
for the ‘‘Corporal’’ valve for the ‘‘Navaho”’ 


Gyroscope After-burner actuator 
for the ‘‘Corporal’’ for jet engines 


AUTOMATIC CONTROLS 
DIVISION 


CLARY CORPORATION 
Dept. J46, San Gabriel, Calif. 
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DESIGNERS 


for MISSILE SYSTEMS 


@ New activities at Lockheed Missile Systems 
Division have created positions for Designers 


VARIED ASSIGNMENTS 


capable of performing creative basic layout 
and design of structural, mechanical, 
electro-mechanical and electronic packaging 


Designers are not limited to 
specific functions under Missile 
Systems Division’s philosophy 


of operation. Diversified 
assignments provide stimulating 
challenge, enable Designers 

to acquire the broadest possible 
background in the field of 
missile systems design. For 
example: it is not uncommon 
for a Missile Systems Designer 
to work on diverse problems 

in structures, controls, 
hydraulics, pneumatics, electro- 
mechanical packaging, fuel 
systems and other fields within 
the span of a few months. 


‘of missile assemblies and components. 


@ Those who will qualify will cope with new 
problems in a field of scientific endeavor that 
grows daily in complexity. A knowledge 
of new materials, finishes, specifications and 
experience on small precision devices will 
prove helpful in meeting the challenge 
of Missile Systems research and development. 


Those possessing a high order of ability applicable 
to these areas of endeavor are invited to write: 


MISSILE SYSTEMS DIVISION 
4 research and engineering staff 


LOCKHEED AIRCRAFT CORPORATION + VAN NUYS, CALIFORNIA 
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ENGINEERS... 
LOOK TEN YEARS AHEAD! 


Take that ten year ahead look. There’s 
a fine career opportunity in the engi- 
neering field you like best waiting for 
you at Douglas. 

And what about the Douglas Aircraft 
Company? It’s the biggest, most suc- 
cessful, most stable unit in one of the 
fastest growing industries in the world. 
It has giant military contracts involv- 
ing some of the most exciting projects 
ever conceived...yet its commercial 
business is greater than that of any 
other aviation company. 

The Douglas Company’s size and 
variety mean that you'll be in the 


DOUGLAS 


Will your income and location 

allow you to live in a home 

like this...spend your 
leisure time like this? 


They can...if you start your 
career now at Douglas! 


A Douglas engineer lives here 


work you like best — side by side with 
the men who have engineered the finest 
aircraft and missiles on the American 
scene today. And you’ll have every 
prospect that ten years from now you'll 
be where you want to be career-wise, 
money-wise and location-wise. 

For further information about opportunities 
with Douglas at Santa Monica, El Segundo and 
Long Beach, California and Tulsa, Oklahoma, 
write today to: 

DOUGLAS AIRCRAFT COMPANY, INC. 


C. C. LaVene, 3000 Ocean Park Blvd. 
Santa Monica, California 
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One development agency for the Air Force’s ground- 
based electronic equipment is Rome Air Develop- 
ment Center, located at Griffiss AFB in Rome, N. Y. 
One of the ten centers of the Air Research and 
Development Command, RADC is concerned with 
the air defense of our nation, with providing equip- 
ment for tactical supremacy, and with developing 
ground complexes for various navigation systems to 
aid all aircraft. In addition, RADC is charged with 
data handling improvements for the Air Force 
intelligence mission. 

RADC is the responsible center for development, 


84 
FORD INSTRUMENT COMPANY 


DIVISION OF SPERRY RAND CORPORATION 
31-10 Thomson Ave., Long Island City 1, N. Y. 


This is one of a series of ads on the technical 


activities of the Department of Defense. 
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Engineers at Ford Instrument Company 
working on a special Air Force project 
in one of the company's laboratories. 


Official U.S. Air Force Photograph 
Friend or foe? Tactical defense officers in MINK control center watch movements of aircraft as reported from 
radar warning sites. This is equipment developed under the direction of Rome Air Development Center. 


ROME AIR DEVELOPMENT CENTER PROVIDES 
AIR DEFENSE EQUIPMENT FOR OUR AIR FORCE 


through its various contractors, of such end products 
as radar sites, including improved tubes, circuits, 
antennae, and shelters; ground communications 
equipment and associated support items; IFF 
(Identification Friend or Foe) environments, and 
electronic countermeasures. 

Bringing complex systems from the written re- 
quirements to the actual hardware items to be used 
in the various Air Force commands is a long and 
tedious business which draws upon the skills of 
RADC’s 500 civilian and military engineers and 
their many counterparts in private industry. 
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WHY ENGINEERS 


FEEL AT HOME 


AT ROCKETDYNE 


First and foremost, RocKETDYNE talks your lan- 
guage—and understands it too. Your associates and 
supervisors here are professional people like you. 
They respect your status, your thinking, your ideas 
and your interest in technical advancement. 

ROcKETDYNE will encourage you to choose the 
field that is most satisfying and rewarding... truly 
best for you. This is possible because its activity 
includes the full range of rocket engine development 
from preliminary design to field testing and produc- 
tion... because its programs include development 
of the largest liquid-propellant rocket engine in the 
Western World... because it has contracts with all 
branches of the Armed Services and the guided mis- 
sile industry for broad variety of rocket engine types 
and sizes. 

It may surprise you to know you can qualify for 
a career at ROCKETDYNE with or without specific 
rocket engine experience ! Engineering experience in 
heating and ventilating, hydraulics, pumps, tur- 


THESE POSITIONS NOW 
DESIGN & DEVELOPMENT ENGINEERS 
Mechanical, Chemical, Electrical, Aeronautical, Stand- 
ards, Structural and Stress. For rocket engine com- 
ponents and systems design or development. Turbine, 
pump, controls and combustion device experience pre- 
ferred. 
TEST ENGINEERS 
Experienced on engine systems, combustion devices, 
turbines, pumps and engine instrumentation. 
EQUIPMENT DESIGN ENGINEERS 
Electrical, mechanical, structural, industrial. For design 
of facilities, specialized test, and handling equipment. 


Write] Mr. Grant Baldwin, Rocketdyne Engineering Per- 
sonnel, Dept. 596-JP, 6633 Canoga Ave., Canoga Park, Calif. 


bines, combustion devices, controls, dynamics, 
structures and instrumentation are just a few of 
the related fields that could open your future at 
ROCKETDYNE. 

ROCKETDYNE’Ss design and manufacturing center 
and its nearby test laboratory house complete, 
advanced facilities... the vital tools you need to 
meet the challenges of rocket engine development. 

RocKETDYNE is North American’s rocket engine 
division. It has just moved into new ultra-modern 
headquarters in Canoga Park, located in the beauti- 
ful West San Fernando Valley of Los Angeles. This 
area is famous for its fine residential sections, 
modern shopping-center convenience, varied recrea- 
tional and entertainment facilities. Any point in the 
San Fernando Valley is just minutes drive from the 
beaches, and the weather is pleasant all year around. 
Many engineers are interested in advanced courses 
offered by fine schools like UCLA, USC and Cal 
Tech, all within a short drive from our headquarters. 


OPEN AT ROCKETDYNE: 


DYNAMICS ENGINEERS 

To analyze rocket engine control systems utilizing elec- 
tronic analog and digital computers, B.S., M.E., or 
B.S.E.E. necessary. Prefer advanced degree. Experi- 
ence in servomechanisms, systems analysis desired. 
THERMODYNAMICIST 

To analyze, design and develop high speed subsonic 
and supersonic turbines. Jet engine or industrial steam 
turbine experience desired. 

SYSTEMS ENGINEERS 

ENVIRONMENTAL TEST ENGINEERS 


STANDARDS ENGINEERS 
TEST ENGINEERS — ELECTRICAL 


ROCKETDYNE 


A DIVISION OF NORTH AMERICAN AVIATION, INC, 
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ae -scale ‘‘on-location” oxygen production is 
economically essential to major rocket research 
and operational facilities. Air Products Oxygen 


oS, Generators provide the most dependable source of 
of oxygen in unlimited quantities at lowest cost. This 
at statement is amply supported by experience based 
on many successful installations. 
er We will provide and install an ‘‘on-location” oxygen 
e, and/or nitrogen generating station to suit your pres- 
to ent and expanding oxygen requirements. This can 
t. be done on a lease basis, without capital invest- 
ne ment by you. Operation and maintenance can be 
n handled by you or by Air Products—at your option— 
i- with price guaranteed at all levels of consumption. 
1S We design and manufacture generators for unlimited 
8, quantities of oxygen and nitrogen regardless of size, 
a purity or cycle. 
és If you will advise us regarding your requirements, we 
i will be glad to furnish you with actual costs and 
‘ engineering information on oxygen generators de- 
] signed for your particular needs. 


~ 
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Cost Analyses  Procés Design Apparatus Manufacture 
Air Products 


INCORPORATED 
Dept. O, Box 538 Allentown, Pa. 
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Precision Mechanics, Optical Devices, Ceramics ee oe 
Electronics oe eee @| 

Hydraulics, Liquids Processing, Heat Exchange ee eee ry.) 

Television: Studio, Theatre, Institutional, Business, Industrial e 

Instruments, Servos, Controls: Hydraulic, Pneumatic, Magnetic, Electronic @| 
Aircraft and Missile Guidance, Control, Simulation eee eee % 

Automatic Computers and Components e| @ 

Radar, Microwave, Ultrasonics @| ee 
Motion Picture and Audio Equipment eee 06 eee 

Nuclear Power Components and Controls ee oee eG 

Systems Engineering: Aeronautical, Naval, Industrial 

1 Kearfott Company, inc.; Little Falls, New Jersey 
lf International Projector Corporation; Bloomfield, New Jersey itt Bludworth Marine; New York 
IV General Precision Laboratory Incorporated; Pleasantville, New York V_ The Griscom-Russell Company; Massillon, Ohio 
Vi Link Aviation, inc.; Binghamton, New York Vil Shand and Jurs Co.; Berkeley, California Vill The Hertner Electric Company; Clevelana 
IX The Strong Electric Corporation; Toledo, Ohio X Precision Technology, Inc.; Livermore, California Xt J. E. McAuley Mfg. Co.; Chicago 
IL Askania Regulator Company; Chicago xi Ampro Corporation; Chicago 
XIV Librascepe, incorporated; Glendale, California 


@ Manufacturing @@ Manufacturing and product oe ig, product 


and research 
OS® Pilot manufacturing, product development and research 
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More than a year ago, one of the operating companies of General Precision 
Equipment Corporation—after nearly 8 years of joint development work with the 
Air Force—placed in scheduled production the most advanced air navigation 
system known to exist. Called AN/APN-66, the system is based on “Doppler effect,” 
a known natural phenomenon, but one which others said could never be applied 
to aviation engineering. Its successful development brought to culmination still 
another notable aviation engineering “first” by a GPE Company. 

AN/APN-66 and related systems adapted for specific purposes compute 
air position with unprecedented accuracy and speed. They work independently 
of any ground-based reference, optical or electronic; are the only operational 
systems that are global. They provide a vitally needed element in the defensive- 
offensive power our country must have to speak softly, yet negotiate from strength. 

For months now, these systems have been functioning as brain and nerve 
center in a substantial number of the Air Force’s newest craft. In millions of miles 
of operational flight they have demonstrated their engineering genius. When put to 
civilian use, they will guide air liners to the remote corners of the globe with equal 

2 accuracy. Then, never again need a plane get lost anywhere...in any weather. THT 
While the GPE Company primarily responsible for this AN/APN equipment 
is General Precision Laboratory, the systems are typical products of GPE Coordi- 
nated Precision Technology. Four other GPE Companies—Askania, International 
4 Projector, Kearfott and Librascope—participated with GPL in providing the com- 
ponents that made them possible. All now manufacture portions of the equipment. 

As the chart on the opposite page shows, all of the producing companies in 
the GPE Group work in‘highly advanced fields. GPE Coordinated Precision Tech- 
nology, a basic GPE operating policy, gives each GPE company access to the skills 
and facilities of the others and coordinates those which are relevant. This coordina- 
tion has played an important part in amplifying the capabilities of all GPE com- 
panies, enabling them to solve highly advanced technological problems with vision. 

Frequently, as with these AN/APN systems, frontiers of science are pushed back 
on the way. 

The GPE Producing Companies serve more than a dozen major industries. 
A brochure, “Serving Industry through Coordinated Precision Technology,” de- 
scribes their fields in some detail. It may well suggest products or facilities appli- 
cable to problems your company is facing. For a copy or other information, write: 


General Precision Equipment Corporation 
92 Gold Street, New York 38, New York 


| 
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BARBER 
COLMAN 


Aireraft Controle 


help keep Boeing IM-99 missile accurately on course 


The Boeing Bomarc, a supersonic pilotless 
interceptor being developed for the USAF, seeks and 
destroys enemy bombers, while still distant from vital 
targets. Electronic devices, including Barber-Colman 
tachometer generators, play an important part in 
controlling its flight path toward the target. These 
generators meet or exceed -all of the requirements 
specified by the missile manufacturer. 


Artist's conception of Boeing Bomare 
The Barber-Colman tach generator used on this missile applica- 
tion is a special adaptation of the widely known Barber-Colman 


line of high quality d-c motors for aircraft. 


The complete line of Barber-Colman aircraft controls includes: 
Actuators; Positioning Controls; Temperature Controls; Small 
Motors; Valves; Ultra-Sensitive Relays; Thermo-Sensitive Ele- 
ments. Write today for literature or consult the Barber-Colman 
engineering sales office nearest you . . . Los Angeles; Seattle; 
Baltimore; New York; Montreal, Melbourne. 
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Here’s how a new era in small jet development was pioneered. 


More than seven years ago, Fairchild Engine Division began 
work on the mighty J44, a power-packed, bantam-weight 
turbojet, to power guided missiles, target drones and pilot- 
less planes. 

Its 150-hour military qualification tests completed, the J44 
can now be put to work as thrust assist on multi-engined 
military cargo and transport aircraft; or as part of auxiliary 
power and boundary layer control systems. 


From the beginning, Fairchild engineers have continually de- 
signed increased performance, reliability and ruggedness into 
the J44. Faster speeds, higher altitudes were achieved, and 
with them, a wider and wider range of new applications. 
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Hypersonic Studies the Leading Edge Effect the 
Flow Over a Flat Plate 


A. G. HAMMITT' and S. M. BOGDONOFF’? 


Forrestal Research Center, Princeton University, Princeton, N. J. 


Using the Princeton helium hypersonic wind tunnel, the 
dlect of leading edge thickness on the flow over a flat plate 
has been investigated. The study was carried out by opti- 
@l measurements of the shock position over a wide range 
@ leading edge thicknesses and surface pressure measure- 
ments for the thicker leading edges. The flow over the 
plate with a thick leading edge was found to be ntially 
inviscid, but viscosity was important for the case of thin 
leading edges. The strong inviscid effect caused by the 
finite thickness of the leading edge seemed to be sufficient 
to explain the discrepancies between viscous theory and 
experiment in the region near the leading edge. 


Nomenclature 


= linear viscosity coefficient in relation p/m = C(7T.,/T;) 

boundary layer parameter depending on surface tempera- 
ture 

Mach number 

pressure 

static pressure increment 

Reynolds number 

leading edge thickness, inches 

temperature 

distance from leading edge 

distance of shock above plate surface 

velocity divided by maximum velocity on expanding to 
zero pressure 

specific heat ratio 

boundary layer displacement thickness 

viscosity 


nou 


Subscripts 
tunnel empty condition 


stagnation 

based on leading edge thickness 
condition on wall 

based on distance from leading edge 


Introduction 


HE analysis of flows at many times the speed of sound are 

complicated by special problems which are unimportant 
at lower speeds. These problems can be broken down into 
two classes: (a) fluid dynamics problems, which depend on 
the conventional parameters of Mach number, Reynolds 
Bumber, Nusselt number, and specific heat ratio; (b) mole- 
tilar problems, which depend on temperature and involve the 
thange of molecular structure at high temperatures. 


Presented at the ARS 25th Anniversary Annual Meeting, 
Chicago, Ill., Nov. 14-18, 1955. 

Research Associate, Gas Dynamics Laboratory. 

* Associate Professor. 


* Numbers in parentheses indicate References at end of paper. 
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The hypersonic fluid dynamic problems are governed by 
the same relations as other fluid dynamic problems. Con- 
siderable study has been given to the predominant fluid 
dynamics problems in the subsonic and supersonic range. 
At hypersonic speeds the relative magnitude of the various 
effects may be considerably different. This change in 
emphasis leads to several new problems which have received 
little attention. Effects of little importance aft supersonic 
speeds may become dominant at hypersonic speeds, and vice 
versa. 

The molecular problems are essentially new in aerodynamic 
applications. These problems only occur at high tempera- 
tures where the molecular structure may be considerably 
different than at conventional temperatures. Some of these 
problems have been encountered at the relatively low speeds 
and high temperatures which exist in rocket motors. 

This paper describes a study of the flow about the leading 
edge of a flat plate and the experimental equipment used for 
this study. This flow involves several hypersonic fluid 
dynamics problems. At hypersonic speeds most shocks are 
strong and the entropy changes through the shocks cannot be 
neglected. The reflected waves, caused by the interaction 
of other waves with the strong shocks, are also important. 
Viscous forces are large and can cause-important modifications 
in the inviscid external flow. All these factors combine in the 
flow about the leading edge of a flat plate. The study of this 
problem was carried out in a hypersonic wind tunnel using 
helium as the working fluid. The tunnel operates in a Reyn- 
olds number range where the fluid is essentially a continuum 
and slip phenomena are limited to very small distances. 
The work reported herein covers an analysis of shock shapes 
and pressure distributions of the hypersonic flow about a flat 
plate over a wide range of Reynolds numbers and leading 
edge thicknesses. 

This work was carried out in the Gas Dynamics Labora- 
tory of the Forrestal Research Center, Princeton University, 
and was sponsored by the United States Air Force, Aero- 
nautical Research Laboratory, Wright Air Development 
Center, Wright Patterson Air Force Base, under contracts 
numbers AF 33(038)-250 and AF 33(616)-2547. More de- 
tailed descriptions of the helium tunnel and research are 
contained in (15)* to (20). 


Hypersonic Wind Tunnel 


General Discussion 


There are several techniques for making fluid dynamics 
experiments which have been proved useful. Wind tunnels, 
shock tubes, and free flight experiments all have important 
uses. The wind tunnel has proved to be very valuable for 
many investigations in the subsonic and. supersonic range and 
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might be expected to be equally useful in the hypersonic 
range. The wind tunnel allows tests to be made over rela- 
tively long periods, under controlled conditions, with ease 
and nominal expense. 

The main problem in operating a wind tunnel hypersoni- 
cally is to heat the air enough so that it does not condense after 
expanding to high velocities. Stagnation temperatures of 
the order of 1300 to 1500 F are required to prevent condensa- 
tion of the air in the test section at Mach number 10. To 
obtain stagnation temperatures high enough to study the 
molecular problems of hypersonic flow, of the order of 4000 to 
6000 F, seems to be impractical in a wind tunnel. For this 
reason hypersonic wind tunnels are limited to a study of only 
the fluid dynamic problems. Another technique for avoiding 
condensation in the working section of a wind tunnel is to 
use a gas which does not condense. Since the wind tunnel 
cannot simulate the full flight problem, there is little loss for 
the purpose of fluid dynamics research in using a gas other than 
air. If a gas with low condensation temperature such as 
helium is used, Mach numbers in excess of 20 can be ob- 
tained at atmospheric stagnation temperatures. The elimi- 
nation of the heater and the accompanying problems of a 
hot tunnel is such an important simplification that the helium 
tunnel seems to be the most satisfactory device for fluid 
dynamics research at Mach numbers greater than 10. The 
results obtained in the helium tunnel do not represent the 
complete answer to hypersonic flight but only lead to an 
understanding of the fluid dynamics problems. The data 
obtained from such a tunnel are directly comparable to the 
data obtained in hot air tunnels at the same Mach, Reyn- 
olds, and Nusselt numbers, corrected for the effect of the 
different specific heat ratio. 


Princeton Helium Hypersonic Tunnel 


The Princeton hypersonic wind tunnel is of the blow-down 
type using helium as the test fluid. Fluid dynamics problems 
may be studied in this tunnel up to Mach 20. A detailed 
description of the tunnel is given in (15). The helium is ob- 
tained at high pressures in tank semitrailers which are used as 
the storage system. From the trailer, the helium is piped 
through a regulating valve to the settling chamber. The 
tunnel may be run at stagnation pressures up to 2000 psia for 
periods of about 10 minutes. A stagnation pressure range of 
800 to 1400 psia has been used in these tests. To obtain the 
high pressure ratios which are needed, an air operated ejector, 
using air at 1500 psi from a second high pressure storage 
system, is used downstream of the hypersonic diffuser. A 
general layout of the piping and tunnel is shown in Fig. 1. 


TO HELIUM SUPPLY 
= 
VALVE 
REDUCING 
GAGE BOARD 
CONTROL 


Fig. 1 Piping and control system for helium hypersonic tunnel 


The first nozzle used, Fig. 2(a), was an axial symmetric 
conical nozzle with no windows and a maximum diameter of 
approximately 3 inches. The second nozzle, Fig. 2(b), is 
similar to the first one except that it has been cut open and 
windows installed. The disturbances from the windows do 
not effect the central core in the region of interest. The 
models and probe are supported in the test section by means 
of the sting and spider shown. The first nozzle was used on 
the initial surface pressure distribution tests and the second 
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Fig. 2. Cross sections of hypersonic tunnel 


Top: axial symmetric nozzle; center and bottom: axial symmetric 
nozzle with windows. 


Pag 


2 B 4 16 17 18 9 20 
Distance from throat, inches 


Fig. 3 Mach number distribution along the center line of the 
hypersonic tunnel 


nozzle for the optical studies and additional pressure distribu- 
tion measurements. 

The Mach number along the center line of the nozzle is 
shown in Fig. 3. Since the nozzle was not contoured, a re- 
gion of constant Mach number was not obtained and the tests 
are made in a region of Mach number gradient. The effects 
of this gradient will be discussed later. The diffuser, Fig. 2, 
is a straight tube with a central body sting support. This ar- 
rangement is probably not an optimum, but gave the best 
results of the limited varieties which were examined. 


Experimental Study 


A test program was undertaken to study the flow about the 
forepart of a flat plate at high Mach number. {If only the 
adiabatic case is considered, there are three parameters which 
might influence the flow. These parameters are the Mach 
number, the Reynolds number, and the geometry. In the 
case of the semi-infinite flat plate, the only geometric parame- 
ter is the distance scale based on leading edge thickness. 
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The experimental model could only approach a semi-infinite 

te since it had finite boundaries. In general, it would be 
expected that the effects of the tips and trailing edges would 
be limited to small regions and would not propagate over the 
entire plate. It would be expected that there should be 
considerable area which behaved as a semi-infinite plate. 
Again judging by lower speed data, the effect of the leading 
edge would be expected to be limited to a few leading edge 
thicknesses. Therefore, a sharp flat plate of reasonable as- 
pect ratio might be expected to behave much like an in- 
finitely sharp semi-infinite flat plate. (This is the model 
treated by several theoretical investigators (1) to (11)). 

The experimental study undertaken was an analysis of the 


‘flow over a flat plate at variable leading edge thicknesses, 


Reynolds number, and Mach number. The effects of the 
tips and trailing edge were investigated in order to find a 
region which was free of these effects. The study was carried 
out by measurements of the shock shape and the surface 
pressure distributions. In order to study the effect of the 
leading edge thickness it was necessary to vary this thickness 
over a wide range. The leading edge thickness and stag- 
nation pressures were varied over the widest ranges possible. 
The leading edge range was limited by the requirement that 
the leading edge had to be thick enough to remain unchanged 
during the run and could not be so thick that the tunnel was 
stalled. The stagnation pressure was varied between 1400 
and 800 psia giving a Reynolds number per inch range of 1.15 
to 0.65 million. To obtain sharp leading edges, a simple 
model with no pressure taps near the leading edge was used. 
This model could be hardened and honed on both surfaces. 
The model is shown in Fig. 4. Interferograms were taken of 
the model at Mach 12.7 at varying leading edge thicknesses 
and stagnation pressures. The leading edge thickness was 
varied from 0.17 to 59.5 thousandths of an inch. The blunt- 
ing was accomplished by cutting off the leading edge per- 
pendicular to the top surface. The measurement of the 
leading edge thickness was made with a microscope and cali- 
brated eyepiece. A typical range of interferograms is pre- 
sented in Fig. 5: ‘“‘a” is for a thin leading edge, “‘b” a medium 
edge, and “‘c’” for a thick leading edge. ‘“‘a’”’ shows a shock 
wave starting at or very close to the leading edge and a pro- 
nounced boundary layer indicated by the large fringe shift 
between the shock wave and surface. “c’’ shows a strong 
detached shock and almost no indication of a boundary layer. 
“b” is an intermediate position showing a blending together 
of the features found in “a” and “ec.” “a” shows the fringes 
approaching the plate perpendicularly. This does not mean 
that the flow is separated, but only that the density is so low 
at this point that the interferometer has not sufficient sensi- 
tivity to show the density changes in this region. A similar 
situation exists in “c.” In this case, the density has been 
reduced to such a level outside of the boundary layer by the 
entropy changes through the strong shocks that the inter- 
ferogram does not show the additional density changes 
through the boundary layer. These interferograms can be 
used only as qualitative indications of the density since the 


Leading edge cut off square for bunt fests 


Fi 


Fig. 4 Fiat plate model 


8. Rey = 177 


c- Re, = 15,010 


Fig. 5 Interferograms of flat plate mode. M = 1.27 


flow field is not completely two-dimensional (the model does 
not span the tunnel) and the fringe shifts are too small to be 
accurately determined. Only the position of the shock 
can be determined and quantitatively analyzed. 

Surface pressures were measured for the thicker leading 
edges. Pressure measurements were made on leading edge 
thicknesses between 1 and 43.8 thousandths of an inch. 
These measurements were made in two series, the first 
series at 0.001, 0.004, 0.007 inches leading edge thickness 
were made in the first nozzle at M = 10.75 and 11.8 and P, of 
800, 1000, and 1400 psia. The leading edges were not meas- 
ured with the microscope and are not known with the same 
exactness as in the other tests. The second series of tests were 
for thicknesses between 0.0176 and 0.0438 inches and were 
made in the second nozzle with the use of the microscope for 
leading edge measurements. Mach numbers of 11.4, 12.7, 
and 13.8 were run at P, = 800, 1000, and 1400 psia. 

The model used for these measurements was similar to the 
one shown in Fig. 4, but contained pressure orifices along its 
top surface. Several orifices were placed off center line to 
check the tip effects. The tip effects were found to be negli- 
gible over the central portion of the plate. A pressure rise 
was found near the back of the plate which could be traced to 
the forward propagation of the shocks around the sting sup- 
port. Only the data which were not influenced by the finite 
dimensions of the plate are presented. For most of the tests a 
square leading edge shape was used. A few of the thicker 
tests were run with a rounded leading edge to find the effects 
of the detailed leading edge shape. 


Analysis 


The positions of the shock waves as shown by the inter- 
ferograms were read on an optical comparitor. These posi- 
tions, reduced by the leading edge thickness, are shown in 
Fig. 6. This figure shows that the thicker leading edge 
shock shapes correlate very well on this plot while the thinner 
leading edge shapes show some deviation. The difference 
becomes pronounced below an Re, of about 4000. These 
measurements were made in the Mach number gradient 
shown in Fig. 3. No correction was made to the shock shapes 
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Fig. 6 Shock shapes for flow about flat plates at variable 
leading edge thickness 
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Fig. 7 Surface pressure increment for flow about a flat plate 
at variable Re; 


so that these should only be compared with each other and are 
not the shape for flow at constant Mach number. The pres- 
sure distributions have also been plotted as a function of 
t/x. These results have been corrected to a first order, for 
the Mach number gradient, by using the difference between 
the pressure on the plate surface and the tunnel empty static 
pressure at the same point. This nondimensionalized dif- 
ference is plotted in Fig. 7 for P, = 1000 psia and M = 11.4 
and 11.8. Only these Mach numbers were used in this plot for 
the sake of clarity, but these results are typical of all the 
Mach numbers. No significant difference is found between 
the square and the round leading edge. From this figure 
it can be seen that all the pressure data, except for that of 
0.001-in. leading edge, correlate quite well. This result 
could have been anticipated on the basis of the shock shape 
studies that showed correlation on the t/z basis down to 
Re, of about 4000. 

Leading edge thicknesses greater than Re, of 4000 give a 
flow which is independent of Reynolds numbers. This 
flow is essentially inviscid and the boundary layer does not 
have a major influence. This conclusion is borne out by the 
interferograms which did not show the boundary layer for 
the thick leading edge, Fig. 5(c). At lower Reynolds number, 
this correlation no longer holds and boundary layer effects 
are more important. This result again checks with the inter- 
ferograms which showed a pronounced boundary layer for 
thin leading edges, Fig. 5(a). 

The shock shape and surface pressure distribution for the 
inviscid flow regime are caused by the interaction of the 
detached shock wave and expansion waves from around the 
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leading edge. The magnitude and downstream persistence 
of the leading edge effects is surprising. Maximum pressures 
of the order of 10 to 20 times free stream are reached at about 
2 leading edge thicknesses from the leading edge followed by 
the pressure decreasing as */ t/z for t/x < 0.25. The follow- 
ing empirical relation has been fitted to the decreasing part of 


the curve. 
t 
4? _ 0.016 
Pi 


This relation is based on the data at all Mach numbers and 
all P,’s. Itis shown plotted in Fig. 7. 

There seems to be no available theoretical treatment of this 
problem. Most of the simplified techniques used to cal- 
culate hypersonic pressure distributions relate the pressure 
to the local slope of the body. Theories such as those based 
on Newtonian flow, shock plus expansion, or tangent-wedge 
are of this type. All such theories will predict constant pres- 
sure over the flat plate surface, so are entirely inadequate for 
this problem. This problem is an extension of the transonic 
problem to hypersonic velocities. For this reason it may be 
called a “tri-sonic’”’ problem and a complete solution must in- 
volve the same considerations as the transonic problem 
complicated by the rotational terms introduced by the curved 
hypersonic shock wave. 

In (14), Bertram has suggested a method for handling this 
problem. This approach was to use the model shown in 
Fig. 8. This model avoids the subsonic part of the problem 
by assuming an attached shock with M = 1 directly behind 
the shock. The results calculated for these assumptions 
are shown in Fig. 7. The numerical values and shape of the 
curve differ from the experiments, but at least the model 
predicts a varying pressure. 


Fig. 8. Simplified flow model 


x 


A better understanding of the problem may be had by a 
study of the hodograph plane shown in Fig. 9 along with the 
same configuration on the physical plane. If the assumption 
is made that the flow becomes sonic just before the sharp 
corner, the contour A represents the flow along the solid 
boundary. A 90 deg expansion gives infinite Mach number 
for helium, so zero pressure would be expected just behind the 
corner. The conditions of the flow along the shock wave 
are given by the shock polar. The sonic point on the shock is 
the intersection of the shock polar and sonic circle. Some of 
the expansion waves from the corner will intersect the sonic 
line and some will intersect the shock wave. Those inter- 
secting the sonic line will reflect as compressions, while those 
intersecting the shock will reflect as expansions. The limiting 
characteristic, the one dividing the waves which intersect the 
sonic line and those which intersect the shock, is given by 
contour B. Therefore all the waves starting before B would 
reflect as compressions and all beyond B as expansions. 
This reasoning would lead to surface pressure distribution 
that starts at zero behind the sharp corner, rises through the 
reflected compressions, and then decreases through the re- 
flected expansions. These compressions and expansions are 
represented by the part of contour A along the 0 deg axis. 
The surface condition first moves toward the origin through 
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at Prandtl number of.1 for the adiabatic case d, = 0.599- 


nee 
out The slope of the boundary layer displacement thickness is 
by 
Ow- 
t of 
This slope is taken equal to the deflection of the free stream. 
The pressure-angle relation may now be taken from the 
{1} Prandtl-Mayer or tangent-wedge relation for the free stream. 
Both of these relations give the same results for the first 
nd terms if expanded in series. 
dé* 
al- dz 
Substituting the relation for dé*/dz into Equation [3] 
ige Ap _p-p 
es- yao VR [ 
for 
nic For the adiabatic case in helium, y = °/3 
be 
Ap M/C M/C 
— = 0.599 (y) (vy — 1) —— = 0.66 —“— 
1 V Rez V Rez 
If this expression is assumed to be correct for the infinitely 
Lis thin leading edge plate, the correct expression for the finite 
in leading edge plate must approach this solution for small Re, 
m and the inviscid expression for large Re,. The simplest rela- 
id tion which meets these requirements is a linear combination 
ns of Equations [1] and [4]. 
he = 
Ap le [c 
el —— = 0.016 M3 ~ 0.66 M3 a 
Pi Rez 
Fig. 9 Flow about a blunt flat plate in the physical hodograph = = 
_ 0.016 (: + 40 JE) (5) 
2 
the compression and then away from the origin to the free 
stream condition. Qualitatively this agrees with the pressure This linear combination neglects the effects of the inviscid 
measurements which show a peak value and then a decrease flow on the viscous flow, and vice versa. It should be ex- 
for high values of t/z. pected to be fairly accurate at the two extremes where one 
It is interesting to note that the model suggested in (14) term predominates. From this relation some idea of the 
gives a considerably different picture on the hodograph plane. relative magnitude of the two terms may be found. The 
In this instance, the conditions on the surface are given by value of Re, at which both effects are equal is a function of 
contour C. This model lacks the overexpansion and re- M(C ~ M**). For Mach number of the order of 12, both 
compression from the sonic line found for the more correct effects are equal at Re, = 400. All of the surface pressure 
case. measurements are at considerably higher Re,. If the pressure 
This description has neglected the vorticity of the flow be- data from the thinner leading edge cases, the only ones in 
hind the shock wave. The vorticity will introduce reflections which +~/C/Re, is significant, are reduced by the factor 
of the expansion and compression waves from the vortex 1 
e sheets. Therefore, Fig. 9 can be used only as a qualitative 1 + 40\/C/Re. and plotted against ¢/x, Fig. 10 results. 
picture and an aid in understanding the flow. Fig. 10 shows that the 0.001-in. leading edge data correlate 
Several attempts have been made to account for the effects fairly well with Equation [5]. The data at 0.004 and 0.007- 
of vi f the leading edge of a sharp flat plate. 
; Ritely sharp leading edge, The general approach to this cous effect at 0.004 and 0.007-in. leading edge is overestima 
, problem has been to take into account the effect on the ex- 
; ternal flow, and subsequent feedback into the boundary een 
layer, of the change of effective body shape caused by bound- Pe, 


ary layer displacement thickness. The method used for.a 
first approximation in (2) is relatively simple and shows the 
general method of attack. The boundary layer displace- 
ment thickness for free stream conditions on the surface of 
the plate is taken as 
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Fig. 10 Reduced surface pressure increment for flow about a 
flat plate at variable Re; 
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by the viscous correction. This result is not unexpected 
since the Mach number near the surface for the thicker lead- 
ing edge is much less than free stream. Therefore, the viscous 
theory, based on free stream Mach number, will overestimate 
the boundary layer buildup. The relative magnitudes of the 
viscous and inviscid effects are surprising. In most of the 
previous work on the leading edge effects on sharp bodies at 
hypersonic flow, the leading edge inviscid effect has not been 
recognized. If the Mach number dependence of C is sub- 
stituted into the relation for surface pressures, the relation 


becomes 
Ap t 1.45 
— = 0016 M?. {1 + [6] 


This result indicates that the inviscid leading edge effect is the 
truly hypersonic one (predominant as M — ©) and the vis- 
cous effect seems to be significant at lower Mach numbers and 
very sharp leading edges. Because these results have only 
been derived from experiments over a limited range of Mach 
number, the accuracy of these extrapolations is open to 
question. However, the large effect on surface pressures of 
the leading edge thickness, even at distances greater than 
1000 leading edge thicknesses, is most significant. 

There are little other experimental data to compare with 
these results. Some data have been presented in (12) and 
(14) at Mach numbers of around 6 or 7 in air. These data 
appear to show much less leading edge inviscid effect than 
the data herein presented. This result confirms the idea that 
viscous effect is more important at intermediate speeds. The 
specific heat ratio probably has some effect on the relative 


magnitude of these effects and may be responsible for some 
of the difference between the helium and air experiments, 


Conclusions 


1 Two hypersonic fluid dynamic flow phenomena, on 
inviscid, the other viscous, occur in the flow about a semi. 
infinite flat plate. 

2 Thick leading edges, M > 10 and Re, > 4000 give an 
essentially inviscid flow in which the effect of the leading edge 
is felt for several thousand leading edge thicknesses. 

3 A linear combination of the experimental inviscid flow 


relation and viscous interaction theory correlates with the § Re: 
sharper leading edge data. vind 
tions 
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Research results in the Berkeley low density supersonic 
yind tunnel are described, covering a range of flow condi- 
tions corresponding to flight at altitudes up to eighty 
miles. 


Introduction 


T VERY high altitudes, of the order of 20 miles or more, 
the earth’s atmosphere becomes so rarefied that it no 
longer behaves as a continuous fluid. Proper account must be 
taken of the basic molecular structure of the air in the predic- 
tion of heat transfer and aerodynamic characteristics. The 
general field of rarefied gas dynamics, i.e., the mechanics of a 
gas so rarefied that mean free path effects become important, 
was discussed by Tsien (1).? The terms “slip flow” and “free 
molecule flow” were introduced to characterize phenomena 
associated with moderate and extreme rarefaction respec- 
tively. A slip flow is a flow in which the mean free path is a 
small but not negligible fraction of the boundary layer thick- 
ness and thus corresponds to values of Mach number M and 
Reynolds number Re in the range 0.01 < M/ V/ Re < 0.1, 
since M/ / Re is essentially the ratio of free mean path to 
boundary layer thickness. A free molecule flow is one in which 
the mean free path is large compared to the body size, and 
corresponds to M/Re > 10, i.e., to a ratio of mean free path to 
body dimension of approximately 10 or larger. The slip flow 
range would apply approximately to flight of a one-foot di- 
ameter missile, to take a specific example, at altitudes of the 
order to 20 to 50 miles, while free molecule flow would obtain 
at altitudes above 80 miles. Important applications to in- 
strumentation occur at very much lower altitudes. The gen- 
eral regions of flow, together with the ranges of Mach and 
Reynolds covered by the various experimental programs de- 
sribed subsequently, are indicated in Fig 1. 

Two experimental research groups—one at the Ames 
laboratory of the NACA under the direction of Stalder and 
Goodwin, and one at Berkeley under the initial direction of 
Folsom and Kane—have been carrying on investigations in 
this area for the past nine years with particular emphasis on 
applications to high altitude, supersonic flight. Most of the 
work at Ames has been centered around problems of free 
molecule flow, while at Berkeley, most of the emphasis has 
been on the slip flow range. A fairly complete summary re- 
port, together with an extensive bibliography, of the progress 
of research in the field through the summer of 1952 is con- 
tained in (2). It is the purpose of this report to describe the 
principal results obtained at Berkeley during the period 1952- 
1955 which have application to high altitude aerodynamics. 
These investigations were carried out in the low density super- 
sonic wind tunnel located at Berkeley and described in (2). 
Only those investigations which have not yet been published, 
or which have been published only in the form of project re- 
ports, will be discussed here. Other investigations during 
this period are listed in the bibliography of the present report. 
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Fig. 1 Mach and Reynolds numbers covered by tests 


Of these, special attention should be called to (3) and (4), 
which present data on the surface pressure distribution on 
slender cones at a Mach number of 4, and on the skin friction 
on flat plates at zero angle of attack over a range of Mach 
numbers from subsonic to 3.5. 


Cone Drag 


In an inviscid supersonic flow, the surface pressure on a 
slender cone is constant, providing the Mach number is large 
enough to insure an attached bow shock wave. In a real 
viscous fluid, a boundary layer forms, which is particularly 
thick under the low Reynolds number conditions associated 
with slip flow. This boundary layer effectively increases the 
surface curvature of the cone, displacing the stream lines out- 
ward. It is to be expected that the surface pressure on a cone 
will thus be higher near the vertex and decrease along the 
slant surface, somewhat like the inviscid pressure distribution 
on an ogive. 

In addition to producing a rise in the surface pressure, the 
thick boundary layer which forms on the surface of a cone 
under low density conditions also gives rise to a propor- 
tionately high skin friction. Ipsen (5, 6) has experimentally 
measured the total drag of the fore portion of 15 deg half- 
angle cones over the range 70 < Re < 1600, M ~ 2; 500< Re 
< 7000, M~ 4. Here M and Re are the Mach and Reynolds 
numbers based on the free stream flow and the cone slant 
length. A force-sensitive element mounted inside a 2-in. 
diam cone-cylinder housing was used to measure the total drag 
force on the floating fore portion of the cone, as indicated in 
Fig. 2. This total drag had components due to skin friction 
and also due to the surface pressure rise induced by the 
boundary layer. The technique did not permit separate 
evaluation of the two components. The results are presented 
in Fig. 3, in terms of the increment in the total drag coefficient 
over that corresponding to inviscid flow. The curves indi- 
cated in Fig. 3, are the theoretical predictions obtained by 

Probstein and Elliott (7), which take account of the induced 
effect due to the boundary layer. As will be noted, the agree- 
ment is quite good, except at the lowest Reynolds numbers 
where the effects of slip might be expected to begin to be 
of importance. 
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Fig. 3 Drag coefficients for cones with 15-deg semivertex angle 


Flat Plate Lift 


Normal forces on flat plates at angle of attack have been de- 
termined by Tellep (8). The experiments covered the range 
228 < Re < 2412, 2.42< M < 3.69 at angles of attack up to 8 
deg. Re is the Reynolds number based on the free stream 
flow and on the length of the plate; M is the free stream Mach 
number. The results are presented in Fig. 4, in terms of a 
comparison with the theoretical pressure distribution to be 
expected on the basis of the hypersonic interaction theory of 
Lees and Probstein (9). It will be observed that the agree- 
ment between experiment and theory is good as regards trend, 
but that the theory in general underestimates the actual pres- 
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Fig. 4 Ratio of experimental normal force curve slope to theo- 
retical normal force curve slope vs. Reynolds number 


sure disturbance induced by the boundary layer. This is jp 
accord with experiments at higher pressures conducted by 


Ferri (10). 


As with the cone drag data of Ipsen, there is little indica. 
tion of any slip effects. Both the cone and the flat plate are 
characterized by inviscid constant pressure distributions. To 
first order, therefore, the pressure gradient at the seam of the 
boundary layer is zero. Since the effect of slip to first order 
has been shown (11, 12) to be proportional to this pressure 
gradient, it is not surprising that slip effects are negligible for 
these two geometries. Experiments are now under way with 
a geometry for which the inviscid pressure gradient is not 
zero, with the expectation that slip effects will be more 


apparent. 


Base Pressure on Cone Cylinders 


The pressure on the base of typical cone-cylinder configura- 
tions in supersonic flow has been investigated in (13, 14) for 
the range Re > 40,000, where Re is the Reynolds number based 
on the free stream and the cone-cylinder length. Kavanau 
(15) has now extended these results to lower Reynolds nun- 
bers, 160 < Re < 800, at M ~ 2.10 and 920 < Re < 7400 for 
M ~ 4.0. He finds considerable variation in the pressure over 
the base, the pressure at the center being as much as twice as 
high as the pressure near the edges. His results, in terms of 
the area mean of the base pressure are presented in Fig. 5, to- 
gether with the higher Reynolds number results indicated 
above. No adequate theory now exists for predicting the base 
pressure under these flow conditions. It is therefore difficult 
to determine the relative importance of slip-flow effects upon 


the observed values. 
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Fig. 5 Base pressure coefficient for cone-cylinder configuration 


The thermal characteristics of a series of cones at zero angle 
of attack have been determined in the slip flow range by Drake 
and Maslach (16). The data covered the ranges 87 < Re < 
3270; 2.16 < M < 3.54, where M and Re are the Mach and 
Reynolds number based on free stream flow conditions and 
the cone slant length. The results are presented in Figs. 6 and 
7 in terms of the recovery factor r and the Nusselt number 


Nu defined by 


where 7, T¢ and T> are the stream temperature, the cone 
equilibrium temperature, and the adiabatic stagnation tem- 
perature, respectively, Q is the heat flux to the cone, p, U, and 
Cy free stream density, velocity and specific heat at constant 
pressure, and A is the cone surface area. The increase in the 
recovery factor with decreasing Reynolds number is inter- 
preted as being due to the large dissipation of energy by fric- 
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is in Survey of Boundary Layer Leading Edge Region 
d by -16< 3 <3. 

The thermal recovery factor of any convex body in free 
dics. ar a molecule flow is a rather sensitive function of the Mach num- 
ber over the range 0.5 < M < 3.5. The equilibrium tempera- 
To ne S ° ture of a wire whose diameter is small compared to the mean 
f the . free path and which is stretched across the test section of a 
order wind tunnel thus depends on the Mach number and the local 
ssure stagnation temperature. A “free molecule probe”’ of this sort 
le for ad rT] has been used by Sherman (18) to determine the internal 
with #. | | | structure of a normal shock wave. The probe has been used 
; not : | by Laurmann (19), to explore the flow field in the immediate 
mie || vicinity of a flat plate in supersonic flow. Only a qualitative 

: | | picture can be obtained, since both the Mach number and the 
; | | | | stagnation temperature vary in an interrelated way. How- 
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ever, the general, over-all nature of the flow field can be in- 
ferred. A serious question for the hypersonic boundary layer 
interaction theory is the strength of the shock at the leading 
edge. The existence of a region of slip here could have a very 
great influence far downstream. The results are presented in 
Figs. 8, 9, and 10 in the form of contours through the flow field 
of constant wire temperature. These clearly show the shock 
and boundary layer regions. The possibility of slip at the 
higher Mach number is strongly suggested. These tests were 
carried out using “flat plates” in the form of ribbons only a 
small fraction of a.mean free path thick. This work is now 
being extended to higher Mach numbers. 


Conclusions 


Empirical data are now available for predicting aerody- 
namic and heat transfer characteristics of basic geometrical 
components such as cones, flat plates, spheres, and cone 
cylinders in the slip flow range. These characteristics are, in 
general, quite different from the corresponding values ob- 
tained at higher pressures. Theoretical analyses, however, are 
not yet sufficiently precise to provide an adequate basis for 
determining how much of this difference is ascribable to non- 
continuum phenomena such as slip, and how much to low 
Reynolds number supersonic continuum flow behavior. 


References 


1 Tsien, H. S8., “Superaerodynamics, Mechanics of Rarefied 
Gases,’’ Journal of the Aeronautical Sciences, vol. 13, December 
1946, pp. 643-664. 

2 Schaaf, S. A., “Theoretical Considerations in Rarefied 
Gas Dynamics,’’ chapter IX, Heat Transfer, University of Michi- 
gan Press, 1953, pp. 237-260. “Experimental Methods and Re- 

_ sults in Rarefied Gas Dynamics,’’ chapter X, Heat Transfer, 
University of Michigan Press, 1953, pp. 261-282. 

3 Talbot, L., “Viscosity Corrections to Cone Probes in 
Rarefied Supersonic Flow at a Nominal Mach Number of 4,”’ 
NACA TN 3219, 1954. 

4 Schaaf, 8S. A., and Sherman, F. S., “Skin Friction in Slip 
Flow,’’ Journal of the Aeronautical Sciences, vol. 21, no. 2, Febru- 
ary 1954, pp. 85-90. 

5 Ipsen, D. C., “Cone Drag in a Rarefied Gas Flow,’’ Uni- 
versity of California Institute of Engineering Research Projects 
Report HE-150-114, August 1953. 

6 Ipsen, D. C., “(New Experiments on Cone Drag in a Rare- 
fied Gas Flow,’”’ University of California Institute of Engineering 
Research Projects Report HE-150-128, April 1955. 

7 Probstein and Elliott, “The Transverse Curvature Effect 
in Compressible Axially Symmetrical Boundary-Layer Flow,’’ 
Princeton University, Department of Aeronautical Engineering, 
Report 261, April 1954. 

8 Tellep, D. M., ‘Lift on Flat Plates in Low Density Super- 
sonic Flow,’’ University of California Institute of Engineering 
Research Report HE-150-131, August 1955. 

9 Lees, L., and Probstein, R. F., “Hypersonic Viscous Flow 
Over a Flat Plate,’’ Princeton University Aero. Eng. Lab. Rep. 
no. 195, April 1952. 

10 Ferri, A., “Experimental Results with Airfoils Tested in 
the High Speed Tunnel at Guidonia,’” NACA TM 946, 1940. 

11 Lin, T. C., and Schaaf, S. A., “The Effect of Slip on Flow 
Near a Stagnation Point and in a Boundary Layer,’? NACA TN 
2568, 1951. 

12 Maslen, S. A., ‘Second Approximation to Laminar Com- 
pressible Boundary Layer on Flat Plate in Slip Flow,’’ NACA TN 
2818, 1952. 

13 Bogdonoff, S. A., “A Preliminary Study of Reynolds 
Numbers Effect on Base Pressure at M = 2.95,’’ Journal of the 
Aeronautical Sciences, vol. 19, no. 3, p. 201, March 1952. 

14 Kavanau, L. L., “Results of Some Base Pressure Experi- 
ments at Intermediate Reynolds Numbers with M = 2.84,’ 
Journal of the Aeronautical Sciences, vol. 21, no. 4, pp. 257-260, 
April 1954. 

15 Kavanau, L. L., “Base Pressure Studies in Rarefied Super- 
sonic Gas Flows,’’ University of California Engineering Project 
Report HE-150-125; also a Ph.D. Thesis, 1954. 


250 


16 Drake, R. M., and Maslach, G. J., “Heat Transfer from 
Right Circular Cones to a Rarefied Gas in Supersonic Flow,” 
University of California Engineering Project Report HE-150-9), 
1952. 

17 Drake, R. M., and Backer, G. H., Trans. ASME, vol. 74, 
1952. 

18 Sherman, F. 8., ‘A Low Density Wind Tunnel Study of 
Shock Wave Structure and Relaxation Phenomenon,’’ NACA TN 
3298, July 1955. 

19 Laurmann, J. A., “Experimental Investigation of the Flow 
About the Leading Edge of a Flat Plate,’’ University of California 
Institute of Engineering Research Report HE-150-126, October 
1954. 


Hypersonic Studies of the Leading Edge Effect 
on the Flow Over a Flat Plate 
(Continued from page 246) 


5 Lees, L., “Influence of the Leading Edge Shock Wave on 
the Laminar Boundary Layer at Hypersonic Speeds,’’ GALCIT 
TR no. 1, July 15, 1954. 

6 Shen, S. F., “An Estimate of Viscosity Effect in Hyper- 
sonic Flow Over an Insulated Wedge,’’ Journal of Math. and 
Physics, vol. 31, no. 3, October 1952, pp. 192-205. 

7 Li, Ting-Yi, and Nagamatsu, H., “Shock Wave Effects 
on the Laminar Skin Friction of an Insulated Flat Plate at Hy- 
personic Speeds,’’ Journal of the Aeronautical Sciences, vol. 20, 
no. 5, May 1953, pp. 345-355. 

8 Pai, Shih-I, “A Note on Hypersonic Viscous Flow Over a 
Flat Plate,’’ Reader’s Forum, Journal of the Aeronautical Sci- 
ences, vol. 20, no. 7, pp. 502-503. 

9 Pai, Shih-I, “On Strong Interaction for the Hypersonic 
Boundary Layer on Inclined Wedge,”’ Journal of the Aeronautical 
Sciences, vol. 20, no. 11, November 1953, pp. 796. 

10 Pan, L. J., and Kuo, Y. H., “Compressible Viscous Flow 
Past a Wedge Moving at Hypersonic Speeds,’’ Cornell University, 
Aero. Engng. Dept. 

11 Li, Ting-Yi, and Nagamatsu, H., “Hypersonic Viscous 
Flow on a Non-Insulated Flat Plate, GALCIT no. 25, April 1, 


12 Lees, L., “Hypersonic Flow,’’ IAS Preprint no. 554, June 


13 Becker, J. V., ‘(Results in Recent Hypersonic and Unsteady 
Flow Research at the Langley Aeronautical Laboratory,’’ Jour- 
nal of Applied Physics, vol. 21, no. 7, July 1950, pp. 622-624. 

14 Bertram, M. H., “Viscous and Leading Edge Thickness 
Effects on the Pressures on the Surface of a Flat Plate in Hyper- 
sonic Flows,’’ Reader’s Forum, Journal of the Aeronautical Sci- 
ences, vol. 21, no. 5, June 1954, pp. 430-431. 

15 Bogdonoff, 8. M., and Hammitt, A. G.: “The Princeton 
Helium Hypersonic Tunnel and Preliminary Results Above M = 
11,’’ Princeton University Aero. Engineering Dept., Rep. 260, 
June 1954. (Also published as Wright Air Development Center 
TR 54-124.) 

16 Hammitt, A. G., “The Development of a Helium Hyper- 
sonic Wind Tunnel and the Investigation of Hypersonic Flow 
About Simple Bodies,’’ Princeton University Aero. Engineering 
Dept. Rep. 253, January 1954. 

17 Hammitt, A. G., and Bogdonoff, S. M., “A Study of the 
Flow About Simple Bodies at Mach Numbers From 11 to 15,” 
Wright Air Development Center TR 54-257, October 1954. 
(Also published as Princeton University Report no. 277.) 

18 Hammitt, A. G., Vas, I. E., and Bogdonoff, 8. M., “Lead- 
ing Edge Effects on the Flow Over a Flat Plate at Hypersonic 
Speeds,’’ Princeton University Aero. Engineering Dept. Rep. 
326, September 1955. (Also to be published as WADC TN 55- 
537.) 

19 Hammitt, A. G., “A Preliminary Study of the Details of 
the Flow Field About a Flat Plate at Hypersonic Speeds,’’ Prince- 
ton University Aero. Engineering Dept. Report 327, October 
1955, (Also to be published as WADC TN 55-538.) 

20 Bogdonoff, S. M., and Hammitt, A. G., “Fluid Dynamic 
Effects at Speeds from M = 11 to 15,”’ Journal of the Aeronautical 
Sciences, vol. 23, no. 2, February 1956, pp. 108-117. 


JET PROPULSION 


S 


e 


i 
sur’ 
tail 
the 
Bot 
cut 
was 
tud 
| of 1 
195: 
195 
| 
Su 
or 
C0 | 
— 
los 
de 
| 
de 
| M 
{ 
| 


A semianalytic method is derived for calculating the 
range of one-stage, ballistic, rocket propelled, surface-to- 
surface missiles. The method consists of obtaining cer- 
tain analytic solutions, supposing flight in vacuum, and 
then correcting these solutions for the effect of drag. 
Both the velocity at thrust cutoff and the range following 
cutoff are handled in this manner. An analog computer 
was used to compute the drag corrections and also the alti- 
tude and range at thrust cutoff for a large range of input 
parameters. The results are presented as working formu- 
las and plots suitable for use in preliminary design studies 
of missiles having ranges of from 50 to 500 nautical miles. 


Nomenclature 


area of cross section 
A/W, 
drag coefficient based on area of cross section 
drag coefficient for M = 2 
drag 
function plotted in Fig. 5 
function plotted in Fig. 7 
acceleration due to gravity 
altitude 
7/r = specific impulse 
Mach number 
—(1/W,)-(dW/dt) = burning rate 
radius of the earth 
time 
thrust 
velocity 
loss of cutoff velocity due to drag 
weight 
range after cutoff for vacuum trajectory over spherical 
nonrotating earth 
range loss due to drag after cutoff 
angle shown in Fig. 2 
flight path angle with local horizontal 
air density 
W./Weo 
T/W, 
COS — COs 0, 
= 


QF 


” 


~ 
= 


Subscripts 


orefers to conditions at launching 
¢o refers to conditions at the end of burning 


Introduction 


NE measure of the performance of a ballistic surface-to- 
surface missile is its maximum range with a given pay- 
load. A simple method for calculating this range has been 
devised. The type of missile considered is single stage and 
powered by a rocket motor. 
The method is intended principally for use in preliminary 
design work. It should be of particular use in design com- 
parison studies, economic studies, and investigation of the 
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Ballistic Missile Performance 
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trade-off between payload and range. In such applications 
the effect of changes in the missile parameters is of greater 
interest than the exact magnitude of the range itself. This 
calls for a simplified procedure for calculating range. In 
order to meet this objective, it is necessary to introduce 
several approximations into the work. Sufficient accuracy is 
to be maintained for the intended uses. n 

For the purpose of calculating range the missile may be de- 


scribed by a set of parameters. The basic parameters con- 


sideredare: 1. thrustratio,7; 2. massratio,»; 3. specific 
impulse, J; 4. arearatio,a; 5. drag function, Cp vs. M. 

The ballistic missile derives the major part of its range after 
thrust cutoff. For the case of trajectory in vacuum over 
spherical nonrotating earth, a solution in closed form is ob- 
tained for maximum range after cutoff.‘ 

The range after cutoff is a function of the cutoff velocity 
and altitude. The dominating factor, cutoff velocity, can be 
solved analytically if drag is excluded. 

The drag terms are not susceptible to analytic treatment. 
This difficulty was overcome by using an analog computer. 
The computer also provided cutoff altitude and range. 

The ballistic trajectory following cutoff is considered first, 
in order to establish a certain boundary condition needed for 
programming the burning phase, which is the flight path 
angle at cutoff. 


Assumptions 


The following principal assumptions are made: 


1 The earth is spherical and nonrotating. 

2 Launching and impact are at sea level. 

3 The missile follows a particle trajectory. 

4 Thrust and burning rate are constant. 

5 The missile is launched vertically, and during the burn- 
ing phase d@/dt = const. 

6 The drag corrections are small compared with the 
quantities to which they are applied. 

7 Earth curvature may be ignored during the burning 
phase and also in obtaining the drag correction to range after 
cutoff. 


Thrust and specific impulse, connected through the rela- 
tion J = 7/r, are to be represented by chosen mean values in 
recognition of the fact that rocket thrust varies somewhat with 
altitude. Any secondary propellants must be accounted for 
in r, the total discharge rate, and consequently in the specific 
impulse as defined here. 


Ballistic Trajectory Following Cutoff 


Trajectory and Range 


The vacuum trajectory of an object moving within the 
gravitational field of spherical nonrotating earth, and having 
less than the energy required for escape, is a Keplerian 
ellipse (3). One of the foci of the ellipse is at the center of 


‘ It has been brought to the attention of the authors that simi- 

= solutions have been obtained (1, 2) for the limiting condition 
= 0. 

% Numbers in parentheses indicate References at end of paper. 
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Fig. 1 Vacuum trajectory over spherical earth. Ellipse with 
foci at C and D 


the earth (see Fig. 1). The length of the majorax is L is 
given by 


gk? 
gR — ['/2Veo* + gheo/(1 + heo/R)] 


R + heo 
1 = (1 + heo/R)Vco?/2gR 


Let it be assigned to maximize the range following cutoff, 
X = Rr, for particular values of V., and h... This range is a 
maximum when the length AB (Fig. 1) is maximized. The 
sum of the focal radii AC and CB is given by 


AC + CB = 2L — 2R — ho 


which is fixed for chosen values of V., and h... This being so, 
the length AB is obviously a maximum when the second focus 
C lies on a line joining points A and B. 

The three sides of the triangle ABD are now known for 
the condition of maximum range, and from the cosine law the 
angle max can be solved. The resulting equation can be 
written in the form 


Qmax = 2 sin=! 1 heo/cR 


l-e 


in which * 


Flight Path Angle at Cutoff 


Through use of the property that the focal radii intersecting 
at a point on an ellipse make equal angles with the tangent at 
that point, it can be demonstrated that the angle DAC, 
Fig. 1, is equal to 20,.. The cosine law may be applied to the 
angle DAB to provide a solution for the optimum value of 
6.0; that is, the value that maximizes range for given V, and 
he. The solution yields 


fl—c(2+h./R) 


et 


“app 

For ranges not exceeding 500 nautical miles, 7/2 < 5 deg, s 
that with sufficient accuracy, sin 7/2 = /2 may be used to 
simplify Equation [3]. The maximum range following cutoff, 
for trajectory in vacuum, is then given with good accuracy by 


2c. 
= V1 t+ heo/cR 


Furthermore, since h../R << 1, this term may be neglected 
in Equation [4] if desired. 

Equation [5] may be approximated upon expanding the 
equation and ignoring terms of the order c?, c(heo/R), (heo/R)* 
and higher. Then, since h../R << 1, ¢ is replaced by 
V .0?/2gR, with the result 


sin? Broom, 


Equation [7] may be compared with the expansion of sin? ¢ 
about 2/4, retaining only the first order terms in the expansion, 
and it is seen that 


An alternative form is as follows 
Beoont 


These last equations put the desired cutoff angle at some- 
thing less than 45 deg. In the cases covered by the work 
that follows it would fall between 39 deg and 43 deg. To 
facilitate the burning phase analysis a cutoff angle of 41 deg was 
chosen to handle all cases. 

The remainder of the work will be directed toward the cal- 
culation of a sensible maximum total range. The calculation 
of a true maximum would involve extensive burning phase 
analysis plus the use of variational techniques. Such de- 
tailed analysis is not in line with, and would not contribute 
significantly to, the present objective. 

The range X is to be calculated from Equation [6] regard- 
less of the fact that the burning phase is assumed to be ter- 
minated in all cases at a flight path angle of 41 deg. If the 
total range is near maximum, the combined range-producing 
potential of V.o, heo, and X,. is insensitive to minor changes in 
the burning phase history. An example of this statement can 
be seen in Fig. 2, which is discussed in the next section. 


Burning Phase Trajectory 


Programming 


Since ballistic missiles are being dealt with, one would be 
inclined to assume a gravity turn path for the burning phase. 
However, this would complicate the analysis. For every set 
of missile parameters several trials would have to be made, 
each time varying the initial launch angle, in order to arrive at 
the prescribed boundary condition on flight path angle at cut- 
off. It was necessary to avoid this. 

Instead, the missile is considered to be launched vertically, 
and caused to turn at a constant rate during the burning 
phase. The effect of this programming was checked against 
a maximum range gravity turn trajectory. The result, 
shown in Fig. 2, was obtained with an analog computer. 
Since the difference in range could not even be detected, the 
chosen programming is evidently a good compromise in so far 
as it simplifies the analysis. 

It is intended that the total range to be calculated in the 
present method be understood to be essentially the maximum 
range of a missile following a zero lift, or ballistic, trajectory 
throughout flight. 
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CONSTANT RATE OF TURN 
DURING POWERED FLIGHT 


50 


GRAVITY TURN(ZERO LIFT) 
DURING POWERED FLIGHT 


50 100 150 
RANGE - N. Mi. 
Fig. 2 Comparison of trajectories having different programming 
during powered flight 


Cutoff Velocity 


The cutoff velocity is obtained upon integrating the equa- 
tion of motion over the time of burning. The equation of 
motion in the direction of flight is 


T — D— Wsin = (W/g)(dV/dt)......... [10] 


Under the assumptions that thrust, burning rate, and rate of 
turn are constant, the integration furnishes a formula for the 


cutoff velocity. 


gT W. COS — COS feo 
Veo =. log — — — —— gla — 
Weo — Bo 0 W.(1 t) 


The burning time is obtained from 


ll 


Computer Operations 


The unsolved term in the equation for cutoff velocity must 
be solved. The cutoff range and altitude must be calculated. 
Finally, a correction for air resistance must be applied to the 
range after cutoff. 

This was all done with the aid of an analog computer. The 
computer was wired to solve the complete trajectory of a bal- 
listic missile in the presence of the earth’s atmosphere. The 
curvature of the earth was neglected. The atmospheric 
functions, speed of sound and air density, were approximated 
as analytic functions of altitude, so that they could be gener- 
ated using ordinary machine components. 

The air density was represented by two third-degree poly- 
nomials in altitude. These were obtained as a least-squares 
fit to data for the NACA Standard Atmosphere. The first 
polynomial supplied the atmospheric density ratio from sea 
level to 100,000 ft. At 100,000 ft the second polynomial was 
switched in, supplying the density ratio between 100,000 ft 
and 200,000 ft. Above 200,000 ft the density was assumed 
zero. 

The speed of sound was approximated as decreasing lin- 
early with altitude, from a value of 1116 ft/sec at sea level 
to a value of 972 ft/sec at 35,000 ft above which it was as- 
sumed to remain constant. 


Drag Function 


A drag curve made up of analytic elements, as shown in 
Fig. 3, wasassumed. This curve is believed to be a reasonable 
approximation to the typical drag curve of a cylindrical body 
with pointed nose and stabilizing fins. 

Despite the arbitrarily assumed burning phase program- 
ming, neither lift forces nor drag due to lift were considered 
in the analog computer operations. Reasons for this were 
covered in the discussion of Fig. 2. The drag was computed 


through the equation 


1956 


Cy = .09+.5/M 
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tet tes 


TIME - SECONDS 

Fig. 4 Hand calculation (left) and REAC calculation (right) 

for flat earth performance of a hypothetical one-stage missile 

having the characteristics: J = 235 sec; v = 2.71; + = 2.12; 
a = 0.00082 ft.?/Ib.; Cp vs. M from Fig. 3. 


As a check on the operations of the computer, Fig. 4 
presents for comparison the results of an analog computer 
calculation and a hand calculation made using the same in- 
put parameters. 


Parameter Bounds 


Using the analog computer, a large number of runs were 
made covering the burning phase. The parameters that 
were varied, and their bounds, are 


15<7r1< 4.0 

Specific impulse ......... 200 <I< = sec 
Arearatio ..........7... 0.0005 < a < 0.0015 ft?/lb 


Burning Phase Drag 


The unsolved term, AVp, in Equation [12], which is the 
decrease of cutoff velocity due to drag, is given by 


The integral in Equation [15] was computed over a wide 
range of parameters. Function F, defined by 


pV? 
F = — dt... 16 
af 1-—rt 


in which Cp* is the drag coefficient for a Mach number of 2.0, 
can be plotted as a function of thrust ratio 7 and mass ratio » 
alone, since the calculated results were found to be insensitive 
to changes in the remaining parameters. The C'p* behaves as 
a mean value of drag coefficient, which was the desired effect. 

The plot for function F is shown in Fig. 5. This and the re- 
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Fig. 5 Drag function for correcting velocity at cutoff 


maining figures provide working plots for use in the present 
method. The value of F obtained from Fig. 5 is inserted in 
the following equation to solve the final term in Equation [12]. 


AV p (ft/sec) = aCp*F.............. [17] 


Cutoff Altitude and Range 


The burning phase calculations made with the analog com- 
puter also provided cutoff altitude and range, shown plotted 
in Figs. 6 and 7 as a function of cutoff velocity for various 


Z 


x 
\ 


CT 

*o 2 4 6 8 0 


Ye 1000 FT. /SEC. 
Fig. 6 Altitude at cutoff 


values of parameter r. The plotted data were observed to 
be insensitive to changes of parameters » and a. A further 
dependence of h,, and X.. on specific impulse is not reported in 
the plots, because of its rather small consequence. An ap- 
proximate correction based on the original plotted data may 
be applied to the values of h.. and X.. read from Figs. 6 and 7, 
through the formula Ah.o/heo = AXco/Xeo = 0.8AI/I., in 
which J, = 220 sec. 


Drag After Cutoff 


It remains to correct the range following cutoff for the 
effect of air resistance. This was also done using the analog 
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Fig. 7 Range at cutoff 


computer. A method somewhat similar to that handling 
burning phase drag wasemployed. The range was computed, 
both in atmosphere and in vacuum, for a number of combi- 
nations of cutoff velocity and altitude, and three values of area- 
to-cutoff weight ratio, va. It was observed that the difference 
between these range computations, which is the loss of range 
due to air resistance after cutoff, is proportional to va. A further 
factor, dvag coefficient just after cutoff, Cp.., is introduced. 
The reason for this particular choice of drag coefficient is that 
the drag losses immediately following cutoff account for the 
greater part of the range reductior. 

Function G was calculated through an inversion of the equa- 
tion 


and then plotted as shown in Fig.8. The range correction can 


Ned 1OODFT. 


rd 


° 


8 


200 


0 2 4 8 
Veo - 1000 FT./ SEC. 


Fig. 8 Drag function for correcting range after cutoff 


now be reconstructed upon inserting into Equation [18] the 
value of G obtained from Fig. 8. In determining Cp* and 
Cdeo, account may be taken of any change in drag due to the 
application or termination of thrust. 


Application and Evaluation 
Parameter Flexibility 


The range of specific impulse values admissible under this 
method is somewhat wider than the bounds that were shown 
above. The reason is that the primary effect of specific 
impulse is handled explicitly in the propulsion term of 
Equation [12]; secondary effects are unimportant. 

There is also reason to believe that the procedure for 
handling drag permits the use of drag functions differing in 
detail from the curve shown in Fig. 3. This is found in the 
results of the burning phase calculations using various values 
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of the area-to-weight ratio,a. Changes in this area parameter 
cause similar changes in the drag. As the over-all drag level is 
changed, it causes the burning-phase histories of velocity, 
altitude, and thus dynamic pressure to be altered. The effect 
igsimilar to distorting the drag curve itself. It was observed, 
however, that the function F, needed to reconstruct the drag 
correction to cutoff velocity, remained insensitive to this 
eect. This leads to the conclusion that moderate changes in 
the shape of the drag curve can be tolerated. 


Range Limits 

At decreasing ranges the effect of drag rises in relative im- 
portance. If the major part of the drag loss occurs at tran- 
sonic or especially subsonic speeds, the method as presently 
constituted cannot accurately predict the drag effects. For 
this reason caution should be used in applying this method to 
the calculation of ranges less than 50 miles. 

For ranges exceeding 500 miles, the optimum cutoff angle 
begins to swing too far below 41 deg for the present results to 
be valid. 


Accuracy 

The original data have been examined to provide an idea of 
the accuracy of the method. Within the parameter bounds 
previously tabulated it is estimated that the errors contained 
in the values read from Figs. 5 to 8 will seldom exceed the fol- 
lowing: 

Function F, Fig. 5 

heo, Fig. 6 

Function G, Fig. 8 

It is estimated that use of Cp* as directed will not usually 
result in an error in excess of +10 per cent for drag curves 
other than that assumed here. 

The total range is affected directly by errors in X.,. and AX p. 
The effect on range of errors in V.,. and h,. can be determined 
by writing the first order approximation for Equation [6] 


+(0.25 X 10° units error 

+5Q% error (when corrected for J) 
+59% error (when corrected for J) 
+ 10% error 


2 
g 
from which 
AX = AVco/g)dVeo + dheo.... [20] 


The error in h,, thus produces an equal error in range, while 
the range error produced by error in V,, is 2V../g times the 
error in Vo. 


Calculation Procedure 


In making use of the present method the parameters r, », 
I, a, and the function Cp vs. M must be known. The steps 


in calculating range are: 
1 Calculation of velocity at thrust cutoff: 


I 
V.. (ft/sec) = 32.27 log — 28.4 - — aCp*F..... [21] 
v 


2 Determination of altitude and range at cutoff. The 
values read from Figs. 6 and 7 can be rendered more ac- 
curate by correcting them for specific impulse, as previously 
noted. 

3 Summation of the components of the range: 


total range = X.. + V1 + heo/cR — vaCnecG. . [22] 


in which 


Numerical Example 


The hypothetical missile for which the parameters are given 
in Fig. 4 may be used to illustrate the application of the 
method. From Fig. 3, Cp* = 0.34; and from Fig. 5, F = 
2.75 X 108. The cutoff velocity is now solved using Equation 
[21]: 

Ves = 4780 ft/sec 


Using Figs. 6 and 7, and making corrections for specific 
impulse: 
heo = 99,000 ft 
Xeo = 72,000 ft 


The Mach number at cutoff is 4780/972 = 4.9. Figs. 3 
and 8 are entered to find Cp.. = 0.19; G = 28 X 10° 
Before calculating the range, the function c must be solved: 


2cR = 714,000 ft; c = 0.0170 
The total range is found from Equation [22]: 
Total range = 882,000 ft or 145 n.mi. 


The optimum value of flight path angle at cutoff is, from 


Equation [8], 41.0 deg. 
The hand calculations for the performance of this missile 


yielded the following results: 


Veo = 4800 ft/sec 
heo = 100,800 ft 
Xoo = 69,400 ft 
Total range = 146 n.mi. (corrected for earth curvature) 


The maximum expected errors in the example worked out 
above are as follows: Vo, 70 ft/sec; Aco, 5000 ft; Xo, 3600 ft; 
AXp, 1200 ft. The range error produced by the 70 ft/sec 
error in V,. is 20,800 ft, which is here by far the most im- 
portant error. 

Under the conservative assumption that the errors in h.. 
and X,. are mutually consistent, the maximum expected total 
range error for the example, taken as the rms sum of 8600, 
1200, and 20,800 ft, comes to 22,500 ft or 3.7 n.mi. Had an 
allowance for 10 per cent error in C’p* been made to provide for 
different drag characteristics, the error in V.. would be 


J (70)? + (77)? = 104 ft/sec, and the total expected range 
error would be increased to 5.3 n.mi. 
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Effect of Chemical Reactions in the Boundary Layer 
on Convective Heat Transfer 


DAVID ALTMAN? and HENRY WISE*® 


Jet Propulsion Laboratory, California Institute of Technology, Pasadena, Calif. 


An analysis is made of the relative change in heat- 
transfer rate resulting from chemical reactions in the 
boundary layer for either turbulent or laminar flow. 
Consideration is given to reactions both in the gas phase 
and on the wall. It is shown that, for gas phase reactions, 
if D = x (where D is molecular diffusivity and x thermal 
diffusivity) or if reaction occurs in which both D and x are 
much less than ¢ (the turbulent-exchange coefficient), the 
relative change in heat flux is given merely in terms of 
thermodynamic properties alone. For this case, the use 
of a “‘nondissociated’’ temperature representing the total 
enthalpy of the system is found to be useful. 


Nomenclature 


isobaric specific heat 

isobaric molar heat capacity 

molecular diffusivity 

molar reaction rate per unit area of surface 

molecular weight 

molar concentration 

mole fraction 

Prandtl number 

heat flux per unit area 

heat of reaction per mole of diffusing species (— AH of 
reaction ) 

molar reaction rate per unit volume 

universal gas constant 

Schmidt number 

temperature 

coordinate of heat flow 

dimensionless distance parameter = (y/v)+\/7/p 

rate of molecular collisons per unit surface 

accommodation coefficient for surface reaction (cf. Eq. 
14) 

dimensionless constant (cf. Eq. 15) 

turbulent-exchange coefficient 

thermal diffusivity 

thermal conductivity 

viscosity 

kinematic viscosity 

density 

shear stress 

wall (y = 0) 

edge of laminar sublayer (y* = 5) 

outer edge of intermediate layer (y* = 30) 

property including turbulent contribution 

diffusing species 

log mean of property in film 

condition without chemical reaction 

average property in boundary layer 

fictitious free-stream property calculated in the absence 
of chemical dissociation 


q 
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I Introduction 


N THE combustion processes which yield a high concen- 
tration of chemically active species such as free radicals 
and atoms, the heat transfer between the free stream and the 
wall may be significantly altered by chemical reaction within 
the system. A simplified analysis is given of the relative 
change in heat flux from a hot moving gas stream to a cooled 
surface as a result of chemical reactions in the boundary layer 
and on the wall. 

The assumptions and approximations incorporated in this 
treatment are as follows: (a) The gradients of the quantities 
considered along the flow line parallel to the wall are negligible 
with respect to gradients normal to the wall (y direction). (b) 
The coefficients of turbulent exchange for viscosity, diffusiv- 
ity, and conductivity (eddy quantities) are equal. (c) The 
boundary layer thickness does not change significantly as a 
result of heat release; i.e., the Prandtl and Schmidt numbers 
remain essentially constant (1) and (2).4 (d) Steady state 
has been established, and pressure is uniform within the sys- 
tem. (e) The mole fraction of reactive species is low if more 
than one such species is considered. (f) The cold-wall tem- 
perature 7» is fixed at a relatively low value by coolant flow. 
(g) Thermal diffusion of species is neglected. For species of 
equivalent molecular weight, this effect is known to be small. 
In the case of species with molecular weight much different 
from that of the mixture (such as H atoms) the inclusion of 
thermal diffusion would still act only as a correction in the 
present case and is therefore neglected to the degree of ap- 
proximation employed in this paper. (h) Radiation effects are 
not considered. 


II Derivation of Equations 


Using von Kaérman’s treatment of the analogy between fluid 
friction and heat transfer (3), it is assumed that the boundary 
layer contains a laminar sublayer of thickness y* = 5 and an 
intermediate layer (buffer layer) in the interval 5 < y* < 30. 
In the sublayer, the transport properties are given solely by 
the molecular quantities » = u/p, x = d/cp, and D. In the 
intermediate region, these quantities increase linearly with 
distance as a result of turbulence. The coefficient of turbulent 
exchange ¢ (3), also called exchange coefficient (4), represents 
the turbulent contribution to the transport quantities v, x, and 
D (Fig. 1). If the effective transport parameters are denoted 
by ke = x + eand D, = D + eto represent appropriate values 
in boundary layer, the conservation equations for heat and 
mass may be written as follows 


4 
dy dy 
d( Dan aN; 
dy 1-—WN i dy 
where Q is the heat released in the recombination of 1 mole of 


the labile species (—AH of reaction), ® is the reaction rate 
expressed in moles per unit volume per sec, p is the density, 
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Fig. 1 Variation of exchange coefficient and effective diffusivi- 
ties in the boundary layer 


is the molar density, and N; is the mole fraction of free radical. 
Elimination of ® between Equations [1] and [2] results in 


( aT +Q ( Dan aN; 0 

dy _ dy dy \l — N; dy 
Since the heat of reaction Q is essentially constant, a first inte- 
gral of Equation [3] between the limits y = 0 and y = y, gives 
directly the heat-transfer rate per unit area of wall surface 
dT dN; 


dy 


Division by cpxe and integration between y = 0 and y = y» 
corresponding, respectively, to JT = T) and T = T, and to 
N;= No and N; = Nie (ef. Fig. 1) yields 


d Ni Dd, 
af, 


dN... 
Nio (1 — N;)Cke 5] 


where cp/n = C is the average molar isobaric heat capacity of 
the gas mixture. If no chemical recombination occurs, Equa- 
tion [5] reduces to the following equation, in which the primes 
denote appropriate values subject to the initial assumption of 
constant and 


Since, as was stated in Section I, chemical heat sources have 
only a small effect on the Prandtl number (2), the boundary 
layer thickness yo’ may be taken to be the same as ya. Further- 
more, since the temperature limits over the boundary layer 
are the same, the integrals on the left-hand side of Equations 
[5] and [5’] may be equated within a reasonable approxima- 
tion. The ratio of heat flux with and without recombination 
for identical boundary conditions in temperature is therefore 


q Q Dt+e 
q’ — Ty JS Nio (1 — N,)C(x + €) 


If more than one recombining species exists and if their com- 
bined mole fraction is low, Equation [6] can be easily general- 
ized by summation of all the species 7 which react as follows 


dN;.... [6] 


1 j Nia D; +e 
| 
+ Nio (1 — + €) 


where i = 1, 2,3,...,7. Dz is assumed to be an appropriate 
binary diffusion coefficient in the mixture.' In the three cases 
of Equation [6] that follow, it is assumed that this generaliza- 
tion can be made. 
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A_ Reaction in Boundary Layer Where D. = x. 


In this first case a simplified expression for the relative 
change in heat flux is obtained if either of the following two 
conditions is satisfied: 

1 D = « with reaction occurring anywhere in the entire 
boundary layer. 

2 Reaction occurring only in the hot portion of the inter- 
mediate layer where D and «x are much less than e. 

In either case, Equation [6] can be integrated directly to 
yield 


q Q(Niz — Noo) 


where N,, is the log-mean mole fraction of the remaining 
species. For the recombination of various labile species of low 
mole fraction 

2Qi(Ni2 — Nio)/(1 — Nim) 

C(T: — To) 
where C(= ¢M) is the average molar heat capacity of the gas 
mixture in the temperature range of recombination.® It is in- 
teresting that, in this case, the increase in heat transfer is no 
longer dependent on the transport properties but only on the 
ratio of recombination heat to sensible heat. Although Ni2 
and T. do not refer exactly to free-stream properties, such an 
identification can be made since, if the conditions of this first 
case hold, D. = x. beyond ys, and integration of Equation [6] 
can be made into the free stream. 

If the recombination rate is fast, chemical equilibrium is 
closely approached in the boundary layer. Since the recom- 
bination heats for most of the free radicals are usually high 
(some 50 keal/mole and higher), it follows that recombination 
occurs largely in the high-temperature region of the inter- 
mediate layer where « >> x. This condition is probably ap- 
proached for many reactions at high flame temperatures. If 
T> is 1000 K or so less than 7», it is generally a good approxi- 
mation to set Ni equal to zero if chemical equilibrium is 
maintained, whereas Niz and T, may be identified with equi- 
librium flame conditions. However, it is important to note 
that, regardless of whether or not chemical equilibrium is fol- 
lowed in the diffusion process, Equation [7] or [7’] is applica- 
ble. Hence the increase in heat transfer is independent of the 
manner or location of heat release except for the minor change 
of the molar heat capacity C, which is averaged over the tem- 
perature range of chemical recombination. 

Another simplification of Equation [7’] proves of practical 
utility for the case where N jo is negligible. If T.° denotes the 
calculated flame temperature for a combustion reaction where 
the products are assumed to be only those stable at ambient 
temperatures (i.e., H:O, CO, Ne, etc.), then the 
energy-balance equation yields | 


where C° is the average molar heat capacity for the stable 
products between the temperatures JT, and T°. At the high 
temperatures where chemical dissociation is important, 
equipartition of energy is closely obeyed, and C° and C differ 
only slightly, mainly as a result of a small composition dif- 
ference and anharmonicity effects of the internal-energy 


states. Under the conditions that C° = C, that Nin <1 [so 
that 1/(1 — Nim) ~ 1 + Nim], and that the Ni values are 


5 The accurate treatment for multicomponent parallel diffusion 
through a stagnant layer is quite intractable. Sherwood (5) has 
given a solution for two parallel diffusing species which was ob- 
tained by Gilliland. Wilke (6) has indicated an approximate 
method of obtaining D; in multicomponent mixture which may 
be suitable for this case. 

6For recombination at or near chemical equilibrium, this 
average corresponds to a temperature close to 7, because the de- 
gree of dissociation rises exponentially with temperature. 
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zero, Equation [8] may be incorporated into Equation [7’] to 
yield 
q TS ZQiNi2Nim 


q’ To 

It may be noted that the second term in Equation [9] is small 
with respect to the first and may be neglected for approximate 
purposes. In using Equation [9], 72° is calculated very 
simply by writing the energy balance between combustibles 
and stable products for the single net stoichiometric equation 
neglecting dissociation. 


B Reaction Restricted to Surface 


In the second case there are several types of reactions which 
are catalyzed by surfaces. Among such reactions are atomic 
and free-radical recombination, the endothermic dissociation 
of certain species such as ammonia, and certain chain reac- 
tions. The particular case of a reaction occurring solely on the 
surface can easily be treated in this scheme. Equations [1] 
through [6] are still valid in this second case, and if the re- 
striction D = x quoted in the first case is satisfied, then 
Equation [7] can be employed. However, it is incorrect to 
employ Equation [7] if e > D and « and if D = x, since the 
integral in Equation [6] extends over the entire boundary 
layer. In order to treat the general case of D # « (as occurs 
with H atoms, for instance), it is convenient to solve Equa- 
tions [1] and [2] independently with & = 0 since it is assumed 
that no reaction occurs in the bulk phase. In this treatment 
the heat liberated at the wall by surface reaction is entirely 
absorbed there by an appropriate coolant flow which holds the 
temperature T, fixed. This condition implies that the occur- 
rence of surface reaction has a negligible effect on the tempera- 
ture contour in the boundary layer and that therefore the con- 
vective flow of heat q’ is the same with and without surface 
reaction. The second integral of Equation [1] can be ex- 
pressed as 
dy 
0 cp(x + e) 

Denoting ¢ and j as the mean values in the film, the resulting 
integral can be solved according to the assumption given by 
von Kérmén (3). The dimensionless distance y* = (y/v) 
(v. t/p) is employed in which it is assumed that « = 0 for y* 
< 5and ¢/v = (y*/5) — 1 for 5 < y* < 30. Since »/x = Pr, 
the result of the integration yields 

dy 5[Pr + In(1 + 5Pr)] 

A corresponding treatment can be made for solution of the 
conservation-of-mass equation (Equation [2]) with ® = 0. 
Let F; be the rate at which the diffusing species reacts per unit 
wall surface; then it is easily shown that 

De (*) _ (Nix Nuo\n/Nm)V 7/6 
*"1-Ni\dy 5[Se + In(1 + 5Sc)] 


The relative heat flux is therefore 


_ _ 1+ — Nio) E + In(1 + 
q NwC(T: — To) + In (1 + 58c) 


| 


Clearly, if D = «(Sc = Pr), Equation [13] reduces directly to 
Equation [7] as just stated. However, the surface recombina- 
tion of H atoms, for which the-Pr > Sc, causes a magnification 
of the heat release effect since the diffusion of H atoms carry- 
ing recombination heat exceeds that of the remaining 
molecules.carrying thermal heat. 

The remaining problem is the evaluation of Nw» in terms of 
the accommodation coefficient a, which is the fraction of 
molecular collisons with the wall leading to reaction. If Z; 
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is the total number of surface collisons per unit time made by 
the diffusing species, the following equation’ for F; defines q 


F; = aZ; = awN; 


Equations [12] and [14] may now be used to evaluate Vu in 
terms of known quantities. Defining 6 as 


nV 7/6 
5no [Sc + In (1 + sso 


B= 


RT, 
27M; 


the solution for Nw is 


and Equation [13] reduces to 


1+ E +In(1 + 
[1 + (8/a)|NmC(T: — To) | Se + In(1 + 5Sc) 


For small values of 8/a, N io > 0, whereas for B/a > 1, Nin > 
Ny. As a typical example, consider a flame at 3000 K with a 
velocity of 100 fps at y2 in a tube with the inner-wall tem- 
perature at 1000 K. With Sc = 1 and M; = 20, the value of 
8 is about 3.2 X 10-5, and any value of a > 10-3 leads to 
Nw 0; i.e., complete reaction at the surface. 


C Reaction Restricted to Laminar Layer 


In the third case, if reaction occurs only in the laminar 
layer where « = 0, Equation [6] can be integrated between 
the limits N; = Nw and N; = Nx provided D/x = Pr/Scis 
independent of N; and T as predicted by kinetic theory for 
rigid spherical molecules. Under these restrictions, the result 
is 


QPr 
q’ CSc Nn T, —T> 


For laminar flow where My and T; are close to free stream 
properties, Equation [18] is useful for direct computation. 
For the general case including convective flow, an accurate 
solution can be obtained for N; and 7’, by the separate solu- 
tions of Equations [1] and [2] in the regions y* < 5 and 5<y* < 
30. In the sublayer where y* < 5, Equations [1] and [2] must 
be solved using an appropriate expression for R = f(Ni, T). 
The unknown values 7; and Nx are then determined by 
matching with the solutions obtained from Equations [1] and 
[2] in the region 5 < y* < 30 where R = 0. However, it 
may be noted that, if D = x, the solution reduces to Equation 
[7] even though reaction occurs only in the laminar layer. 
Such a situation may occur in hot flames containing easily 
dissociated species such as Br or I atoms in low concentration. 


III Specific Application 


It may be of interest to examine the application of the equa- 
tions derived to a specific flame system transferring heat to a 
cold wall. In most cases dealing with dissociated gases in 
high-temperature flames, the rates of chemical recombination 
are fast relative to the rate of transport through the boundary 
layer; thus, to a good approximation, the mole fraction of 
labile species at the wall may be considered negligible if the 
wall temperature is sufficiently low. The data represented in 
Tables 1 and 2 refer to the combustion products resulting from 
the reaction of liquid ammonia with gaseous oxygen at 


(Continued on page 269) 


7 Strictly s ing, m, Nio, and 7 should be evaluated one 
mean free path from the wall. However, for most practical cases, 
the magnitude of one mean free path is much less t tj conse- 
quently the simplification employed is valid. 
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(16 This paper deals with two limiting cases of laminar heat 
transfer over blunt-nosed bodies at hypersonic flight 
speeds, or high stagnation temperatures: (a) thermo- 

) dynamic equilibrium, in which the chemical reaction 

] rates are regarded as “‘very fast’? compared to the rates of 
diffusion across streamlines; (b) diffusion as rate-govern- 

[17] ing, in which the volume recombination rates within the 

boundary layer are “‘very slow’? compared to diffusion 

h 6 across streamlines. In either case the gas density near the 

my surface of a blunt-nosed body is much higher than the 

a density just outside the boundary layer, and the velocity 

i: and stagnation enthalpy profiles are much less sensitive 
to pressure gradient than in the more familiar case of 
moderate temperature differences. In fact, in case (a), 
the nondimensionalized enthalpy gradient at the surface 
is represented very accurately by the “‘classical’’ zero pres- 

dar sure gradient value, and the surface heat-transfer rate 

en distribution is obtained directly in terms of the surface 

r 1S pressure distribution. In order to illustrate the method, 

for this solution is applied to the special cases of an unyawed 

ult hemisphere and an unyawed, blunt cone capped by a 
spherical segment. 

In the opposite limiting case where diffusion is rate- 

8) controlling the diffusion equation for each species is re- 
duced to the same form as the low-speed energy equation, 

m except that the Prandtl number is replaced by the Schmidt 

n. number. The simplifications introduced in case (a) are 

te also applicable here, and the expression for surface heat 

u- transfer rate is similar; the maximum value of the ratio 

< between the rate of heat transfer by diffusion alone and by 

st heat conduction alone in the case of thermodynamic 

). equilibrium is given by: (Prandtl no./Schmidt no.)*/*. 

ry When the diffusion coefficient is estimated by taking a 

d reasonable value of atom-molecule collision cross section 

it this ratio is 1.30. Additional theoretical and (especially) 

n experimental studies are clearly required before these 

r. simple results are accepted. 
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Nomenclature 


The subscript w refers to gas properties evaluated at the local 
surface temperature, while the subscript e denotes quantities 
at the outer edge of the boundary layer, and the subscript © refers 
to ambient properties just ahead of the bow shock wave. The 
subscripts 7 and j represent the ith and jth chemical species. The 
subscript o generally denotes quantities at the forward stagnation 
point. A prime denotes differentiation with respect to 7. 


a = speed of sound 

C; = concentration by weight of the ith species 

C., Cp = specific heats at constant volume and constant pres- 
sure, respectively, per unit mass (C, = =C%C>,) 

C(n) = pp/pete Chapman-Rubesin function 

D = diameter 

Di; = coefficient of diffusion of species ‘‘i’’ into species ‘‘j’’ 
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U/Ue 

static enthalpy 

standard enthalpy of pene : 

total enthalpy per unit mass 

coefficient of thermal conductivity 

Mach number u/a 

molecular weight 

Nusselt number, hD/k, where G, = h/Cy, (hs, — hw) 

static pressure 

Lewis number, 

Schmidt number 

stagnation pressure behind normal shock wave 

Prandtl number, Cpu/k 

local heat transfer rate, per unit area, per unit time 

radius of cross section of body of revolution 

distance between centers of two molecules or atoms 

nose radius 

universal gas constant 

distance along body surface, measured from forward 
stagnation point 

transformed coordinate along body surface (Eq. [7]) 

absolute temperature 

rT = kT'/e, where k is the Boltzmann const = 1.380 X 10~"6 
erg/°K 

components of velocity parallel and normal to surface 

net rate of production of ith species, gm/cm*/sec 

distance along body axis, measured from forward stag- 
nation point 

distance normal to surface 

nondimensional concentration 

velocity gradient parameter, (2%/u.) (du./d3) 

(28/M.) (dM./d8) 

ratio of specific heats, C,/C, 

boundary layer thickness 

boundary layer displacement thickness, 


So” (1 — u/we] dy 
maximum attraction energy, erg 
nondimensional coordinate normal to body surface 
(Eq. [7]) 
= angle between flight direction and radius vector from 
center of curvature of nose 
angle of inclination of local surface element with 
respect to flight direction 
cone half-angle 
absolute viscosity 
kinematic viscosity, 
gas density 
zero-energy distance, or collision diameter 
¢(r) interaction potential for intermolecular force field 
y stream function 
™* (T*) and ™)* (7'*), collision integrals 
We ReT 


-. 


>> 


PR 


i] 


i] 


§ 


= 
I 


1 Introduction 


-T HYPERSONIC flight speeds all bodies must be blunt- 
nosed to some extent in order to reduce the heat transfer 
rates to manageable proportions and to allow for internal heat 
conduction. Since the maximum local heat transfer rate in 
most cases will occur in the vicinity of the forward stag- 
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nation point, the fluid-mechanical problems associated with 
this region are of considerable practical, as well as theo- 
retical, interest. Because of the strong falling pressure gra- 
dient and large heat transfer to the surface, the boundary 
layer in this region is expected to be laminar and therefore 
amenable to analytic treatment. 

At hypersonic flight speeds the detached bow shock wave 
just ahead of a blunt-nosed body converts most of the kinetic 
energy associated with the flight velocity into thermal and 
chemical energy. At the “satellite escape speed” of 26,400 
ft/sec, for example, this energy amounts to about 14,000 
Btu per Ib of air, which is sufficient to dissociate almost all 
of the nitrogen and oxygen molecules. The mechanism of 
heat transfer in the boundary layer under these conditions 
may differ markedly from the purely molecular heat conduc- 
tion process of conventional fluid mechanics. For this reason 
it is desirable to examine certain limiting cases in order to be- 
come familiar with the relative importance of the various 
possible physical processes, and to provide a basis for com- 
parison with experiment. Two important limiting cases are: 

1 Complete thermodynamic equilibrium, in which the 
gas properties and concentrations of atoms, molecules, and 
ions are identical with their equilibrium values appropriate 
to the local temperature at every point of the flow. Here the 
chemical and electronic reaction rates are regarded as ‘‘very 
fast” in comparison with the rates of convection along 
streamlines, or the rates of diffusion across streamlines. 

2 Diffusion as rate-governing for heat transfer. In this 
opposite limiting case the volume recombination rates within 
the boundary layer are “very slow” compared to diffusion 
across streamlines, because of the relatively low temperatures 
near the surface. Heat transfer from the boundary layer to 
the body surface is accomplished partly by ordinary heat 
conduction, and partly by the heat released by catalytic re- 
combination of those atoms that manage to penetrate the 
screen of molecules diffusing away from the surface.” 

The actual physical situation is expected to lie somewhere 
between these two limiting cases, at least at flight altitudes 
low enough so that the gas can still be regarded as a “con- 
tinuum,” and heat transfer to the body surface is still of 
major importance. 

Even in these idealized limiting cases the laminar boundary 
layer problem is not simple, unless certain special approxi- 
mations are introduced. Cohen and Reshotko (1)? extended 
Thwaites’ approximate momentum-integral method based on 
the “similar solutions’ to compressible flows with constant 
specific heats. Also, S. Lal (2) developed straight-forward 
techniques for applying the momentum and energy integral 
methods to such flows. In a recent paper Hayes‘ shows 
how the Stewartson-Illingworth transformation, which re- 
duces the problem to an equivalent low-speed flow, is general- 
ized (at least approximately) to the case of an arbitrary im- 
perfect gas. He also discusses the various methods and ap- 
proximations available for solving the reduced problem. By 
following Hayes’ approach, these methods could now be 
modified to include the case of temperature-dependent specific 
heats. 

Fortunately, at hypersonic flight speeds, as we shall see in 
Section 2.2, the problem is simplified considerably by recog- 
nizing that the gas density near the surface of a blunt-nosed 
body is much higher than the density just outside the bound- 
ary layer. The ratio p../p, is large, i. e. 10-20, not only because 
T./T.< 1, but also because the molecular weight of the gas 
near the surface is roughly twice the average molecular weight 
just outside the boundary layer. Since the pressure gradient 
in the direction parallel to the surface depends on the ‘‘ex- 
ternal” density, p,, and is uniform across the boundary layer, 


? The present paper does not consider the problem of radiative 
heat transfer from the hot gases outside the boundary layer to 
the surface. 

* Numbers in parentheses indicate References at end of paper. 

4 Private communication. 
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the velocity profile near the surface is much less sensitive to 
the pressure gradient than in the more familiar case of moder- 
ate temperature difference across the.boundary layer. The 
stagnation enthalpy distribution is even less sensitive. Ip 
fact, for pw/p. > 1, we will show that an excellent first ap- 
proximation is obtained by ignoring the direct effect of pres- 
sure gradient entirely in comparison with the influence of the 
local pressure on the gas density.5 The velocity profile is then 
identical with the Blasius distribution in properly transformed 
coordinates, and the heat transfer rate distribution over the 
surface of a blunt body is obtained directly in terms of the 
surface pressure distribution. 

Since the nondimensionalized enthalpy gradient at the sur- 
face is independent of surface temperature in this approxi- 
mation, the method is applicable to nonisothermal surfaces, 
so long as p./p.e >> 1. The calculations are particularly 
simple when the modified Newtonian flow approximation is 
employed, and they can be carried out for a blunt-nosed 
cylindrical body as well as for a body of revolution. In 
Section 3 the heat transfer rate solution of the present paper 
is applied to the special cases of an unyawed hemisphere and 
an unyawed blunt cone capped by a spherical segment. 

The methods developed here are also applicable to the 
second limiting case where diffusion is rate-controlling. The 
diffusion equation for each species is reduced to the same 
form as the low-speed energy equation, except that the 
Schmidt number u/pD,2 replaces the Prandtl number (Section 
4). For p./p, > 1, the maximum value of the ratio between 
the heat energy that would be transferred to the surface by 
diffusion alone, and by heat conduction alone is given by 
(pD12C,/k)u'/*. When the diffusion coefficient Dy, is esti- 
mated by taking reasonable values of atom-molecule collision 
cross section, this ratio is found to be 1.30. Additional 
theoretical and (especially) experimental studies of this 
problem are clearly required before these simple results can 
be accepted. 


2 Laminar Boundary Layer Equations for an 
Atom-Molecule Mixture of Perfect Gases 


2.1 Relative Magnitude of Heat Transfer by Diffusion and 
by Heat Conduction: The Parameter (pD,:2C,)/k 


Before examining the laminar boundary layer equations 
including diffusion in some detail, it is helpful to have some 
estimate of the relative magnitude of the rate of heat transfer 
by diffusion and by heat conduction. Because of the close 
similarity between the transport properties and atomic 
weights of oxygen and nitrogen, it is sufficient to treat the 
gas in the boundary layer as a binary mixture of atoms and 
molecules, so far as diffusion is concerned, in which case one 
deals with a single diffusion coefficient, D2. Also the con- 
centration of nitric oxide is generally low enough so that 
its contribution to the energy transfer is negligible. For the 
same reason one is justified in ignoring the presence of the 
ions and electrons in a discussion of boundary layer heat 
transfer at flight Mach numbers less than about 25. 

With these simplifications the heat transfer rate is given by 


oT oT ac; 


where V;, the diffusion velocity, is equal to —(D12/C;) (OC./dy) 
according to Fick’s law, and 


h; + hf 


5 This approximation was first employed by the present author 
in treating the hypersonic boundary-layer shock-wave interaction 
problem (3), but it appears to be generally useful. : 

6 The author is indebted to Dr. S. 8. Penner of the California 
Institute of Technology for this suggestion. 
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where h:° is the heat of formation of the ith species. Now, 
dh = CAT + Oh; where C, = 80 that 

In the limiting case of local thermodynamic equilibrium the 
atom concentration vanishes some distance away from the 


surface, provided the surface temperature is below the dis- 
sociation limit at the local pressure, p,. At the surface 


oc; oh oh, hs, — hw 
so that 
(Gu equil (hs hw) 


In the opposite limiting case where the volume recombina- 
tion rates in the gas phase are slow compared to the diffusion 


rates 
oy Jw 6 


for the atoms, on the assumption that the atom concentration 
is zero at the surface, while for the molecules 


oy 6 


Therefore 
conduction 
k w (hs ie hw) 
diffusion 


In other words, the relative magnitude of the heat transferred 
by the two processes is governed by the parameters 


Cy hilCie — Civ) 
( fu ~ 


At sufficiently high flight speeds this last parameter is nearly 
unity, and diffusion is responsible for practically all of the 
heat transfer. In this case, the surface heat transfer rate at 
a given flight speed (given h,,) is larger than the value for 
thermodynamic equilibrium roughly by the factor (pDi2C,/ 
bk)». When (pDy2C,/k)» = 1, the surface heat transfer rate is 
independent of the mechanism of heat transfer at any flight 
speed, and is exactly equal to the value for thermodynamic 
equilibrium.” 

According to the kinetic theory of gases the quantity 
eD,:C,/k is of the order of unity, where Dj; is the coefficient 
of self-diffusion. In fact, for nitrogen molecules u/pD;: = 
0.74 and c,u/k = 0.71 at T = 273°K and p = 1 atm, so that 
(eD,:cp)/k = 0.96.8 But here we require the value of Dj, 
which involves atom-molecule collisions rather than collisions 
between molecules; an estimate of this coefficient is made 
in Section 4.2. One expects Dj, to be larger than D;; if only 
because lighter and faster atomic particles are involved; 
therefore (pD,2C,/k)» > 1. For this reason, it is desirable 
to obtain a somewhat closer estimate of the relative magni- 
tude of the two modes of heat transfer by examining the 
boundary layer equations. Our procedure is to deal first 
with the case of thermodynamic equilibrium in some detail; 
many of the techniques developed here are immediately ap- 
plicable to the case where diffusion is rate-controlling for heat 
transfer (Section 4). 


_' A similar conclusion was reached by R. Bromberg using a 
different approach (private communication). 
® Table 1.2-3, page 16, Ref. (4). 
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In order to bring out the main features of the problem as 
simply as possible, consider the steady flow over an unyawed, 
blunt-nosed body of revolution (or a cylindrical body). The 
basic equations of laminar boundary layer flow are as follows 

Continuity Equation 


+ 2 (sont) = (1] 


where k = 0 for a planar body and k = 1 for a body of 
revolution.® 


Momentum Equation 
ou dp. ou 
Energy Equation 


+» 2 (4 4 
oy dy \Pr dy 


ra) 1\ 0 u? 1 ra) 


[3] 
Continuity Equation for Each Species 
oc; oc; 
10 
pu + oy (. ) [4] 
Equation of State 
p= T, where M = [= [5] 


h, = + u?/2, it includes the chemical enthalpy 
Chi. 
Clearly there are two cases in which the energy equation 
takes the familiar form: : 

1 Lewis number = pD,,C,/k = 1. 

2 C; = constant across the boundary layer. 

In the first case, as shown previously, the heat transfer 
rate is given by the expression 


C, dy 


independently of the mechanism of heat transfer. In the 
second case, either the gas temperatures are low enough so 
that no dissociation has occurred, or no chemical reactions 
take place either in the gas phase or at the surface; i.e., the 
gas composition is “frozen.” 

Apart from these general remarks not much progress can be 
made until one finds the “‘similar solutions,” or one-parameter 
family of velocity and total enthalpy profiles. In this respect 
the situation is not very different from the low-speed boundary 
layer problem. 


2.2 Similarity Considerations for Thermodynamic Equi- 
librium: The Case of the Highly Cooled Surface 


In the limiting case of thermodynamic equilibrium the con- 
centration of each species is uniquely related to the local 
pressure and temperature, and the net rate of production of 
each species w; is supposed to be large enough to balance the 


® For Reynolds numbers based on nose radius that are not 
“too low,” 5/ro << 1 and the continuity equation takes this re- 
duced form. Also, 3é and 6703C/ds are assumed to be “small,” 
where 3 is body surface curvature. The boundary layer and 
shock wave are distinct when 50/5. ~ (pe/po)/VWualo/vo < 
0.10 (say), or > 104. (Here 6, is the 
separation distance between the bow shock and the forward 
stagnation point. ) 

10 These equations are readily derived from the general equa- 
tions given by Hirschfelder, et al. (Ref. 4). Thermal diffusion is 
considered to be negligible compared with mass diffusion. 
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convection and diffusion terms. In other words, Equation 
[4] is superfluous in this case. Also, tables or graphs of static 
enthalpy and molecular weight vs. temperature for air are 
either available or can be constructed by methods similar to 
those described in (5b). By means of the relation p/p, = 
(T./T) (g/m.) the density ratio is then uniquely deter- 
mined by the enthalpy ratio h/h., for each value of h, and 
the pressure p, across the boundary layer. 5 

Our procedure is to assume that “similarity” exists and 
then to determine the conditions (if any) under which this 
assumption is in fact justified. For this purpose the momen- 
tum and energy equations are reduced to the more convenient 
equivalent “low-speed” form by utilizing Levy’s transforma- 
tion (6), which he introduced for the special case of constant 
specific heats."1 In order to simplify the discussion the 
Lewis number is taken to be unity; this restriction is not an 
essential one, and in no way affects the validity of the general 
arguments. 

Let 


Pelle y 
ro’ p/pe dy 
and 

The transformation from s, y coordinates to §, » coordinates 
is carried out by means of the relations 

_ 2 

dy —(28)'/* 


[2 


The over-all continuity equation is automatically satisfied 
by introducing the stream function y, which is defined by the 
usual relations 


and 


Let (8, ») = (28)'/* f(n); then by Equations [7] and [8a], u/u. 
= f'(n), where the prime denotes differentiation with respect 
to ». Also define h,/h,, = g(n). The momentum equation 
then takes the form 
25 du. Pe 
Tr +7" + 


and the energy equation is 


r’) + fo’ + = =0....[10] 


Pebte 


The boundary conditions are as follows: 


f(0) = f'(0) = 0, and f(y) > asy 
9(0) = gu(s), or g'(0) = 0, while + 1 asn > ~. 


where C= 


For an isoenergetic flow outside the boundary layer the 
conditions for similarity are: 


(a) (25/ue) (du,/d8) — (f’)*] = F(n), or a const. 
{b) C = C(n), or a const. 


1 Because the gas is treated as a binary mixture of atoms and 
molecules, the mass ratio of oxygen (atoms + molecules) to 
nitrogen is constant along every streamline and equal to the 
stoichiometric ratio for the standard atmosphere. 

11 Here we combine Levy’s transformation with the Mangler 
transformation for bodies of revolution. 
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(c) Either (a) Pr = 1 or (b) Pr = Pr(n), ora const, while either 
u,*/2h,, < 1 (low local Mach number) or u,? — 2h,, (high 
local Mach number). 

(d) gw(s) = const. 


Except for the special case of a uniform external flow, all of 
these conditions for similarity are never satisfied simul- 
taneously at hypersonic flight Mach numbers, even if the 
surface temperature is uniform. The density ratio p,/p, 
the Prandtl number, and the ratio pu/p.u. are functions not 
only of h/h., but also of h, itself, and of p, through its effect 
on the composition of the gas. Of course it is always possible 
to find a “locally similar” solution in the neighborhood of the 
forward stagnation point that is “exact” so long as the 
boundary layer and bow shock are distinct. The boundary 
layer thickness h, and p, are all very nearly constant along 
the surface in this region, and the velocity and enthalpy 
ratios are functions of a single variable of the form )/R, 


V/ where Y = (p/p.) dy. The question na- 


turally arises as to whether one can obtain approximate 
solutions extending over most of the body surface. Such 
solutions can be developed in at least two cases of interest: 
(a) slowly varying external flow properties; (b) highly cooled 
surface, or go 1. 

In the first case the approximation of “local similarity” is 
employed, in which the velocity and enthalpy profiles at any 
station are supposed to be identical with the “similar profiles” 
corresponding to constant values of 


B= ( ) and 
loca 


all along the surface, and to the local relations between p,/p, C, 
Pr, and the enthalpy ratio h/h.. This approximation is 
utilized by Stine and Wanlass (7) (for example) in the special 
case of constant specific heats and ordinary stagnation tem- 
peratures, where p./p = h/h, and C and Pr are both practically 
constant. 

In the second case, where the surface temperature is much 
lower than 7, or p./p» < 1, the influence of the pressure 
gradient term in Equation [9] should be small according to 
the qualitative physical argument in Section 1, and the prob- 
lem is simplified considerably. This argument is supported 
in two ways: (a) rough estimates; (b) examination of known 
solutions for the special case of constant specific heats. At 
hypersonic flight speeds 


whe 


when the difference in molecular weight across the boundary 
layer is taken into account. Now 
9 = + (1 — = 

when g»<1.!2 Therefore, the term in brackets in Equation 
[9] is represented approximately by [u/u. — (u/u.)*], which 
is zero both at the surface and at the outer edge of the bound- 
ary layer, and has a maximum value of 1/, within the bound- 
ary layer. Thus, at most, the effect of the pressure gradient 
of the velocity profile is equivalent to the case of an iso- 
thermal low-speed flow with an “effective” value of 6 = 
(28/u.) (du,/d8) equal to one fourth of the actual value of 
this parameter.}* 

12 This relation is “exact’’ when the pressure gradient is zero, 
and either Pr = 1, or (u,*/2h,,)< 1. 

13 The first-order effect of the pressure gradient term on the 
skin friction, for example, is given by 


Pe 


where = (p/p) — (f’)? — (u/u.)* and f(») is the 
Blasius function. 
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At the forward stagnation point the parameter (23/u,) 
(du,/ds) is equal to '/, for a body of revolution and unity 
fora cylinder. The maximum “effective” values of 8 when 
bw < 1 are and 1/4, respectively, and the influence of 
the pressure gradient on the velocity profile is much reduced, 
compared to the case of an insulated surface. But the effect 
of the pressure gradient on the heat transfer rate is even 
smaller, because the change in g’(0) from itz zero pressure 
gradient value depends upon an integral of the change in the 
yelocity profile across the boundary layer. In fact 


Iw 
foe dn 
These conclusions are borne out by the results of the nu- 
merical calculations of Cohen and Reshotko (8) for the 


“similar solutions,’ with the particular assumption of con- 
stant specific heats and Pr = 1,C = 1. In that case 


pe/p = 9 + — 1)/2 M.? [g — (f')*) 
and Equations [9] and [10] become 


9'(0) = 


and 


a”’ + fg’ = 0 


with 


B = (23/M.) (dM./ds) = const 


The surface shear stress and heat transfer rate functions 
fe’ and gw’/1 — gv = —8w'/8» as calculated by Cohen and 
Reshotko (8) are shown as functions of @ in Figs. la and 1b, 
for the cases g(0) = h./h,, = 0, 0.20, and 1.0.14 Clearly the 
approximation that the direct influence of pressure gradient 
is negligibly small for p./p» << 1 is highly accurate for the 
surface heat transfer rate, and somewhat less accurate for 
skin-friction."*5 Similar conclusions can be drawn from the 
numerical results obtained by Levy (Fig. 10 of Ref. 6), and 
by Lal (Fig. 5 of Ref. 2). 

With regard to the function C(n) appearing in Equations 
[9] and [10], extrapolation of existing data on the viscosity 
of air by means of Sutherland’s formula shows that at high 
enthalpy levels (30-100 times ordinary enthalpy), pu = pete 
over a large portion of the boundary layer (see Ref. 5c, for 
example). Therefore so far as the heat transfer rate is con- 
cerned, it is sufficient to take C = 1, especially when we recall 
that C enters the problem as a square root. Also, in the 
present state of our knowledge, it seems reasonable to re- 
gard the Prandtl number as a constant, equal to some aver- 
age value Pr across the boundary layer (Section 4.2). 

With regard to the energy equation (Sp. [10]), the solution 
can always be chosen so that the inhomogeneous last term 
does not contribute to the enthalpy gradient at the surface; 
in fact, this term determines the “recovery” enthalpy, or sur- 
face enthalpy for no heat transfer, as shown by Chapman 
and Rubesin (9). Of course this term vanishes at a stagna- 
tion point, in which case the recovery enthalpy is identical 
with h,,. In any event this term is numerically small on a 
blunt-nosed body at hypersonic flight speeds, where 
(u.?/2h.,)max & 1/2 and Pr & 0.65-1.00. 

These considerations suggest that the nondimensionalized 
enthalpy gradient g’(0) at the surface of a blunt-nosed body at 
hypersonic flight speeds is represented quite accurately by 
the zero pressure gradient value of 0.47 Pr'/*, independently 
of surface temperature, local Mach number M,, or pressure 
gradient, and is therefore applicable to a very general class 
of surfaces, whether isothermal or not, so long as pw/p,. >> 1." 


'® From Figs. la and 1b it pee that the “effective” value 
of 8 is more nearly one sixth of the actual value. 
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Fig. 1(a) Effect of pressure gradient and ratio of surface tem- 
perature to stagnation temperature on surface shear stress 
function, Pr = 1,C = 1 
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Fig. 1(b) Effect of pressure gradient and ratio of surface tem- 
perature to stagnation temperature on surface heat transfer rate 
function, Pr = 1,C = 1 


If other quantities are required which depend more sensi- 
tively on the distribution of pu and p across the boundary 
layer, they can be calculated by means of a straightforward 
iteration process. In the first step the stagnation enthalpy 
distribution g(n) is obtained from Equation [10] by re- 
peated quadratures, with f(n) given by the Blasius function, 
u,?/2h,, equal to its local value, C = 1 and Pr = Pr. (Here 
the local similarity approximation is employed.) By sub- 
stituting the corresponding density profile p,/p and the 
function 


Cn) = 
Pe Me 


into the momentum equation (Equation [9]), the next ap- 
proximation f(y) is calculated by numerical integration for 
the local value of 

25 du, 

ds 


and the process is repeated, if necessary. Finally, the actual 
distance normal to the surface is computed by means of the 


quadrature 
93)'/2 
= { pe/p dn 
0 


k 
To Pelle 


and the displacement thickness (for example) is given by the 
expression 


16 This approximation is clearly not applicable near boundary 
layer separation. 
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With the aid of modern computing machines and tabulated 
air properties, Equations [9] and [10] can be integrated 
simultaneously step-by-step for any constant values of 8 and 
2.3 Surface Heat Transfer Rate Distribution 


In the limiting case of thermodynamic equilibrium, all of 
the heat transfer at the surface is accomplished by conduction, 
and 


kw 


for pw/pe >> 1, taking round numbers. Now pupw = pele, to 
be consistent, and 


(Pette po’ 


where w, = u./&.T,. Therefore 
s 
28 = pette p/po' “* ds 
0 Ue Wey 


and 


Ge 0.50 (Pr)—*/* the F(8). [12] 
where 


Lf (p/po') (we/we,) 
0 


Here k = 0 for a planar body and k = 1 for a body of revolu- 
tion, and Pr is some average Prandtl number [see Section 
(4.2)]. In particular at the forward stagnation point 


& 0.50 2*/2 (Br)—*/* W/(pettelo V theo hee 
V (du./ds) 


Some years ago Squire (10) worked out the solution for the 
convective heat transfer near the forward stagnation point of 
a cylinder on the assumption of constant fluid properties 
(small temperature differences). More recently, Sibulkin (11) 
obtained an analogous solution near the forward stagnation 
point of a blunt body of revolution. Their results are ex- 
pressed in the form 


Nup = = 0.570 (Pret (4/3) 
ke uz ds v2 


where go = h(T., — T.), and the subscript 2 denotes quan- 
tities evaluated in the gas just behind the normal shock, fol- 
lowing Sibulkin’s suggestion.” For purposes of comparison, 
define h in our case by the relation (gu)o = (h/Cy,) (hu — 
hw); then our result is expressed as 


Nup,, = AD = (0.50) 2*/2 (Pr)'/* (2 
uz ds /o Ve 
As expected, the present result has the same form as the 


F(s) = 


1 Cohen and Reshotko (12) have applied their similar solutions 
for a compressible fluid in the ial case of constant specific 
heats and C = Pr = 1 to the forward stagnation point, with 
— results. They also investigated the effect of fluid injec- 

ion. 


“low-speed,” constant-property solutions, except that now 
the value of [(D/ue)-(du./ds)]o is obtained from the New. 
tonian flow approximation at hypersonic speeds, and not 
from the irrelevant low-speed potential flow over a cylinder 
or a sphere. 

The expression for Nup, . at the forward stagnation point 
is of the same form as the result obtained by Beckwith (13), 
who utilized the momentum and energy integral method, 
including an approximate density-enthalpy relation for dis- 
sociated air in thermodynamic equilibrium. 

Also, R. Mark (14) worked out the differential equations 
for the similar solutions in the neighborhood of the forward 
stagnation point, and integrated these equations numeri- 
cally, utilizing the equilibrium thermal properties of disso- 
ciated air given in (5c). (Lewis number = 1.) When his 
result for the stagnation point heat transfer rate is expressed 
in the same form as our Equation [12a], then his quantity 
\[(Eq. (14) and Fig. 5 of Ref. 14)], multiplied by puu./p.u,, 
should be compared with the numerical value 0.50 V/2 = 
0.70 (more precisely with 0.47 V/2 = 0.67). For surface 
temperatures of 3000 and 5000 R, and over a range of values 
of stagnation enthalpy from 50-120 times the ordinary 
static enthalpy of air, the value of \(pupw/peue) differs from 
0.67 by 5 per cent at most. For 7, = 1000 R, the value of 
this quantity is about 10 per cent lower than 0.67. This re- 
sult tends to justify the simplifying approximation py = 
Pete Utilized in this paper. 5 

As discussed in Section 4 of (15), a careful examination of 
the experimental data on pressure distributions over hemi- 
sphere cylinders and blunt cones with spherical noses shows 
that the modified Newtonian law in the form C,/C,,,,. = 


sin? @ is confirmed for M, > 2 (approximately). Therefore, 
in Equation [12], p/po’ = sin? 6, + (p../po’) cos? %, and 


Ue? 1 M,,? Po’ 


where 7 is the mean ratio of specific heats behind the bow 
shock wave. (At high temperatures, 7 = 1.10-1.20.) Near 
the forward stagnation point sin? 6, = 1 — sin? @, so that 


~ 1 M,,! 


where p,./po’ 1/7..M,,? and 


The stagnation point heat transfer rate is then 


7; sind +... 
YoU w 


0.50 2/2 (Pr)- V (pette)o Vu. hee*G(M ) 


VR 


(Gw)o = 


where 


y¥-1 “ 2 1 


1 
( 


In flight at hypersonic speeds h,, = u,2/2; in the wind 
tunnel h,, is usually the reservoir stagnation temperature; 
in the shock tube h,, = u,? for sufficiently strong shocks, where 
u, is the velocity of the shock wave in the straight section. 


8 This expression agrees with the experimental data (16) ob- 
tained on a hemisphere cylinder to within 3 per cent in the range 
1.735 M. 5, with 7 = 7. = 1.40. 
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3 Surface Heat Transfer Rate Distribution over 
Unyawed Hemisphere and Blunt Cone 


According to the modified Newtonian flow approximation, 
for the particular case of an unyawed body of revolution with 
constant radius of curvature in the meridian plane, p/po’ = 
cos? 0 + (1/7.M.,2) sin? 6, where @ is the angle between the 
radius vector and the flight direction. By retaining the 
second term in this last expression the heat transfer rate 
distribution obtained is quite accurate even for Mach numbers 
as low as 2.0, and joins smoothly with the distribution over 
an afterbody of arbitrary shape. 

The velocity at the outer edge of the boundary layer can 
always be obtained from the isentropic flow relation with an 
efiective 7. But the predicted velocity over a hemisphere 
varies very nearly linearly with @ up to @ = 80 deg, and the 
experimental data for air in the range 1.97 < M,, < 5.8 show 
a remarkable linearity (15, 16). In the present ‘calcula- 
tions, thereore, we take u./u. = (1/w.) 9. On 
the conical skirt 


Because of the slow variation of the quantity o. = u./ 
R.T, with enthalpy h, at high enthalpy levels it is sufficient 
to take w & ,,. With this simplification the ratio of the 
local surface heat transfer rate to the heat transfer rate at the 
forward stagnation point is given by 


Go _ __(1/2) (p/po') 
(Geo) (p/po') (Ue/U as| 
VR 
V (1/tq) (due/d8)o 


Case I Hemisphere: 
For this body ds = Rod6, ro = Ro sin 6; therefore 


qe 1/.@ sin [cos? + (1/7,,M,,”) sin? 6] 
ee } f, o (cos? + (1/7_M sin?6] ao} 


or 
Ge _ Wein {[1 — cos? 
where 
1 — cos 40 
= (1 - 75) (#- )+ 
6? — sin 20 + [15a] 


Near the forward stagnation point 


qu 
—— = 1 — { 0.722 
(Gu)o ( 

In Fig. 2 the ratio g./(gu)o is plotted as a function of @ for 
M., = 2, 3, 5, 10, and ~, with y, = 1.40. The quantity 
V p/p! for M ,, = 2 is also shown in Fig. 3 for comparison. 
Evidently the heat transfer rate drops off somewhat faster 
than V/p/po! ~ V pate. By dividing by we 
obtain a quantity proportional to the local heat transfer 
parameter 


Yo 


Nu/ / Re = ———— 
V petted / te 
employed by Korobkin (16) and by Stine and Wanlass (7). 
For M,, = 2, the values of Nu/~/Re calculated by the 
present method for 7../7. <1 are about 12 to 15 per cent lower 
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Fig. 2. Laminar heat transfer rate distribution over isothermal 
hemispherical nose, with (T'w/T.)<1 
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Fig. 3 Laminar heat transfer rate distribution over blunt cone, 
with (Tw/T.)<1 M,, sin 1 


than the theoretical values obtained by Stine and Wanlass 

for T/T. & 1 (solid curve in Fig. 4b of Ref. 7). This dif- 

ference is almost entirely accounted for by the fact that the 

nondimensionalized enthalpy gradient at the surface for 

T../T. = 1 is 12 to 15 per cent higher than the Blasius valya_ 
(See Fig. 1b, 8 = '/2; go = 1, for example.) Unfortunately- 
the scatter in the experimental data (7, 16) so far obtained is. 
too great to perrizt any detailed comparisons to be made, 
with the theoretical results. 

Except for the portion of the hemisphere aft of @ = 60 deg 
(approx.) the heat transfer rate distribution is close to the 
“frozen” hypersonic distribution when M,, > 2. For M, > 
5 the heat transfer rate is close to its hypersonic values over 
practically the entire hemisphere. Actually the static pres- 
sure is somewhat higher over the region near @ = 90° than 
predicted by the Newtonian approximation, so the predicted 
values of g./(Gw)o there should be regarded mainly as indica- 
tive of relative magnitudes. 

Case II Blunt Cone Capped by a Spherical Segment: 

For this body the surface heat transfer rate distribution 
over the “nose” is of course identical with that found for the 
forward region of the hemisphere [9 < +/2 — 0@,, where 0. is 
the cone half angle (see Fig. 3)]. On the conical skirt all 
physical quantities are regarded as constants, and ry = 
s’ sin 6., where s’, the distance along the suriace measured 
from the “virtual tip,” is related to the actual distance s by 
the expression (Fig. 3) 


By substituting the expressions for p/po’, u./u,, and 7 into 
Equation [14], one obtains 


8’/Ro 
[B(8.) + (s’/Ro)*]'/* 


qu 
— = A(@ 
for s’/Ry = cot 6., where 
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3/16 
sin? @, {[1 — (1/yoMu*)] sin®@, + 
9 jo=(x/2)—0, 
and D(@) is defined by Equation [15a]. 
At the junction of the nose and the conical skirt, s’/Ry = 


cot 6., and this expression for gu/(Gw)o agrees with Equation 
[15] for @ = 1/2 — 6. “Far” from the junction (s’/Ry > 1) 


Ge > V3 (0.35) V peste Vue 1/V8’) 


which is exactly the value of the surface heat transfer rate 
on a sharp-nosed cone with an “equivalent velocity” 


due 
Ue = — @) 


= 


— cot? 6, 


When sin 6. >> 1, A(0.) & VW3/2 sin 6.--W(x/2) — 6. 
and 


1 
B(6.) = */16 ante cot? 6, 


The numerical values of these functions for @ = 10 deg, 20 deg, 
30 deg, and 40 deg are tabulated below. 


.° A(6.) B(6.) 
10 0.181 172.0 
20 0.326 +4.2 
30 0.443 —0.20 
40 0.520 —0.26 


Evidently for 6. = 30 deg and 40 deg the surface heat 
transfer rate distribution over the conical skirt is practically 
identical with the distribution for the “equivalent sharp-nosed 
cone.” For slender cones, however, the heat transfer rate is 
smaller than the value for a sharp-nosed cone for a distance 
of several nose radii aft of the junction of nose and skirt. 
These conclusions are illustrated quantitatively in Fig. 3. 


4 Diffusion as Rate Controlling for Heat Transfer 


4.1 Surface Heat Transfer Rate 


When the volume recombination rates in the boundary 
layer can be disregarded in comparison to the rate of diffusion 
across streamlines, the diffusion equation for each species 
takes the following form (Equation [4], Section 2.1, with w; = 
0): 


ac; eC: ac; 
-2 (ons [16] 


By introducing the transformations of Equation [7], and de- 
fining z;(7) = C./C;, for the atoms, and 


2i(n) = 


for the molecules, Equation [16] is reduced to an equation 
for z;(n) of the same form as the energy equation (Equation 
{10]) with the term proportional to u,?/2h,, suppressed, except 
that the Prandtl number is now replaced by the Schmidt 
number, or the parameter u/pD,,.. By the arguments em- 
ployed in Section 2.2, it follows that when pw << 1, 2;’(0) & 
0.50 (u/pD,2).'/*, and the rate of heat transfer to the surface 
by diffusion is 


(Gu )aite = (pD12)w (2), [> hifCi, 


By comparing this expression with Equation [11], one obtains 


hin (Cig — Cin) 
Nee 


(Gu 

equi k 
where (Gw)equit is the rate at which heat is transferred to the 
surface by conduction in the limiting case of thermodynamic 
equilibrium. In other words, at a given flight speed the mazi- 
mum difference in surface heat transfer rate between two limiting 
cases is measured by the factor (pD,2C p/k) w’/*. 

For the particular case in which the surface temperature is 
high enough so that all the C,,’s are nearly the same (while 
Pw/ pe > 1), the enthalpy transferred by conduction can be 
treated separately from the full enthalpy transfer. The full 
static enthalpy equation for the laminar boundary layer is as 
follows: 


= +0 = o(v + = 


re) oT ou \? dp, 
(: or) + tun [18] 


By utilizing Equation [16] and the fact that }>pc:V; = 0, 
Equation [18] is reduced to the familiar form 


oT oT re) oT du \? dp 
pu Cy + = 2 (027) 


where 
r 


Without going into details it is clear that in this case the sur- 
face heat transfer rate by conduction alone is given by 


(Ge)eond. 0.50 (Pre)'/* Kells — (cis — 


and 


(Gw)eond. + (Guw)ditt ~ Cie 
(Gw)equit 


‘Se 


In this particular case the differential equation for the 
partial enthalpy 


w ha, 


+ u?/2 
hre + Ue?/2 


is identical with Equation [10], and the temperature dis- 
tribution is obtained by a method similar to that described 
in Section 2.2. The concentration profiles are obtained by 
numerically integrating the equation 


) = gr(n) 


a” + (+) fo'=0 with f=sf%e 
fa(n) and = 1" 

and the boundary conditions z;(0) = 0, z:(n) ~ 
Once z;(n) is known, the density profile is obtained from the 
relation 

_ mT 

m T. 
where = Of course iteration schemes 
similar to the method described at the end of Section 2.2 can 
be developed in this case also. 


When the cp,’s are not all nearly the same the problem is 
somewhat more complicated and the details will not be given 


19 fp() is the Blasius function. 
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here. In any event the ratio of the total surface heat transfer 
rates in the two limiting cases considered will lie between the 
yalue given in Equation [19] and (pD,:C,/k)u'”. 


42 Estimate of the Diffusion Coefficient Dy. and the 
Parameters and €,/k) According to the 
Kinetic Theory of Gases 


In recent years several different groups of investigators have 
been remarkably successful in predicting the transport proper- 
ties of dilute, nonpolar gases” and gas mixtures with the aid 
of approximate interaction potentials for the intermolecular 
force field. Hirschfelder and co-workers have employed the 
two-parameter Lennard-Jones potential 


where r is the distance between the centers of the two mole- 
cules, ¢ is the depth of the potential “well,” or the maximum 
energy of attraction, and o is the zero-energy distance, or 
separation distance at the point of closest approach for low- 
energy collisions (xo? is the collision cross section). The 
inverse sixth-power term takes into account the long-range 
attractive forces, while the inverse twelfth-power term ap- 
proximates the strong, short-range repulsive forces that are 
called into play when the electron clouds surrounding the 
two molecules begin to interpenetrate.*! 

For any potential function of the type g(r) = «f(r/c) the 
first approximation to the diffusion coefficient of a binary 
mixture is given by 


TX™M; + Me) 1 [20] 
2MiM2 


in cm?/sec, where p is the pressure in atmospheres; 7’ is the 
absolute temperature in °K; ™;, Mare the molecular weights 
of species 1 and 2, respectively; 712* = kT'/e2, where k is the 
Boltzmann constant; 042, €12/% are the interaction parameters 
for collisions between species 1 and 2, in A and °K, respec- 
tively, and Q,.%s* (7).*) is the collision integral for mass 
transport; this integral is a measure of the departure from 
the rigid-sphere model, for which 2).* = 1.0.22 The higher 
approximations to Dy, are given by (Diz)k = (Diz): fou, 
where fp,“ depends on the mixture ratio, but for most 
binary mixtures 1.00 < fp,‘ < 1.03 for the Lennard-Jones 
potential. In other words, D,2 is virtually independent of the 
mixture ratio. 

For low energies [7:2* = 0(1)] the collision integral for the 
Lennard-Jones potential is larger than unity, while for high 
energies Q;24* < 1, The Lennard-Jones molecule has a 
larger cross section than the rigid-sphere model at low 
energies, because of the importance of the attractive forces, 
and a smaller cross section at high energies to allow for the 
interpenetration against the repulsive forces. The function 
2", 0*(7'"*) is tabulated for 0.30 < 7* < 400 in Table 1-M, 
pp. 1126-1127 of the Appendix to (4). 

In principle, measurements of the diffusion coefficient pro- 
vide the best source of information about interactions between 
unlike molecules, but such measurements are usually very 
limited. In general the interaction parameters 2 and «2 
are approximated by the combining “rules” 12 = (0; + o2)/2 
and = V «1 €2, Where oj, and are determined from 
measurements of the viscosity of the pure components of the 
binary mixture. As we shall see, the value of oi: for atom- 
molecule diffusion determined in this manner is too low, be- 


(Dh = 


_* Nonpolar molecules are those for which the potential func- 
tion does not depend on the orientations of the two interacting 
molecules, but only on the distance between centers. 

*1 An excellent account of the whole modern development of 
the molecular theory of liquids and gases is given in (4). 
22 For the rigid-sphere model g(r) = 0 for r > a, and ¢(r) = 

forr<o. = 1 Angstrom unit = 10-*cm. 
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cause this method does not take into account the important 
physical differences between atom-atom and atom-molecule 
interactions. 

Fortunately a considerable amount of theoretical work has 
been done on the hydrogen atom-hydrogen molecule inter- 
action, and also on the molecule-molecule interaction in 
hydrogen. Experimental data on the viscosity of molecular 
hydrogen provide a check on the theoretical results obtained 
for the second case. By studying these results for the 
simplest element, we are able to make a reasonable estimate 
of the difference between the known collision cross section 
for molecule-molecule interactions in air, and the required 
collision cross-section for atom-molecule interactions. 

Even for the “simple” hydrogen molecule the calculation of 
the short-range repulsive forces in the H,-H, interaction is so 
complicated that reasonable approximations must be intro- 
duced (see Ref. 4, pp. 1083-1092). Four possible orienta- 
tions of the molecules are considered, as follows 

45° 
(d) 


When the spatial average of the potential over all four orien- 


(a) (b) (c) 
tations is calculated, it is found that a Lennard-Jones function 
fits this potential very closely, provided « = 2.928 A and 
e = 5.107 X 10-“ erg. It is remarkable that these values 
agree to within 1 per cent with the values of o and ¢ obtained 
from viscosity data on molecular hydrogen. 

On the basis of rough physical considerations, one might 
expect that the “collision diameter,” ox-n,, would be smaller 
than ou,-n, by a distance of the order of one half the equilib- 
rium separation r, of the atoms in H:, or about 0.37 A. 
Margenau (17) considered the atom-molecule configurations 


H 
H 
fo and 
H 


(A) (B) 

He calculated the repulsive part of the potential in each case 
by a first-order perturbation method, while the attractive 
portion of the potential is approximated quite accurately by 
the sum of an inverse sixth-power and an inverse eighth- 
power term. For configuration A, on-n, = 2.65 and 
en-p, = 5 X 10-* erg; for B, cu-n, = 3.03 A and en-z, = 
3.5 X 10-“ erg. Configuration A appears to be much more 
probable, and the potential is surprisingly similar to the 
potential for configurations (b) and (c) in the H,-H, inter- 
action, at least in the range 2.50 A < r < 4.0 A. Thus the 
collision diameter for atom-molecule interaction in hydrogen 
is about 0.28 A smaller than CH,-H,, OF OH,-H, — CH-H, 
3/4 (r./2), while the maximum attraction energy in the 
two cases is virtually identical. 

Now the best experimental value of the distance between 
the two nitrogen atoms in the N; molecule is about 1.1 A. 
On the basis of the results of the study of the H-H, and H.-H; 
interactions, it seems quite reasonable to suppose that the 
collision diameter for atom-molecule collisions in air is smaller 
than the value « = 3.69 A®* for molecule-molecule interaction 
by about 3/4(r,/2), or about 0.4 A. In other words, 01. & 
3.69 A —04A =33A. Also = = 84 K, at 
least for 300 K < T < 1900 K. The diffusion coefficient D,. 
is to be obtained from Equation [20] with these values of o12 
and e2/k.?4 

Suppose we consider first the viscosity and conductivity of 
the gas at the surface, where the concentration of atomic 


23 Table I-A, og Ref. (4). 
e 


24 In Ref. (18) the value of o. is obtained from the rough ap- 
proximation oi. = (0; + o2)/2, and is taken as !/2 i.e., 
oi = */, o: & 2.77. This value seems to be far too low. 
(Here the subscript ‘1’? denotes the atomic and the subscript. 
“2” denotes the molecular species.) 
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species is zero. Experimental data on the viscosity of 
molecular air are reproduced to within 1 per cent in the tem- 
perature range 300 K < 7 < 1500 K, and to within 3 per cent 
in the temperature range 1500 K < T < 1900 K, by the fol- 
lowing expression based on the Lennard-Jones potential: 


Vv mT 
2)* (T2*) 


where 0, 2 *(7',*) is the collision integral for momentum 
transport. By Equations [20] and [21] the Schmidt number 
for atom-molecule diffusion in air is 


M2 D*(T2*) 1 
2 {1 + M:/™:] 


where 0.833 is the value of the Schmidt number for the 
idealized case of self-diffusion in a monatomic gas according 
to the rigid-sphere model. But €12/k & e2:/k = 84 K, so that 
T.* = T,*, and the ratio of the collision integrals = 1/ 
A*(T*), where A*(7"*) is a very slowly varying function of 
temperature given in Table 1-N, Appendix, Ref. (4). Also 
Ms = 2M, and o;2. = 3.3 A, while o2 = 3.69 A, so that 


( ) 0-543 

~=A*(T*) 

A few calculated values of the Schmidt number are tabulated 
below: 


X 10" = 266.93 in gm/cm-sec. . . [21] 


T, °K 
252 
420 
840 

1680 

3360 


Clearly the variation in the Schmidt number over this tem- 
perature range is less than the probable inaccuracies in the 
estimated value itself. 

For a monatomic gas or a mixture of such gases the thermal 
conductivity is predicted by the kinetic theory with the same 
accuracy as the viscosity coefficient. But for polyatomic 
gases the question of the exchange of energy between the 
translational and internal degrees of freedom enters the 
problem. For dry air the vibrational relaxation times for 
Nz and O, are known to be long at temperatures below 
1000 K; thus, one might expect the thermal conductivity to 
behave as if only the translation and rotational degrees of 
freedom were involved. In other words the Prandtl number 
should be practically constant at its value for room tem- 
perature up to about T = 1000 K. According to the NBS- 
NACA Tables (Table 2.44, Ref. (5a)) the Prandtl number 
for dry air based on the best fit to experimental values of 
uw and k, and spectroscopic data on C,, varies from 0.708 at 
300 K to a minimum of 0.680 for 500 K < 7 < 600K, and 
then rises again to 0.702 at 7 = 1000 K. This variation in 
Prandtl number is less than the indicated scatter in the ex- 
perimental data on thermal conductivity in Fig. 1, p. 5 of 
Table 2.42, Ref. (5a). Of course at somewhat higher tem- 
peratures the vibrational mode is much more readily excited, 
and the Prandtl number should increase gradually to the 
undissociated Eucken value of 0.783. Unfortunately no 
experimental data exist at temperatures above 1000 K, and 
the Eucken formula is known to give values of Prandtl num- 
ber at 7 = 300 K that are about 5 per cent too high. For 
the-present, we take Pr = 0.71. 

Finally the parameter (pD,2C,/k)» governing the relative 
importance of heat transfer by atom-molecule diffusion and 
by thermal conduction is given by 
(Cou / oy 071 1:45 

k Jw ‘(u/pDi2)e (0.49 


and the maximum value of the ratio of the rate at which heat 
is transferred by diffusion and by heat conduction is (1.45)*/, 
or about 1.30. 

In Sections 2.2 and 4.1 it was pointed out that the effect 
of the variation in Schmidt number, Prandtl number, and the 
quantity C(n) = pu/peu. across the boundary layer can be 
taken into account by means of straightforward iteration 
schemes, or by simultaneous integration of the equations of 
motion with the aid of computing machines. For the case of 
thermodynamic equilibrium the properties of air are given 
in (5c). However, the calculations of air properties in this 
report are based on values of atom-atom collision diameter 
~1.5 A, which are far too large. According to the quantum- 
mechanical treatment of the N-N interaction by Kopineck 
(19), the equilibrium separation distance is 1.109 A, which 
agrees very well with the experimental value of 1.095 A, while 
the collision diameter o; = 0.80 Also «,/k = 1.14 x 10, 
On the other hand the values of atom-molecule collision 
diameter utilized in (5c) and (18), namely o42 = (0: + o2)/2 = 
2.8 A, seem to be too small, if the estimate of the present 
paper is correct. While the effects on Prandtl number of 
the difference-in these two sets of estimated values of o; and 
o12 is small, the effect on the Schmidt number and C(n) may 
be significant. 

For the second limiting case in which the volume recom- 
bination rates are regarded as negligible in comparison with 
the diffusion rates, the viscosity of the gas considered as a 
binary atom-molecule mixture is calculated from Equation 
8.2-30, page 533 of (4), while the Prandtl number is obtained 
by the method described by Hansen (18), except that the 
correct values of 01, o12, and o2 should be employed. The 
diffusion coefficient D,. is obtained from Equation [20] of 
the present paper. The following values are recommended 
for the interaction parameters. 


a/k = 1.14 X 10K 
e2/k = e/k = 84K 


“Hj = 0.80 A, 
on = 3.34 


= 3.69 


References 


1 Cohen, C. B., and Reshotko, E., “The Compressible 
Laminar Boundary Layer with Heat Transfer and Arbitrary 
Pressure Gradient,’’ NACA TN 3326, 1955. 

2 Lal, S., “Heat Transfer in Compressible Laminar Bound- 
ary Layers,’’ Ph.D. Thesis, GALCIT, June 1955. 

3 Lees, L., “On the Boundary Layer Equations in Hyper- 
sonic Flow and Their Approximate Solutions,’’ Journal of the 
Aeronautical Sciences, vol. 20, no. 2, February 1953, pp. 143- 
145. 

4 Hirschfelder, J: O., Curtiss, C. F., and Bird, R. B., 
“Molecular Theory of Gases and Liquids,’’ John Wiley and Sons, 
Inc., New York, 1954. 

5 (a) The NBS-NACA Tables of Thermal Properties of 
Gases, Heat and Power Div., NBS, Wash., D. C., February 1, 
1951, and January 1, 1952. 

(b) Wooley, H. W., “Effect of Dissociation on Thermo- 
dynamic Properties of Pure Diatomic Gases,’? NACA TN 3270, 
1955. 

(c) Romig, M. F., and Dore, F. J., “Solutions of the Com- 
pressible Laminar Boundary Layer Including the Case of a Disso- 
ciated Free-Stream, ‘Convair Report no. ZA-7-012, Aug. 4, 
1954, San Diego, Calif. 

6 Levy, 8., “Effect of Large Temperature Changes (Includ- 
ing Viscous Heating) upon Laminar Boundary Layers with 
Variable Free-Stream Velocity,’’ Journal of the Aeronautical 
Sciences, vol. 21, no. 7, July, 1954, pp. 459-474. 

7 Stine, H. A., and Wanlass, K., “Theoretical and Experi- 
mental Investigation of Aerodynamic-Heating and Isothermal 
Heat-Transfer Parameters on a Hemispherical Nose with Lami- 
nar Boundary Layer at Supersonic Mach Numbers,’? NACA 
TN 3344, 1954. 

8 Cohen, C. B., and Reshotko, E., ‘Similar Solutions for the 
Compressible Laminar Boundary Layer with Heat Transfer 
and Pressure Gradient,’’ NACA TN 3325, 1955. 


JET PROPULSION 


1 
1 
4 
‘ 
| 
7 
a 
{ 
| 
~ | 
© 
| 
268 
j 


9 Chapman, D. R., and Rubesin, M. W., “Temperature and 
Velocity Profiles in the Compressible Laminar Boundary Layer 
with Arbitrary Distribution of Surface Temperature,’’ Journal 
of the Aeronautical Sciences, vol. 16, no. 9, Sept. 1949, pp. 547-565. 

10 Goldstein, S., “Modern Developments in Fluid Dynamics,”’ 
Vol. 2, pp. 631-632. 

11 Sibulkin, M. J., “Heat Transfer Near the Forward Stag- 
nation Point of a Body of Revolution,’’ Journal of the Aero- 
nautical Sciences, vol. 19, no. 8, August 1952, pp. 570-671. 

12 Reshotko, E., and Cohen, C. B., “Heat Transfer at the 
Forward Stagnation Point of Blunt Bodies,’’ NACA TN 3513, 
1955. 

13 Beckwith, I. E., “The Effect of Dissociation in the Stagna- 
tion Region of a Blunt-Nosed Body,’’ Journal of the Aeronautical 
Sciences, vol. 20, no. 9, Sept. 1953, pp. 645-646. 

14 Mark, R., “Compressible Laminar Heat Transfer Near 
the Stagnation Point of Blunt Bodies of Revolution,’’ Convair 
Report no. ZA-7-016, April 11, 1955, San Diego, Calif. 

15 Lees, L., ‘Hypersonic Flow,’’ paper presented at Fifth 
International Aeronautical Conference (IAS-R.AeS), Los Ange- 
les, Calif., June 20-24, 1955. I.A.S. Preprint no. 554. 

16 Korobkin, I., “Laminar Heat Transfer Characteristics of a 
Hemisphere for the Mach Number Range 1.9 to 4.9,’’ U.S. Naval 
Ordnance Laboratory, White Oak, Md., NAVORD Report no. 
3841, October 10, 1954. 

17. Margenau, H., “The Forces Between a Hydrogen Mole- 
cule and a Hydrogen Atom,’’ Physical Review, vol. 66, nos. 11 and 
12, December 1 and 15, 1944, pp. 303-306. 

18 Hansen, C. F., “Note on the Prandtl Number for Disso- 
ciated Air,’’ Journal of the Aeronautical Sciences, vol. 20, no. 11, 
November 1953, pp. 789-790. 

19 Kopineck, H.-J., ‘“Quantentheorie des N2-Molekiils,’’ 
Zeitschrift fiir Naturforschung, Teil a, Band 7a, 1952; part I, pp. 
22-23; part II, pp. 314-324. 


(Continued on page 274) 


Effect of Chemical Reactions in the Boundary Layer on 
Convective Heat Transfer 


(Continued from page 258) 


stoichiometric mixture ratio and a total pressure of 20 at- 
mospheres. Inasmuch as recombination is presumed to occur 
at equilibrium, the stoichiometry of the various reactions 
shown in Table 2 was arbitrarily chosen to yield stable end 
products, i.e., Ne and H.O. The total percentage increase in 
heat transfer calculated in Table 2 is 20.7 per cent. The use 
of the approximate simple formula given by Equation [9] em- 
ploying a ‘“‘nondissociated”’ flame temperature T.° of 3555 K 
shows a heat transfer increase of 20.6 per cent, which is in good 
agreement with the more detailed analysis in Table 2. The 
small discrepancy is accounted for in the expansion 
(1 Nim)™ =1+ Nim. 

The actual heat flux increase may be somewhat larger than 
shown in Table 2 since two of the reacting species, H and Hp, 
occur at appreciable mole fractions and have values of (D/x) 
> 1. Hence the chemical kinetics may not be sufficiently fast 
to maintain equilibrium for these species, and it is clear from 
Equation [6] that, if some of these species react where D ~ e, 
the heat flux should be increased somewhat. A detailed in- 
vestigation of this correction can be made by solving Equa- 
tions [1] and [2] with suitable expressions for the chemical 


reaction rates. 


IV Conclusions 


Heat transfer from flames to cooled surfaces generally can be 
expected to show an increase because of the heat release re- 
sulting from the recombination of dissociated species in the 
boundary layers and on the surface. For most chemical 
species the reaction tends to completion in the intermediate 
layer where D and x are small with respect to «. Thus, even 
if D x, (D + + 1; and the chemical reaction 
rates are fast because of the high temperatures encountered 


Aprit 1956 


in this region. Under these conditions, a simplified expres- 
sion may be used for the increase in heat flux which applies 
regardless of the location of the heat release as long as D + 


Table 1 Equilibrium flame composition for representative 


system? 
AF tormation 
Species Mole fraction kcal/mole 
H 0.0100 51.6 
OH 0.0377 10.0 
O 0.0039 58.6 
N 0.0002 112.5 
Hz 0.0582 0 
NO 0.0080 21.5 
Oz 0.0162 0 


« Ammonia (l)-oxygen (g) at 20 atm and stoichiometric 
mixture ratio. 


Table 2 Theoretical increase in heat flux for representa- 
tive system* 
Heat 


Reaction (A +B Mole fraction of release Q, Agq/q’, 


= products) —reactants’—. kcal/mole % 
A B 

H + OH = 

H,0(g) 0.0100 0.0100 119.4 4.2 
OH + '/2H: = 0.0277 0.0139 67.8 6.7 

H,0(g) 0.0039 0.0039 58.6 0.8 
0 + = (g) 
NO + Hz 

H:,0(g) + 1/2N2 0.0080 0.0080 79.3 2.2 
N+N=N,2 0.0002 0.0002 225 0.2 
Hz + 1/202 = 

H.0(g) 0.0324 0.0162 57.8 6.6 

Total 20.7 


* Ammonia (l)-oxygen (g) at 20 atm and stoichiometric 
mixture ratio; flame temperature 3053 K; wall tempera- 
ture 600 K; average heat capacity 11.7 cal mole“! °C™ 
at 2700 K. 

> Cf. Table 1. 


¢™« +. The theoretical analysis considers also the effect 
of surface reaction on heat transfer. In this case, certain 
reactions such as the heterogeneous dissociation of NH; may 
show a decrease in the heat-flux ratio due to the endothermic 
nature of the reaction. 

For the usual case of heat transfer from flames, the concept 
of the “nondissociated” flame temperature T;° is introduced 
since it permits a simple, yet reasonably accurate calculation 
of the change in heat transfer. This concept should be useful 
in predicting high-temperature, heat-transfer coefficients 
from low-temperature measurements. 
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On Laminar Boundary Layers With Heat Transfer 


WALLACE D. HAYES! 
The Ramo- Wooldridge Corporation, Los Angeles, Calif. 


In part I, a transformation is given which reduces the 
laminar boundary layer equations approximately to the 
incompressible form in the case of an imperfect gas. The 
transformation is related to the Stewartson-Illingworth 
transformation, and involves fitting an approximate rela- 
tion to the density profile. In part II, the problem of 
utilizing similar solutions in a momentum-integral 
method is treated with an approach in which the boundary 
layer thickness is characterized by a single function of dis- 
tance along the boundary layer; the momentum, dis- 
placement, and enthalpy-defect thickness are expressed 
in dimensionless forms directly calculable from the similar 
solutions. It is noted that a basic inconsistency appears 
in the energy relation in that two different expressions for 
the heat transfer at the wall are obtained if the classic 
momentum-integral method is used. Various alternate 
methods are discussed. 


Nomenclature 
u,v = velocity components 
z,y = Cartesian coordinates 
= transformed velocity components 
%,9 = transformed coordinates 
p = density 
B = viscosity 
h = enthalpy 
H = total enthalpy 
o = Prandtl number 
a(x) = transformation function 
N(§) = viscosity ratio function 
dimensionless parameters 
Cp = specific heat at constant pressure 
K = coefficient of thermal expansion 
T = absolute temperature 
g(x) = boundary layer thickness function 
B = similarity parameter 
f = reduced stream function 
n = reduced normal coordinate 
A = dimensionless momentum thickness 
é = dimensionless displacement thickness 
A = dimensionless enthalpy-defect thickness 
A = ratio of displacement to momentum thickness 
A* = logarithmic derivative of A with respect to subscript 
A* = logarithmic derivative of A with respect to subscript 
B = enthalpy similarity parameter | 
T(z) = kinetic energy correction parameter 
4, G2 = vorticity parameters for exterior flow 
Subscripts 


1 = free-stream conditions 
w = wall conditions 
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I APPROXIMATE METHOD FOR AN 
IMPERFECT GAS 


Introduction 


NALYTICAL methods for treating laminar boundary 
layers in a compressible fluid are generally based upon 
the assumption that the fluid is a perfect gas with constant 
ratio of specific heats. The first part of this paper presents an 
approximate method for treating the case of a gas which is not 
perfect; this method utilizes a velocity transformation similar 
to the Stewartson-Illingworth transformation (1, 2).2 The 
basic Howarth transformation which replaces normal distance 
by normal mass is included, of course. 
The boundary layer equations are taken in their two-dimen- 
sional form; the modification required for the axisymmetric 
case is indicated later. 


Basic Equations 


The basic boundary-layer equations for a two-dimensional 
laminar boundary layer may be expressed 


d 10 
y p dx p Oy 


ox oy p Oy Lo \ Oy 2 oy 


ox oy 


In these equations the total enthalpy H is given by 


and the viscous shear and conductive heat transfer are 


u(Ou/dy) and (u/c)(Oh/dy), respectively. 
The quantity N is introduced, defined by 


Using this terminology, the quantity p,u, may be taken out- 
side the derivative in the last terms of Equations [I-la] and 
[I-1b]. 

A variant of the Stewartson-Illingworth transformation is 
now introduced, through a function a(x) which is as yet unde- 
termined. A reduced longitudinal velocity component 1 is 
introduced 


in terms of which the basic equations are rewritten 


dz a? dz a dz ap p 
Or =a Oy ap OyLep \ dy 2 oy 
{I-5b] 
Capt pv 
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The quantity « is a function of z, in general. A dimensionless 
yariable @ is introduced to replace H, through a parameter a 
which will be assumed to be constant 


1 


It may be noted that 6, = 1. 
An approximate relation connecting the density ratio and 
the variable @ is now postulated 


2 
p 


From Equations [I-6], [I-7], and [I-8] may be obtained the 
relation 
pi UW du; a? da di, 


The transformation is now completed with the introduction 
of the transformed quantities 


j= af p {I-10b] 
1 
AP ox 


The basic equations now take the form 


Except for the term in (¢ — 1) in Equation [I-11b], these 
equations are now in the low-velocity form for a perfect gas. 


Fitting the Density Ratio 


The application of the foregoing transformation requires the 
establishment of the approximate validity of Equation [I-8]. 
This relation is automatically satisfied under free-stream con- 
ditions. It is now required that it be satisfied at the wall and 
in its first derivative at the wall. 

Using Equation [I-7], the condition at the wall is expressed 


HA, K Pw 


Eliminating « yields 
a {I-15} 


whence «x may be expressed 


It is presumed that a as given by Equation [I-15] is ap- 
proximately constant. If this is not so, a suitable value must 
be chosen and the fitting accomplished by obtaining « from 
Equation [I-12]. 

An important special case is that for which the gas is im- 
perfect in general, but under the conditions at and near the 
wall behaves like a perfect gas with constant ratio of specific 
heats. For such a perfect gas K = T-, and C,/K = h. In 
this case the parameter a@ as derived from Equation [I-15] is 
zero, and 


The case for which the gas is perfect throughout serves as a 
check on these results. In this case x = H,/h, and Equation 
[I-6] may be integrated to give a proportional to h,'*. The 
transformation in this case is identical with that of Stewartson 
(1) and Illingworth (2), for which Equation [I-8] is exact. 


Axisymmetric Case 


For a body of revolution at zero incidence for which the 
body curvatures times boundary-layer thickness are small, the 
well-known Mangler transformation may be used to reduce 
the equations to the two-dimensional form. 

A shape function r(x) appears as a factor within the dif- 
ferentials of the continuity equation Equations [I-Ic], 
[I-5c]. The Mangler transformation may be included in the 
transformation of Equations [I-10] by replacing a by ar, 
Prolty DY Pui”, and v by or in those equations. The set of 
Equations [I-11] then follow from the altered set Equations 
{I-5]. 


Forward Stagnation Point 


The density ratio only plays a part in determining the 
boundary layer solution if the pressure gradient is different 
from zero. To exemplify the method presented here, the case 
of a forward stagnation point is considered. The function g 
and the parameter 8 are defined 


g? 2 fa dk, {I-18a] 
diy 
[I-18b] 


and a stream function gf and a reduced normal coordinate 7 
are introduced such that 


oat [I- 19a] 
oy 
I-19b 
[ ] 


The (o — 1) term of Equation [I-11b] is dropped, 8 and @,, are 
assumed to be constant, and N and o are assumed to be func- 
tions of 7 alone. Under these conditions, similar boundary 
layer profiles may be obtained, corresponding to the solution 
of the ordinary differential equations 


(Nfan)n + + — fy?) =0........ [1-202] 


These are essentially the similar solutions found by Stewartson 
(1) and treated in detail with extensive calculations for the 
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The quantity « is now defined, relating the functions u(x) and Pe — 
= H, — he 
du, a du, 
: 

a=0.............[H7] 
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ae” 
he 1 — a) p 
The coefficient of thermal expression K is introduced a oe 

j 

1/p) 

i, 
in terms of which the condition on the derivative is expressed ae 
1956 = 


case (o = 1, N = 1) by Cohen and Reshotko (3). These solu- 
tions have also been studied by Li and Nagamatsu (4) and by 
Levy (5). Cohen’s quantity S is equal to @ — 1 here. 

In the vicinity of a forward stagnation point the quantity 
x may be assumed constant and the relations 


established locally. The resulting similar profile is one for 
which 


instead of 6 = 1, as in the perfect-gas case. 
In the axisymmetric case the influence of the Mangler trans- 
formation with r = z is to change Equation [I-21b] to 


leading to the value [I-24] 


in place of 8 = 1/2, as in the perfect-gas case. 

In either the two-dimensional or the axisymmetric case, 
then, the fact that the gas is imperfect may be approximately 
taken into account for a forward stagnation point through a 
change in the similarity parameter 8 of the boundary-layer 
profile. 


Shear Stress and Heat Transfer 


In order to reconnect the physical quantities of interest with 
the analysis presented, the shear stress u(Ou/Oy) and the heat 
transfer (u/a)(O0h/Oy) are expressed in terms of the quantities 
introduced. The shear stress is given 


The heat transfer is given 


dy dy ay og 


o Oy oy 


Tw 


(* *) (1 — a)a (1 — 


In the axisymmetric case the factor r appears in these equa- 
tions. 


II MOMENTUM-INTEGRAL METHODS 


Introduction 


The calculation of laminar boundary layers with an arbi- 
trary pressure distribution is generally carried out with one of 
various assumptions as to the profiles which may be present, 
utilizing equations involving integrals taken across the bound- 
ary layer. The primary equation considered is that de- 
rived from the momentum equation which involves the 
momentum and displacement thicknesses. It is convenient to 
assume that the boundary layer profiles at any point are the 
same as those for one of a family of similar solutions for which 
separate calculations are available. This second part of the 
paper presents a modification of the classical momentum- 
integral method in which the various integral thicknesses are 
used in the dimensionless forms which arise naturally from the 
similar solutions. 

Application of similar solutions to calculations with arbi- 
trary pressure distributions was first made by Falkner and 
Skan (6) in a method which did not use the momentum-inte- 
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gral concept. The two concepts were combined by Thwaites 
(7) who assumed a correlation connecting the shear stress at 
the wall, the local pressure gradient, and the ratio of displace- 
ment to momentum thickness, and whose selection of a par- 
ticular correlation was influenced by available calculations of 
similar solutions. The approach of Thwaites has been ex- 
tended by Rott and Crabtree (8) who treated in particular 
the case ¢ = 1, 0, = 1, and by Cohen and Reshotko (9) who 
treated in particular the case o = 1, 6, = const. For detailed 
discussion and extensive references, the reader is referred to 
the authors cited. 

The approach of this note is similar in many respects to the 
Thwaite approach; it differs in that the basic boundary layer 
thickness used as a reference value is a function g(x) appearing 
in the separation-of-variables procedure underlying the similar 
solution analysis rather than the related quantity of Thwaites 
defined in terms of the momentum thickness. 


Basic Equations 


The laminar boundary layer equations will be considered in 
their fully reduced form, as given in Equations [I-11]. The 
bars have been dropped here, so that the fact that wu, v, x, and 
y are reduced quantities must be kept in mind. 


ou ov 
= by {II-1c] 
II-2 


is a parameter equal to twice the ratio of kinetic energy to 
total enthalpy in the free stream. The boundary conditions 
are 

y=0: u=v=0, 6=80,(z), (N = 1).. [II-3a] 


y= 


A function g(x) which is as yet undetermined is introduced, 
together with a stream function for Equation [II—1c] equal to 
gf. The variable y is replaced by a variable n defined by 


[11-4] 
g 
With these quantities Equations [II-1] become 
9'9 g?u;' 
(Nfon)n mom + (6 — fy?) (fanz Sanfz) 
[II-5a] 


N 
n 1 1 


{II-5b] 
under the boundary conditions 
n = 0: f =090, fn = 9, 6=0,(z), (N = 1) 
[II-6a] 


An integral of Equations [II-5] from 7 = 0 to © isnow taken, 
giving 
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where f, (IT-8a] 
(0 — fa) dy = [II-8b] 


are dimensionless momentum, displacement, and enthalpy- 
defect thicknesses, respectively. Equations [II-7] may be re- 
written 


+ + (1 + A)? A = - [II-9a] 
uw dz 


Similar Solutions 


A solution to the set of equations 
(Nfonn + ffan + B(0 — fn?) = 0....... (II-10a] 


+ fey + = 0... {II-10b] 


considered as ordinary equations in 7 is here defined as a 
similar solution. Such a solution corresponds to a physically 
realizable similar solution which is an x-independent solution 
to the boundary layer equations only if N and ¢ are functions 
of y alone and I is constant or ¢ = 1. It is now assumed that 
the variables A, A, andf,,, in Equation [II-9a] and A and 6,, in 
Equation [II-9b] are the same as they would be if the actual 
solution were a similar one with the local values of 6,, and T 
and with a suitable value of 8. This is the underlying assump- 
tion of this method and corresponds to Thwaites’ correlation 
assumption. With this assumption, the relations 


(1 + (1 + A)BJA = fon,.......... {II-1la] 


1 


are valid, and Equations [II-9] may be put in the form 


g da 


.. {II-12a] 
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The quantities A and A are now considered to be functions of 
8 and @, alone, for the purpose of simplifying the analysis. 
Dependence upon I or o,, may be treated in the same way as 
dependence upon @,,, while variations in N(y) or in ¢/¢, with x 
would entail some additional complication. Logarithmic par- 
tial derivatives are defined 


B dA 6, OA 
A* —; At, ="—....... 
B OA OA 
*,=-—; A*, =——....... - 


with respect to 6 and 6,; similar notation may be used with 
respect to T' and o,. Equations [II-12] now becomes 


B 6 
1 
A )+ 1 = 0... 
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The functions u(x) and 6,(x) are known functions, while g(x) 
and A(z) are as yet unknown. Some latitude is available in the 
manner in which g and £ are determined. 


The Classical Method 


In the classical method of applying the momentum-integral 
method, the enthalpy equation is ignored and an assumption 
is made regarding the relation between the basic parameter 
and the velocity gradient. As interpreted here, Equation 
[II-14b] is ignored and 8 is determined by the relation 


(II-15] 


Equation [II-14a] gives the differential equation for g 


, ” 2 
UW 2 Uy w 


Pw 
{II-16] 
or the equivalent differential equation for 8 
{II-17] 


Either of these equations may be integrated numerically to 
obtain the quantities g and £. 

The heat transfer may be calculated from the momentum 
solution in two different ways. The relation Equation [II-11b] 
may be used, corresponding to the slope of the assumed 
local enthalpy profile corresponding to the velocity profile; 
or, the relation Equation [II-9b] may be used, corresponding 
to the heat transfer necessary to maintain the assumed 
corresponding enthalpy profile. These give the heat transfer 
equations 


g? 2 
Ony = (A*g — A*s) — — 
1 Yw 
.. (II-18b] 


These equations cannot be consistent in general unless A 
equals a constant times A over the entire range of the parame- 
ters. It is not completely clear which of Equations [II-18] 
should be preferable, although it appears to the writer that the 
simpler form Equation [II-18a] should be more trustworthy. 

A basic inconsistency in the calculation of heat transfer from 
the enthalpy profile corresponding to the assumed velocity 
profile thus appears. That this is so is not surprising in light 
of the fact that the upstream temperature profile has no effect 
on the heat transfer calculated by this method, except in so far 
as it may affect the velocity profile. 


Alternate Methods 


A. In this method the relation Equation [II-15] connect- 
ing 6 and g is dropped, and Equations [II—14] are to be solved 
as a pair of simultaneous equations for determining both 6 and 
g. The derivatives of 8 and g may be separated, giving the 


(A*s — A*s) + 
(A%, — + (1 + A) a) =0 
- 1) — A*s) + 
(A*pA*,, — ) — (1 + ay )a 
.. [1I-19b] 
273 


: 
tes 
at } 
of 
lar 
ho 
ed 
to 
he 2,.? 
er 
| uj 
ow 
) 
| 


This set of equations may pose difficulties of a mathematical 
nature, particularly with regard to the question of stability. 
An approximate statement of a stability condition may be 
formulated, that 


This method gives a simultaneous solution for the skin friction 
and the wall heat transfer, in a self-consistent form, provided 
that a mathematically unobjectionable solution can be found. 

B. The classical method as just described may be followed 
except for the computation of heat transfer properties. The 
enthalpy profile is now assumed to be a similar one with the 
appropriate values of 6, and I’ with the function g(x) com- 
puted for the momentum equation, but with a different value 
of the fundamental parameter B, corresponding to a different 
value of A. This gives the equation for fin place of Equation 


{II-21 


The heat transfer is given directly in terms of A. This method 
may be expected to give reasonably trustworthy results if the 
difference between 6 and 8 does not become too great. 

C. Fora very rough estimate all the derivatives of A and 
A in Equations [II-14] may be neglected, giving the relation 
Equation [II-15] for 8 in terms of g, and the relation 


Although rather crude, this method should give an estimate of 
the heat transfer from Equation [II-11b] which is about as 
accurate as that given by the classical method. 

D. The solution obtained by method C may be used to 
calculate the A and A derivative terms in Equations [II-14] 
for use in an algebraic calculation of second approximation. 
Such a successive approximation scheme, although not 
strictly convergent, may have semiconvergent properties that 
make it useful if |A*s — A*g| is small. This should give es- 
sentially the solution of method A by an easier method. 


External Rotational Flow 


One of the basic boundary layer assumptions is that the 
layer is so thin that the structure of the external flow does not 
affect the structure of the boundary layer except through the 
single function u;(x). It is possible that in certain cases, if 
the vorticity in the external flow is sufficiently large, that this 
vorticity may influence the boundary layer structure even 
though the layer is quite thin. Fora similar boundary layer 
profile similarity conditions must be imposed on the external 
vorticity which are, in general, inconsistent with the steady- 
flow conditions on the inviscid exterior flow. However, the 
case of the neighborhood of the stagnation point on a blunt 
body of revolution is one in which finite vorticity at the wall 
does follow the similarity condition. With the type of use of 
similar profiles of this note the requirement of the similarity 
condition on the vorticity is not essential. 

The inviscid profile (no boundary layer present) is as- 
sumed to be of the form 


fn = 1 + ain + aon? + 


The procedure to be followed is first to form the inviscid ex- 
pression for f 


f=nt+ aun? + + 


and then eliminate » between Equations {II-26] and 
[II-27]. The resultant relation connecting f, and f then re- 
places the boundary condition f, = 1 of Equation [II-6b]. 
With only the quantity a; present or important, this relation 
is 
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The absolute magnitude of the parameter a; serves as a 
measure of its importance in affecting the boundary layer 
profile. Unless a; is negligibly small, it should be added as a 
parameter determining the boundary layer solution. 
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Technical Notes 


Sun-Follower for High Altitude 
Sounding Rocket 


D. D. TERWILLIGER! and G. J. GRANROS? 
Aircraft Armaments, Inc., Cockeysville, Md. 


PROGRAM of upper atmosphere research has required 
the direct measurement of solar radiation below 3100 
Angstroms and atmospheric absorption of this radiation. 
Units have been designed, built, and flown which house a 
spectrograph and direct the solar radiation into this spectro- 
graph. These units, known as Sun-Followers, are mounted 
on high altitude sounding rockets. In order to obtain 
sufficient exposure times at the shorter wave lengths, sunlight 
is directed into the spectrograph regardless of the spin and 
precession of the rocket. In the latest equipment this is 
accomplished by means of a two-axis tracker, design and 
performance of which are described in this report. 

The present design of the Sun-Follower has been developed 
as a result of the experience gained on several previous units 
built and tested. The first units were single axis units, track- 
ing the sun in azimuth only, and were mounted on a V-2 
and on an Aerobee. A subsequent pair of units was de- 
signed to mount on a Viking rocket and an Aerobee. These 
units were also single axis types but incorporated design 
improvements. 

In order to obtain more intense spectra the latest units 
track the sun in both axes. This is accomplished by two- 
axis control of a mirror which reflects the sunlight down the 
longitudinal axis of the rocket, on which the spectrograph 
is fixed. An azimuth drive corrects for rocket roll motion 
and an elevation drive corrects for rocket pitch and yaw 
motion. Photoelectric cells mounted on the line of sight are 
used to sense angular pointing errors and, by means of servo 
systems, the error is kept to a minimum. 

Figs. 1, 2, and 3 show some of the construction details of the 
latest Sun-Followers. A block diagram of the control system 
is shown in Fig. 4. 


Photoelectric Cells 


For control of elevation and azimuth axes, two pairs of 
photoelectric cells are employed; one pair being used for 
coarse control to search for the sun, lock on, and to bring the 
Sun-Follower approximately in line with the sun, the second 
pair of photoelectric cells to keep the Sun-Follower accurately 
on the sun to within +0.25 deg. A special switching circuit 
is used to remove the supply voltage to the coarse photo- 
electric cells when the Sun-Follower is operating on the fine 
eye system. 

The coarse eyes for azimuth control are mounted on the 
azimuth ring as shown in Fig. 3. The fine azimuth and 
fine and coarse elevation eyes are mounted in an eye block 
which is driven by the elevation servo to the line of sight 
(Fig. 3). The mirror is geared directly to this eye block but 
is driven at one-half speed in order to obtain the reflection on 
the longitudinal axis. 

Each photoelectric cell is followed by a cathode follower, 
and the pointing errors are obtained as the differential out- 
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Fig. 1 Assembly of Sun-Follower including 
rocket skin sections 
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Fig. 2. Internal details of Sun-Follower 


Eprror’s Nore: This section of Jer Proputsion is open to short manuscripts describing new developments or offering comments on 
inte previously published. Such manuscripts are published without editorial review, usually within two months of the date of receipt. 


uirements as to style are the same as for regular contributions (see first page of this issue). 


Aprin 1956 


aS 8 
om- 
: 
vol, 
ure 
er CONTAINING 
| AND SWITCH 
TUBES 
EYE BLOCKS a 
al 
COARSE 
EYE 
‘ 
i 
y FOR AZIMUTH . 
of 


Fig. 3 Details of Sun-Follower aD eyes 


puts of ‘the pairs of cathode followers. The coarse photo- 
electric cells use an optical system to produce a pickup 
pattern which is essentially a half sinusoid with the maximum 
dead ahead, as shown in Fig. 5a. The two coarse eyes in 
azimuth are mounted with their optical axes located +75 deg 
from the boresight line. In this manner a composite error 
curve as shown in Fig. 5b is obtained. 

This permits a proper servo-error signal to be obtained 
from +165 deg to —165 deg. The remaining 30 deg directly 
opposite from the sun is not covered, but unbalance of the 
servo system, stray light, or slow roll of the rocket should 
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Fig. 5a, b Coarse eye assembly and pickup patterns 


bring the system out of this possible dead zone. The eleva- 
tion coarse eyes are boresighted along the line-of-sight. 
A blinder for each eye is provided to distort the pickup 
pattern to give a large signal on either side of dead ahead. 
This results in a relatively broad null but the fine eyes are 
capable of controlling the Sun-Follower in this region. As 
elevation coverage is only required to be +45 deg, no dead 
zone exists. 

The fine photoelectric cells use an ed system to obtain 
a pickup pattern as shown in Fig. 6. The design of the 
fine and coarse eye assemblies and the eyes used on this model 
Sun-Follower were supplied by the University of Colorado. 
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Fig. 4 Block diagram of control circuitry 


JET PROPULSION 


= 
TELEMETRY, a 
ke 
— a 
i 
a fi 
e 
il 
e 
— al 
4 CF / 
| 
a 
| k 
LE 
8 PHOTO-CELLS 
NETWORKS 
1P42 sh 
8) 
| fal 
4 CF \ al 
{ CF 4 $ er 
68F7. 5702. ch 
276 A 


These eyes were developed by them for use on their Biaxial 
Pointing Control System (Special Report No. 7, “Photo- 
electric Eyes for a Solar Pointing Control”). These are 
arranged so that one photoelectric cell has a pattern as shown 
by the solid line, to the right of center, and the other photo- 
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Fig. 6 Fine eye assembly and pickup pattern 


electric cell of the pair has the mirror image on the left of 
center as shown by the dashed line. Similarly, for elevation, 
the fine eye pair is arranged above and below center. The 
combination of lens and knife edge at the focal point of the 
lens gives the sharp cutoff on the unwanted side of center, 
and the steep linear rise on the other. 


Coarse to Fine Switching Circuitry 


A special switching circuit was developed to perform two 
important functions. The first was to switch from coarse to 
fine eye control. The second was to disable the coarse 
eyes once the fine eyes were functioning. This is required 
in order to prevent coercion of the fine data during tracking. 

The possibility of coercion exists because the coarse azimuth 
eyes are mounted on the azimuth axis while the fine eyes are 
mounted on the line of sight or deflection axis. Any mis- 
alignment of the coarse eyes as a function of elevation will 
cause a coercive error. Also extraneous sources of light 
such as reflections from the earth or clouds may cause errors 
in the coarse system. In one of the previous rocket firings, 
the telemetered data indicated that this actually happened. 

In the present switching circuitry, the fine eye photoelectric 
cell cathode follower outputs go to the grids of a dual triode, 
the plates of which are tied together. If either of the fine 
eye photoelectric cells are illuminated, the grids of the switch 
tube are driven sufficiently positive to saturate the tube and 
to bring its plate potential to a low value. The coarse eye 
photoelectric cell anodes are connected to the switch tube 
plates, so that when the fine eyes are illuminated, the coarse 
eye photoelectric cell anode potential drops to a sufficiently 
low value that the coarse eyes are inoperative regardless of 
the illumination on these cells. 

The outputs of the coarse and fine systems are mixed as 
shown in Fig. 4, by means of adding resistors, with only one 
system operative at any one time. 


Servo Amplifier 


Following the mixing network, the differential signals are 
amplified, passed through a lead network for stabilization 
of the servo loop, and modulated by means of a 400 cps 
chopper. In the elevation system, this a-c error signal is 
amplified to drive the elevation servo motor which in turn 
drives the elevation axis. In the azimuth system, the a-c 
error signal is amplified, demodulated by a synchronous 
chopper, and fed to the grids of pentode tubes which control 
the current in the clutch coils of a Lear Magnadrive. This 
Magnadrive in turn drives the azimuth axis of the Sun- 
Follower. 
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Additional stabilization for the azimuth axis was found 
to be necessary due to the high azimuth inertia and ex- 
tremely high lock-on rates. This was accomplished by in- 
serting a small condenser across the coarse eye adding re- 
sistors, to provide a small additional amount of lead. In 
this manner it was possible to make the system capable of 
searching and locking on the sun while traveling at azimuth 
rates in excess of 720 deg/sec, and following with an accuracy of 
less than +0.25 deg. No additional stabilization was re- 
quired for the elevation servo system due its relatively small 
inertia. Elevation tracking accuracy also is +0.25 deg. 

The 400 cps power required for the synchronous choppers, 
and the fixed phase of the 2-phase induction motor in the 
elevation servo were obtained from a separate oscillator and 
power amplifier. Anode voltages for the circuitry were ob- 
tained from a dynamotor and vacuum tube voltage regulator 
circuit. 

In order to permit conventional electronic circuitry to be 
used on equipment exposed to very high altitudes, the entire 
electronic package is placed in a pressurized container. No 
special precautions are taken to insure that the container is 
airtight. It is ground tested at 15 psig and sealed so that 
leakage does not exceed 1 psi per 10 min. As the maximum 
flight times are less than 5 min, this is considered adequate. 

All circuits not pressurized are limited to 27 volts on any 


exposed contacts. 


Telemetry 


In order to’ obtain a continuous check of the operation of 
the circuitry, various data are telemetered to ground receivers 
in the vicinity of the rocket launching site. These data are: 
Plate voltage, 26.5 volt battery supply voltage, Sun-Follower 
azimuth and elevation pointing angles, azimuth and eleva- 
tion pointing error signals, azimuth magnadrive clutch cur- 
rents. An FM-FM telemetering system using a single sub- 
carrier oscillator together with a commutator is used to trans- 
mit the required data. 


Results of Flight Tests to Date 


The history of the previous four Sun-Follower flights show 
a steady increase in over-all performance and in data gathered 
by the spectrograph. 

The first flight started out with a bang; that is, the V-2 
rocket on which the instruments were mounted blew up and 
was destroyed completely at take-off. 

An Aerobee was used to carry the second Sun-Follower. 
The telemetry records for this flight indicated that the single 
azimuth sensing eye was pointing away from the sun, al- 
though the servo system was nulled. No data were ob- 
tained on the spectrograph camera. 

An Aerobee was used again to carry the third Sun-Follower. 
Usable data were obtained on the spectrograph camera, al- 
though the servo system was oscillating about the null 
position. The oscillation apparently reduced the spectro- 
graph exposure time. 

A Viking rocket was used to carry the fourth Sun-Follower. 
Telemetry records indicated that at 120 sec after take-off 
the rocket began to roll and accelerated to 180 rpm at ap- 
proximately 140 sec after take-off. This roll rate appreciably 
exceeded the Sun-Follower design roll rate of 120 rpm. At 
200 sec the Sun-Follower finally lost track of the sun because 
of the high roll rate. There was little usable spectrographic 
data before 120 sec, due to the atmosphere; however, from 
120 to 200 sec the Sun-Follower tracked sufficiently well to ob- 
tain the best data of the series so far. 

The fifth and sixth Sun-Followers have been built, tested, 
and delivered by Aircraft Armaments Inc., to the Naval 
Research Laboratory where they are undergoing installation 
and check-out tests. The present schedule calls for the fifth 
and sixth Sun-Follower to be fired on an Aerobee rocket at 
White Sands Proving Ground about May 22, 1956. 
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Heat-Up Time of Wire Glow Plugs 


S. K. CLARK! 
University of Michigan, Ann Arbor, Mich. 


Nomenclature 
power input, units per unit time, based on the average re- 
sistance of the wire over the temperature range to be 
encountered 
specific heat of glow plug material < weight of glow plug 
absolute temperature of glow plug wire at any time 
time 
emissivity X surface area of glow plug wire X Stefan- 
Boltzmann constant 
absolute temperature of chamber walls to which glow 
plug loses heat by radiation 
absolute temperature of plug wire at start of heating, 
ie,t=0 
P + BT 
(C/B)'/* 
4By%t/A = 4Ct/Ay 
HE time-temperature relationship in a wire resistance 
glow plug, used as a heat source starter in a reaction 
chamber, can be obtained with a considerable degree of ac- 
curacy by assuming that the primary heat loss of the wire is 
due to radiation effects. This assumption is justified in 
many respects since the chamber walls are usually quite close 
to the wire and almost completely surround it, and since the 
air remaining inside the system under actual starting con- 
ditions is usually stagnant or of negligible velocity. Knowl- 
edge of the time-temperature relationship allows the fuel 
injection to take place at the proper prescheduled time (tem- 
perature) after the power input to the glow plug has begun. 
Alternatively, if the time delay of fuel injection is fixed, solu- 
tion to the problem enables the proper glow plug to be chosen 
for a specified set of conditions, such as a desired wire tem- 
perature at the time of fuel injection. 
Neglecting convection and conduction, the differential 
equation governing temperature is 


where 7 is assumed to remain a constant since the chamber 
mass is almost invariably many times larger than that of the 
glow plug. Equation [1] may also be written 


A\ dT, 


The solution of this differential equation can be obtained by 
direct integration of Equation [la], written in the form 


t 4) dt, 
dt = — (= 
J, rm (To — 
by noting that 
by expansion by partial fractions. The result is 


A (To — + 
t= tan-!— — tan — 


This may be conveniently reduced to dimensionless form by 
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noting that the constant + is the equilibrium or steady state 
temperature of the wire, since from Equation [1] 


= [P/B + Tot]'* = 7 
dt 


Defining the dimensionless ratios 


and assuming further that JT! = To, or that the chamber and 
plug are at the same temperature at the instant resistance 
heating of the plug begins, Equation [3] becomes 
* (6 — 1(A + 1) 

(8 — 1) 
where the limits on the variables \ and 8 are 


O0<A<1 
ASB<1 


a = 2(tan-'g — A) — 


Equation [6] may be shown clearly by a plot of a vs. 8 for 
various fixed values of \. This is illustrated in Fig. 1. It is, 
of course, possible to solve this problem by numerical step-by- 
step integration, in which case the actual change in resistance 
with temperature of the glow plug wire may be used, again 
on a step-by-step basis. Comparison of results obtained 
by the use of Equation [3], using the average resistance of the 
wire over the temperature range encountered, with step-by- 
step calculations shows agreement to within 2 per cent on time 
to rise to a given temperature. 
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As an example of the use of Fig. 1, assume the following 
specific wire plug: 

Calculated equilibrium temperature = 1800 F 

Fixed starting temperature = JT, = T! = 190 F 

Calculated 4C/Ay = 0.5 

To find: Time for wire to reach 1650 F 


Calculate: 
= 650/2260 = 0.287 
= 2110/2260 = 0.934 
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[5) 


6) 


By interpolation, follow along the \ = 0.287 curve until 
reaching 6 = 0.934, at which point a ~ 3.8. Then 
4C 


t = 7.6 sec 


The largest uncertainty in the use of Fig. 1 or Equation [3] 
lies in the choice of the effective emissivity. In regard to this 
it should be pointed out that a single test in which all other 
constants are known except emissivity, and in which wire 
temperature is measured at a specific time before steady-state 
conditions are reached and again after steady-state conditions 
are reached, gives sufficient information to determine the 
average effective emissivity, since 


(6 — + 2] 


A 
B 1g — tan~d) n (@4+ 1a 1) 


Comments on “Manifolds for Solid 
Propellant Rocket Motors” 


HENRI VERDIER? 


Société d’Etude de la Propulsion par Réaction, Villejuif 
(Seine), Paris, France 


Introduction 


N HIS paper “Manifolds for Solid Propellant Rocket 
Motors,” Mr. Geckler gives a very interesting theory for 

computing the adequate size of manifold connecting two 
similar solid propellant rockets. However, the author says 
he does not know of any experimental tests verifying the 
validity of his equations. 

Many static and flight tests have been conducted by 
SEPR with manifolds for double-base solid propellant boosters. 
In static tests, the mechanical and thermal behavior of the 
manifold was experimented by using nozzles with different 
throat areas in each rocket motor. Thus, an important mass 
flow rate of propellant gas was obtained in the manifold. 
From the results of those tests, experimental values of the 
factor CY used in Mr. Geckler’s theory can be deduced as 
stated below. 


Method for Determining CY 


For two rocket motors having the same nozzle throat area 
A,, interconnected by a manifold of cross-sectional area 


Ay, Mr. Geckler writes 
T(1 — n) (Ap; — Apy) 

fl] 
2CY -V/ 2px Ap, 


+1 
2(y — 1) 


Am = At 


y¥+1 


r= 


y = Cp/Cv = ratio of specific heats 

n = pressure exponent in the burning rate rule 

cY = factor in the mass flow rate rule for the 
manifold 

Ap; = pu — px = pressure difference between the two 


chambers, without manifold 
Ap = py — Po = pressure difference between the two 
chambers, with manifold 


For two rocket motors having different nozzle throat 
Received Feb. 28, 1956. 
| “Manifolds for Solid Propellant Rocket Motors,’”’ by Richard 


D. Geckler, Jer Proputsion, vol. 25, Oct. 1955, pp. 540-541. 
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areas, An and Aw we can substitute for Equation [1] the 
following equation 


T(1 — n)((Aupu — — (Aapy — Avpy)] 
Ay = [2] 
2CY-V 2p py — pr) 
The factor CY is given by 
V pulpy — Px) 


The ballistic properties of the double base solid propellant 
used are known 


y = 1.26 r = 0.66 

n = 0.68 1 —n = 0.32 
where A, = propellant burning area 


6 = propellant temperature 
For each test, the different terms of Equation [3] are known 
— n)/2V/2 = 0.074 


An, Am, An, An, and Ay are given; 6 is measured; pi; and 
poi are deduced from Equation [5]; pis and poy are measured; 
a and a, are deduced from Equation [4]; we can then obtain 
numerical values for CY. 


Experimental Results 


Experimental results were obtained for values of: Ap,/ 
poy = (pis — Pos)/Par Varying between 1.5 per cent and 6 per 
cent. 

The precision on numerical values of CY is small because 
(pir — Poy) is not directly measured, but obtained as the 
difference between pi, and po, with an important incertitude. 
Nevertheless, all the results with the type of manifold rep- 
resented on Fig. 1 have given values of CY between 0.38 and 
0.59 

0.38 < CY < 0.59 
CY =0.48+ 0.11 


Moreover, static tests have shown other advantages of 
manifolds than pressure equilibrium. Particularly, the simul- 
taneity of ignition is very good (difference inferior to 0.01 s) 
and so is the simultaneity of extinction. 


Fig. 1 Type of manifold used on SEPR boosters 
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Measurement of Total Emissivity of 
Porous Materials in Use for 
Transpiration Cooling 


E. R. G. ECKERT,! J. P. HARTNETT,’ and 
T. F. IRVINE, JR.® 


University of Minnesota, Minneapolis, Minn. 


Introduction 


RANSPIRATION or porous-wall cooling is one of the 

more promising proposals for the maintainence of toler- 
able temperatures on surfaces subjected to intensive heating. 
The use of such transpiration cooling has been advocated for 
cooling combustion chamber walls, gas turbine blades, and 
external surfaces of high-speed aircraft subjected to aero- 
dynamic heating. The equilibrium temperature of such 
transpiration-cooled surfaces is determined by a balance 
between the convective, conductive, and radiative heat 
transfer. In view of this, the knowledge of the radiative 
characteristics of porous surfaces is of considerable engineer- 
ing importance. Such a study is currently under way at the 
Heat Transfer Laboratory of the University of Minnesota. 
The purpose of this note is to present a description of the test 
apparatus that is being used for the determination of the 
total normal emissivity values of a number of porous ma- 
terials and to discuss its theory. In addition, total emissivity 
data are presented for a sintered bronze material. A more 
comprehensive report covering a number of porous materials 
will be presented at a later date. 


Description of Apparatus 


The test arrangement used in the present emissivity 
measurements is essentially the same as that first described 
by E. Schmidt (1, 2).4 It is illustrated in Figs. 1 and 2. 
The principal piece of equipment is the radiometer which is so 
designed that the radiant energy entering the opening (a) is 
collected by a 5-in-diam gold surface mirror (b) and directed 
to the receiving surface of a thermopile (c) (Kipp and Zonen 
Microthermopile, E-3). The amount of radiant energy 
striking the thermopile is controlled by two orifices located in 
front of the thermopile. The deflection of a sensitive gal- 
vanometer attached to the thermopile gives a measure of the 
energy impingent on the thermopile surface. A black body 
or ideal radiator (d) constructed of copper and heated elec- 
trically serves as a reference for the radiation measurements. 
The porous test sample (e) is held in a balsa wood frame, 
balsa being used for its low thermal conductivity of 0.03 
Btu/hr ft°F. The porous specimen is held tightly against an 
electrically heated copper plate. Both the black body and 
the test sample temperatures are controlled by varying the 
electrical power input to appropriate heaters. A double 
walled container (f) is arranged in front of the test sample. 
Its inner and outer surfaces are painted bjack with a paint of 
high absorptivity for long wave radiation and these surfaces 
are kept at a constant temperature by a stream of water flow- 
ing through the space between the double walls. Its purpose 
is to create a surrounding of well-defined temperature and to 
avoid any reflection of stray radiation. Its cylindrical open- 
ing permits the radiometer to see the test sample and an in- 
dentation of 1!/, in. diam and 4/2 in. long in the surrounding 
surface creates a black body (g) at the surrounding guard 
temperature. 

In using this equipment the thermopile is first removed and 
is replaced by a light source of the same size. The radi- 
ometer may then be directed at the black body (d) by looking 
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REMOVABLE 
OPENING (a) 


GOLD 
SURFACE SIGHTING 
MIRROR (b) 


Fig. 1 Emissivity apparatus 


BLACK BODY -(d) 


x - CALIBRATED 
THERMOCOUPLES 


RADIOMETER 


Fig. 2 Schematic of emissivity apparatus 


through a sighting hole in the back of the radiometer. The 
unit is in proper adjustment when the light bundle, after 
reflection from the gold mirror, goes on through the opening in 
the black body without hitting the adjacent front walls. 
Similar checks are made on the sample and on the guard 
black body (g). The light is then removed and replaced by 
the thermopile. Then the galvanometer deflections oc- 
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mrring when the radiometer points toward the porous sur- 
face (e) and the two black bodies (d) and (g) are recorded. 
These readings, in conjunction with the corresponding tem- 
perature measurements, allow the calculation of the emissivity 
of the sample. 


Analysis of Data 


The radiative energy flux leaving the black body at tem- 
perature 7’, and arriving on a unit area of the thermopile is 
determined by the solid angle intercepted by the two dia- 
phragms ahead of the thermopile and may be designated as 


Fis an angle factor determined by the geometry and position 
of the two diaphragms and ¢ is Boltzmann’s constant. 

Since the sample at temperature 7, is completely sur- 
rounded by a black container at a temperature 7',, the energy 
flux G. leaving its surface and arriving on a unit area of the 
thermopile is composed of the emitted radiation and the re- 
flected radiation 


with e, indicating the emissivity and p, the reflectivity of the 
test specimen. 

With the assumption e, = 1 — p,, the following result is ob- 
tained from Equation [2] 


— 7,9 — [2a] 


The radiant energy flux leaving the small black body in the 
guard surface at temperature 7’, and arriving on a unit area of 
the thermopile is 


Subtracting Equation [3] from Equations [1] and [2a] and 
rearranging results in the following expression for the emis- 


sivity 


The assumption «, = 1 — p, requires that the sample is 
sufficiently thick so that none of the impinging radiation 
leaves it through the back surface and strictly that the emis- 
sivity does not depend on wavelength. Equation [2a] is a 
good approximation even when the last condition is not ful- 
filled provided 7’, is considerably smaller than 7, or when 
both temperatures are nearly equal. The possible error on 
the end result of a variation of the emissivity with wave- 
lengths has to be decided from case to case. 

It now remains to relate the radiation fluxes, G, to the 
thermopile deflection, A. A heat balance may be written on 
the front face of the thermopile for the condition that it is 
irradiated by energy coming from the black body 


where 


loss factor accounting for absorption of energy by CO: 
and water vapor in air and by the gold surface mirror 
= emissivity of thermopile 
T;, = temperature of exposed surface of thermopile 
= radiant energy, absorbed by thermopile, originating 
from surrounding walls 
Gioss = losses from thermopiles by convection and conduction 


° 
It 


We may write similar expressions when the thermopile is 
subjected to radiation from the sample, and from the guard 
black body. 


= oerTy,4 — qrad + [7] 
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If the assumption is made that the decrease in intensity is 
equal in all cases (C; = C2 = C3), and that the environment 
radiation to the thermopile and the conduction and convection 
losses remain the same, we may write, on combining Equations 
[5], [6], and [7] 


G: — G; 47 ;,(T;, T;,) 


The error introduced in Equation [8] by writing it in the 
form on the right is quite small, with the numerator being in 
error approximately by the amount 1.5(7', — T,,)/T,;, and 
the denominator by 1.5(7',, — T,;,)/T;,. The value — T's, 
is of the order of 1 to 2 degrees, thereby giving an error of the 
order of !/2 per cent in the numerator and denominator. The 
error in the quotient is probably less than this amount. 

The assumption C, = C, = C; may not be admissible in all 
cases since the absorption of carbon dioxide and water vapor 
depends strongly on wave length. Provision is, therefore, 
made to flush the radiometer with a nonabsorbing gas. 

The thermopile indication A is directly proportional to the 
differential between the front side temperature, 7',, of the 
thermopile and the unexposed back side temperature, T'*. 

Therefore, when the radiometer is sighted on the black body, 
the indication is 


The thermopile is so constructed that 7'* remains constant 
when the radiometer is directed toward the other surfaces, 
yielding the following relations 


Substituting Equations [9], [10], and [11] in Equation 
[8], and combining with Equation [4], we arrive at the final 
relationship for 


— T,*\ [Ae — As 


Thus the knowledge of the three temperatures and the 
deflections A,-A; and A,—A; allows the calculation of the 
emissivity of the porous sample. 


Total Emissivity Values for Porous Bronze 


Data were obtained for a sintered bronze material currently 
used for transpiration cooling experiments. The test sample 
of 4 in. sq was removed from a large plate 14 in. wide and 
72 in. long X 1/2 in. thick. Its porosity is such that an 8 psi 
differential across the plate allows an air flow of 0.55 lb per 
sq ft per sec. It was fabricated by Powdercraft, of Spartan- 
burg, South Carolina. The surface roughness of the plate 
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Fig. 3 Total normal emissivity of porous bronze 
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was found to have an absolute average roughness value of 
0.0005 in. when it was measured with a Brush Surfindicatoy 
having a stylus with a tip radius of 0.0005 in. 

The test results are presented in Fig. 3 and cover a tempera- 
ture range of 160 F to 350 F. The emissivity value found for 
this material is 0.57 + 0.04. At 350 F the sintered bronze 
rapidly oxidized which resylted in an increase in the emissivity 
value to 0.63. 

The emissivity value of the solid bronze surface obtained by 
grinding the porous test specimen to close the openings was 
measured and found to be 0.04 at 150 F. This is an order of 
magnitude lower than the value for the porous surface and 
emphasizes the possible error involved if emissivity values for 
solid surfaces are used to predict the radiative heat transfer 
from porous materials. 
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The Near-Constancy of Full-Power 
Duration for Unboosted Rocket Vehicles 


WILLARD P. BERGGREN! 
University of Bridgeport, Bridgeport, Conn. 


HE firing duration (at constant thrust) of a rocket engine 

which constitutes the sole propulsion system of a piloted 
craft or single-stage missile may conveniently be expressed 
in terms of the initial thrust-weight ratio F/Wo, the mass ratio 
¢, and the specific impulse /,, 


T =1,,{ Wo/F 


One example of magnitydes, for a piloted aircraft, might be 
I.p = 220 sec, § = 0.5, F/Wo = 1.0, for which combination 
T = 110sec. Thus, a rocket-propelled fighter of moderately 
high performance might operate less than two minutes at 
full throttle. 

A second typical set of magnitudes for insertion into Equa- 
tion [1] would apply to a missile of advanced design where 

" : ee I, = 250, § = 0.9, F/Wo = 3.0; giving a duration T = 75 
POWERING: —— sec. This vehicle would, if launched vertically, reach a speed 
ubsonic Supersonic of the order of Mach 2.5 at 15,000-ft altitude and an accelera- 
—, —— tion of 30 g’s at burnout. These conditions might be in- 
Guided Missiles : tolerable because of drag and stress, respectively, and a 25 
Dcicoblens per cent reduction of design thrust (F/Wo = 2.25) would 
P result in 7’ = 100 sec. 

Convertiplanes It is rather striking to find this latter figure so near to the 
Thrust Augmenters — full-throttle duration of the piloted aircraft considered in the 
Target Drones first example. This near-constancy of full-powered operating 
JATO 4 time for unboosted vehicles of widely varying designs has 
been recognized for some time; however, the elementary 
oe Equation [1] may be an aid to quick visualization of the in- 

= fluence of basic design parameters in this regard. 
Designers & Manufacturers ’ Power plant design questions of a most fundamental sort, 
of Fuel Injection Nozzles a such as pressurization and cooling, depend upon firing dura- 
; tion. Because of the near-constancy of duration often found, 


D E LA Vv A N M f g. C 0. : these questions may have unexpectedly similar answers for 


widely different applications. 
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Jet Propulsion News 


Alfred J. Zaehringer, American Rocket Company, Associate Editor 


Jet Aircraft, Engines 


TRANSPORTS. New version of the Boeing 707 jet 
transport is called the Intercontinental (photo). The 
plane will gross 290,000 lb and will carry 124 or 145 pas- 
sengers. The range will be over 4000 miles. Power will 
come from four P&W JT4 turbojets. Specifications in- 
clude: length, 146 ft, 8 in; span 141'/, ft; payload will be 
36,500 lb; fuel capacity is 21,200 gal. All models will be 
equipped with the new Boeing-developed combined sound 
suppressor and engine thrust reversers. Pan American 
World Airways has ordered 12 jets for delivery in August, 
1959. Other orders include: 4 for Continental Air Lines; 
Air France, 10; Sabena Belgian World Airlines, 3; and 5 
for Braniff Airways. 


Intercontinental jet transport 


@ Orders for the Douglas DC-8 jet are also rolling in. 
Eastern Airlines has ordered 26, Scandinavian Airlines 
7. Orders now total 95. Douglas has stated that al- 
though the J-57 turbojet is presently provided for in the 
DC-8, use of the Rolls-Royce Conway bypass engine is 
also under consideration for foreign use. 

@ Other jet transports have also been making news. The 
Comet 3 on a round-the-world flight is expected to cover 
30,000 miles in a flight time of about 68 hours. The 
Comet 3 is a prototype of the Comet 4 which goes into pub- 
lic service in 1959. Comet 4 will carry 60-76 passengers 
over stage lengths of 2700-3000 miles. Optimum cruising 
speed is Mach 0.74 (about 500 mph). Maximum takeoff 
weight will be around 152,000 lb. 

@ Deliveries of the Britannia, Britain’s latest turboprop 
transport, will be made in 1957 to Canadian Pacific Air- 
lines. The Britannia carries 100 passengers at 400 mph 
and at ranges up to 6000 miles. P 

@ Douglas Aircraft Co. has announced that 25 DC-8 jet 
airliners were sold to Pan American World Airways. The 
Pan Am purchases, first for the DC-8, were in excess of $160 
million. Deliveries are scheduled for 1959. Another order 
for 30 DC-8’s was placed by United Airlines. Deliveries are 
also scheduled in 1959. It is expected that flight times in 
the U. S. will be materially reduced. For example, the 
Chicago-New York run will be around 1'/; hr (about 1 hr 
less than best existing times) while the Los Angeles- 
New York trip will be around 4'/; hr (about 3 hr better than 
existing runs). 
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Norman L. Baker, Indiana Technical College, Contributor 


@ The Air Force has ordered the Convair B-58 supersonic 
bomber into initial production and the building stage is near 
at Forth Worth, Texas. 

JET ENGINES. The 500th gas turbine engine was com- 
pleted by Boeing. The engine was one of the 502-7 
series. The lastest production engines weigh 320 lb 
complete with accessories and oil. It produces 270 hp 
maximum and 240 hp on a continuous basis. Specific 
fuel consumption is about 1lb/hp/hr. Time between over- 
hauls will be about 1000 hours. An automobile equipped 
with a similar gas turbine engine, for example, would travel 
at over 100 mph and would travel over 100,000 miles 
without engine overhaul. It is anticipated that by 1957, 
the fuel consumption of these engines can be expected to 
drop to about 0.7 lb/hp/hr and making it competitive 
with gasoline engines of the same power output. A 
175-hp engine was successfully used to propell a 24-ft 
Navy personnel boat (photo). 


Gas turbine afloat 


e@ Precision forging of titanium jet turbine blades was 
announced by the Canadian Steel Improvement, Ltd. 
Secret of the process lies in protecting the titanium during 
forging and heat treatment. The process was developed 
for Orenda Engines in Malton, Ontario. The blades are 
to be used on turbojets now being produced for the 
Avro CF-100. 


deHavilland 
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FIGHTERS. New drag chutes were recently installed on 
the F-84F jet planes and have ‘given landings in only 
00 ft. This is 1000 ft shorter than the normal landing 
yithout chute. The chute is 16.ft in diameter and trails 
95 ft behind the craft. At- 130 ‘knots, the drag of the 
chute is about 8000 lb. The eatin is being installed 
on present production equipment. ‘ 


¢@ First flight of the F-105A, the USAF-Republic super- 
sonic fighter-bomber has been completed at Edwards 
AFB, Calif. The model is powered by a P&W J-57 turbo- 
jet engine of 10,000-lb thrust. The craft has short, very 
thin swept wings, a long cylindrical fuselage, wing-root.air 
soops, and a one-piece flying tail that sets low on the aft 
fuselage section. Unusual feature of the ship is a vertical 
fin on the bottom aft fuselage section for lateral stability. 
Classed as a weapons system, the F-105A can carry nu- 
clear weapons, as well as conventional bombs and rockets. 


@ Crusader, the F8U-1 jet fighter by Chance Vought is to 
be produced under a $100 million production order. The 
carried based Navy fighter has flown supersonically in 
level flight and gets its thrust from a P&W J-57-P-12. 


@ Two light fighters have been in the news. The British 
Folland Gnat has hit close to Mach 0.95 in flight tests. Ceil- 
ing of the small 7000 Ib craft is above 50,000 ft. Powered by 
a Bristol Orpheus turbojet of about 4000 Ib thrust, specs are: 
span, 22 ft, 2 in; length, 29 ft, 9 in. The French Dassault 
Mystere 22 is a NATO fighter grossing 12,000 lb. With a 
4§-deg sweepback, the craft is powered by two Turbomeca 
jets of 2200 Ib thrust each. Armament is expected to be 30- 
mm cannon, rockets, and external stores. 


@ The Soviets have disclosed that their twin-jet medium 
bomber has a range of about 3000 miles. The bomber has 
also been modified for use as an airliner. 


Facilities 

@ DuMont has consolidated its missile work in a new Missiles 
Engineering Department. Los Angeles may be the site for 
this group. 

@Sperry Rand Corp. is locating its new Air Division in 
Phoenix, Ariz., and will include work on drones and pilotless 
aircraft. 


@ Steel RATO sub-assemblies will be made by the Norris- 
Thermador Corp. in a new $2 million plant at Waco, Tex. 
Work will be concerned with the new T-60 RATO being pro- 
duced by Phillips Petroleum at McGregor, Tex. 


@ Philadelphia has been selected as the new headquarters of 
the Special Defense Projects Department of General Electric. 
The department is currently engaged in research and develop- 
ment work on several large weapons systems in the guided 
missile field. The move from Schenectady will take place in 
mid-1956 and total employment will be about 1000 personnel. 


@ Rocketdyne is the name of the rocket division of North 
America Aviation. Headquarters of the $9.5 million estab- 
lishment is taking shape at Canoga Park in the San Fernando 
Valley of Calif. Rocketdyne will be responsible for the design. 
development, manufacture, and tests of rocket power plants 
for aircraft and missiles and for the design and development 
@ small caliber aircraft rockets. The division supplies the 
ATO rocket for SM-64 NAVAHO and rocket power plants for 
ATLAS and REDSTONE. Downey Division is now work- 
ing on SM-64 (photo). All large rockets, however, are static 
fested at the Simi Hills Field Test Laboratory in the Santa 
Susana Mountains. 


® Two new supersonic wind tunnels were announced by in- 
dustry. Republic Aviation is to build a Mach 4 tunnel cost- 
ing $12 million at its Farmingdale, N: Y., facility. Boeing.is 
putting the finishing touches on its 9 X 9-in. tunnel at Wich- 
ita, Kan. This tunnel is also of Mach 4 design (photo). 


Aprin 1956 


North American 


Simulation tower for flight loads on NAVAHO 


Boeing 

Supersonic wind tunnel. Compressed air stored in tower at 

left passes instrumented test section (center) and expanded and 
discharged through tower at right 


ARS Regular Feature on Detroit T. V. 


‘Rocket Digest’? is a regular Friday night ARS program on 
Detroit Station WTVS. Featured on the show is Alfred J. Zaeh- 
ringer (left), Associate Editor of JET PROPULSION. Recently 
a live firing was made of one of Zaehringer’s medium size gas 
generator units—the smokeless, odorless, variety. The half- 
hour program is moderated by ARS member Donald J. Kenney 
(right), a member of the faculty at the University of Detroit. 
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NEW FRONTIER 


“If we had to put aman on the moon, we could do it.” 
— Overheard at an Institute of Aeronautical Sciences luncheon 


This impromptu statement was not a matter of idle 
conjecture. It was a statement of a positive and scien- 
tific fact — as provable as if he’d said the Aleutian 
Islands—and contingent only upon three prime requi- 
sites: enough time, money, and necessity. And by 
“we” he meant today’s mindpower and facilities oper- 


ating under the most advanced concepts of research 
and development. 

Those concepts as practiced at Martin today would 
be essential to the fastest possible solution of any 
complex flight systems problem now within the capac- 
ity of man to solve. 

It is this that has established Martin as one of the 
prime forces in the coming conquest of the new fron- 
tier—Space itself. 
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Modern 
Falconry 


Hughes 
Missile Bobbins. .. FALCON’s ‘“‘innards”’ 
me fashioned on these multiple-coil 
Winding machines each winding 12 coils 
@ copper wire. Coils will be used in 


Sectromagnets for control system. En- 

ire plant is air conditioned at higher than 

G@imospheric pressure to guard against 
dust 


Cradled... FALCON is about to undergo 

vibration tests suspended in this cradle. 

Shock tests are conducted in the bottom 

ofthe cradle. Note altitude test chamber 
in background 


Hughes 

Missile Row. . . Final assembly and testing 

completed, the birds are boxed for ship- 
ment to the Air Force 


1956 


Falcon Nest. ..Home of the FALCON is at the Tucson, Arizona, plant of Hughes Aircraft 


Hughes 


Co. Designed for the GAR-1 production, the plant, with 12 acres under one roof, was con- 


structed in 1952 at a cost of $8 million. 


Hughes 

Tube Tester...Capable of testing over 

600 vacuum tubes per hour, this auto- 

matic machine is loaded at left. Shock 

and vibration are among some of the tests 
simulated 


It contains over $6 million of production equip- 


Hughes 

Shaping Up...Forward sections, mag- 

nesium shells, begin to take shape on this 
sheet metal assembly line 


(Left) Hughes 
Bird in Hand...Smallest guided missile 
in production, the GAR-1 tips the scale at 
over 100 lb. Power comes from a Thiokol 
solid propellant rocket motor. Although 
having a typical canard layout with small 
forward fins and large aft supporting sur- 
faces, the forward wings are fixed and 
steering is accomplished by control sur- 
faces hinged to the trailing edge of the 
stabilizers. A round nose fairs into the 
cylindrical fuselage which terminates in a 
boat-tail housing rocket exhaust nozzle 


Thiokol 


Falcon Strikes. . . Makes direct hit on tar- 
get drone aircraft in tests 
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& res: a. ; ment and $15 million of special test equipment cle 


British Information Services 
DROOP-SNOOT. . .Second British aircraft capable of super- 
sonic speed in level flight is this Fairey Delta 2, single-seat craft 
with an elongated nose that can be lowered like a drawbridge 
to give pilot better vision for landing, take-off, and taxiing. 
Wing span is 26 ft, 10 in; length is 52 ft, 3 in. . 


British Information Services 
GROUND LEVEL EJECTION. . .British ejection seat has suc- 
cessfully worked at ground level and speeds of 130 mph. Top, 
dummy is fired from aircraft. Bottom, it parachutes to safety 


JPL 

PLUMBERS NIGHTMARE. . .Technician finishes hooking up a 

liquid rocket motor for test. Note observation window in center, 

armor plate at right protecting vital components, exhaust silencer 
at left 


University of Michigan 
HOT CAN. . .Observers watch progress on J-47 turbojet burner 
in test stand 


BUGLIKE. . . English Electric P-1 interceptor was the first Aerojet-General 
British plane to hit sonic speed in level flight. It has been flying AND AWAY WE GO. . .PV-1 airplane equipped with internally 


since 1954 and has made over 100 flights. 20 preproduction mounted rocket engines makes take-off. Two 15KS-1000 units 


models have been ordered by the RAF are being used here 
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fn this time we are particularly 
interested in interviewing 
graduate engineers and scientists 
in the fields of aerodynamics, air- 
craft structures, mechanical engi- 
neering, chemistry, chemical 
engineering, heat transfer, elec- 
tronics, systems analysis, electro- 
mechanical instrument design, 
instrumentation, metallurgy, 
nuclear physics and solid state 
physics. 

These men should be definitely 
interested in scientific research 
and development relating to the 
problems of the future. 


CALTECH 


JPL...an Established Center 
of Research and Development 


The Jet Propulsion Laboratory is a center devoted entirely to scientific 
research and development. Its prime objective is obtaining basic informa- 
tion in the engineering sciences related to missile development — and to 
explore the various phases of jet propulsion. In addition a large share of 
its program is devoted to fundamental research in practically all of the 
physical sciences. 

The Laboratory extends over more than 80 acres in the foothills of 
the San Gabriel mountains north of Pasadena. It is staffed entirely by 
personnel employed by the California Institute of Technology and conducts 
its many projects under contracts with the U.S. Government. 

Exceptional opportunity for original research coupled with ideal 
facilities and working conditions have naturally drawn scientists and engi- 
neers of a very high caliber. These men, working in harmony, are building 
a very effective task force for scientific attack on the problems of the future. 

An unusual atmosphere of friendliness and cooperation is apparent 
at the “Lab” and newcomers soon sense the warmth of their acceptance. 
New advanced projects are now providing some challenging problems — 
and good jobs for new people. 

If you would like to develop your skill and knowledge at the “Lab” 
and, at the same time, help us solve some of our problems — write us today. 


JET PROPULSION LABORATORY 


A DIVISION OF CALIFORNIA INSTITUTE OF TECHNOLOGY 
PASADENA, CALIFORNIA 
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ARS News 


One Thousand Attend Los Angeles 
Conference on Thermal Barrier, Satellites 


ETEEN SESSIONS—including one 

ARS session which was conducted 
twice in order to accommodate demand— 
featured the March 14-16 ARS-ASME 
Aviation Division Conference at the 
Hotel Statler, Los Angeles. Cosponsors 
were the ASME Instruments and Regula- 
tors Division and the Institute of the 
Aeronautical Sciences. 

Two sessions were run under ARS 
sponsorship. One of them, on Earth 
Satellites, attracted a total of some 600 
people in two repeat sessions. 

David Shonerd of Aerophysics Develop- 
ment Corporation was ARS chairman 
of the meeting. The ASME Aviation 
Division Chairman, George T. Hayes, 
Stanford Research Institute, announced 
that total registration approached 1000. 

The Aviation Division offered sessions 
on Air Cargo (market analysis, com- 
mercial applications, military research 
forecast, military applications); on the 
Thermal Thicket (human factors, struc- 
tural materials, optimum airframe struc- 
tures, accessory drives, factor of safety 
concepts, propulsion, controls) and the 
Instruments and Regulators Division on 
aircraft instrumentation. 

Two luncheons were held and the 
featured speaker at the banquet was 
C. C. Furnas, Assistant Secretary of 
Defense for Research and Development. 
Dr. Furnas, Fellow Member of ARS, 
stressed the importance of keeping the 
public well informed on the nation’s 
defense setup and of maintaining U. S. 
leadership in research. 


Leica 


HE PINCH HIT 


Summary of ARS papers: 
EARTH SATELLITES 


Project Vanguard—The IGY Earth 
Satellite (289-56), F. R. Furth, Farns- 
worth Electronics Company. Successful 
construction and launching of the satellite 
will demand going to the present limits— 
or beyond—our engineering knowledge in 
many fields. Admiral Furth discusses the 
problems involved, the techniques to be 
followed, and the roles to be assigned to 
various military and scientific organiza- 
tions concerned with implementing the 
three-year project. 

Scientific Aspects of the Satellite 
Program (288-56), W. H. Pickering, Jet 
Propulsion Laboratory, California In- 
stitute of Technology. Observations 
made in the field of upper-atmosphere re- 
search are briefly reviewed, followed by a 
discussion of the actual scientific observa- 
tions which can be taken from the earth 
satellite. Many experiments can be per- 
formed with no instruments aboard, but 
other experiments require carrying instru- 
ments of a type in which the data can be 
sent to the earth by radio. First flights 
will certainly carry experiments to meas- 
ure the kind of environment existing on 
the satellite because this factor will affect 
the operation of instruments. 


HIGH TEMPERATURE MATERIALS 


Some Aspects in the Development and 
Performance of Pure Oxide Coating: 
(Section 1), R. V. Westerholm and W. M. 
Wheildon, Norton Co.; and Application 
of Rokide A Coating to Ramjet Engine 


Stapp Delayed by Weather, Truax Pinch Hits in Wichita. 


(Wichita Eagle) 


HE WENT ON ANYWAY 


Combustion Components (Section |j) 
(295-56), Alan V. Levy, Marquardt 
Aircraft Company. Designers of rocket 
propulsion units and high performance 
aircraft turbines are approaching the 
limitations of existing materials. Section 
1 of this paper is limited to the broad line of 
approach that if we discount design im. 
provement (such as various cooling 
methods) we can employ existing ma. 
terials to their maximum low temper. 


ture physical properties by providing 
adequate surface protection and insulation 
by means of a coating. 

Section II discusses flame-sprayed alv- 
mina coating as an insulating refractory 
coating for ramjet applications. 

High Temperature—Short Lived Plastic 
Laminates (292-56), Gerald E. Dodson, 
H. I. Thompson Fiberglass Co. New 
design concepts and higher performance 
requirements have shown the need for 
developing plastic laminates capable of 
withstanding extremely high temperatures 
(2000 to 5000 F). Temperature span 
includes the melting points of most metals 
and is far beyond the long time tempera- 
ture resistance of any plastic material 
(including the irradiated materials), 
There are three families of plastic resins 
in use that are suitable both costwise and 
reproducibility-wise for use in the high 
temperature ranges. These resin families 
are discussed in detail together with 
general manufacturing techniques designed 
to utilize the best physical characteristics 
of the laminates. 

Is There a Thermal Barrier for Super- 
sonic Ramjet Tailpipes? (293-56), Walter 
Unterberg, Marquardt Aircraft Company, 
Single-walled metal tailpipes, internally 
ceramic-coated where necessary, encounter 


Attendance: 660 


Commander Robert C. Truax, national Vice President of 
ARS, was slated only to present the charter to the newly forined 
Wichita Section on Jan. 27. On his arrival, however, he found 
that the scheduled speaker, Lt. Col. John P. Stapp, renowned 
aeromedical scientist, had been delayed by weather, and 460 
dinner guests were present to hear him speak. Commander 
Truax, again demonstrating the stuff of which old time rocket 
men are made, went ‘‘on’’ and regaled his audience with the 
“old days’’ in the liquid propellant field. He is shown above 
(third from left) with John P. Gaty, Vice President of Beech 
Aircraft; James F. Reagan, Beech’s manager of missile engi- 


neering and President of the Section; and Brig. Gen. L. ©. 
Coddington, Commander of McConnell Air Force Base. 

On the following day, Col. Stapp did arrive and 200 ‘‘seconé 
nighters”’ — at McConnell Air Force Base to hear his 
lecture on his rocket sled experiments. Above photo shows 
Col. Stapp (left) with Lawrence J. McMurtrey, power plant 
staff engineer at Boeing Airplane Co., and Vice President of the 
Section; and Dr. Reagan. They are inspecting a model space 
vehicle designed and built by David Henderson, an aeronautical 
engineering student at the University of Wichita and member of 
the Section’s board of directors. 
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, thermal barrier only in the high Mach 
number-low altitude region where more 
complex cooling means are required. 
4 high-performance variable-geometry 
“paper” ramjet was chosen as example, 
yith internal gases at 4000 F throughout, 
year the maximum obtainable from hydro- 
carbon-air combustion. Graphs show how 
the tailpipe metal temperature at ther- 
mally critical locations is affected by the 
najor variables of steady state flight 
enditions, ground and air launch tran- 
sents, tail-pipe exit geometry, wall thick- 
ness, and wall material. 


New Mexico—West Texas 
Section Establishes 
Scholarship 


T° stimulate interest in rocketry and 
provide assistance to a qualifying 
junior student during his senior year at 
New Mexico College of Agriculture and 
Mechanic Arts, a $200 scholarship will be 
given by the New Mexico—West Texas 
Section. The award will be made at the 
June Commencement for the best paper 
oa subject related to the field of rocketry. 
Papers must be submitted for certification 
of standing and classification through 
the dean of the appropriate school and 
reeived by the Section by May 1, 1956. 


More Student 
Chapters Forming 


UR more student groups have sub- 

mitted by-laws and petitions to be- 
come chartered as student chapters of 
ARS. The University of Michigan stu- 
dent chapter was first to be chartered, on 
December 8, 1955. 
The new units and their officers are as 
follows: 


Academy of Aeronautics, Flushing, N. y. 
President, Alex Havriliak 


varo 
Secretary, Norman Lange 
Paculty Advisor, Ralph Hautau 


New York University, Bronz, N. Y. 
President, Ned Greenberg 
Vice President, Bertram Kramer 
Secretary, Larry Gottlieb 
urer, Andrew Dapuzzo 
Faculty Advisor, Michael Maccarone 


Polytechnic Institute of Brooklyn, Brook- 
lyn, 
President, Donald Morin: ’57 
Vice President, Mario Cardullo ’57 
Secretary, Bernard Shatz ’57 
Treasurer, J. Codispoti ’58 
Faculty Sponsor, T. P. Torda, Depart- 

ment of Mechanical Engineering 


Under the guidance of Professor Torda, 
the group is developing a liquid-fueled 
motor. Project coordinator is William 
Peschke. Design group members include 
R. Mootchnik (injector systems), Morin 
thamber configurations), I. Reba and 
P. Schuyler (exhaust nozzles), Cardullo 
propellants) and Peschike (feed systems). 
It will have a pressure feed system and is 
being designed for 250 Ib thrust. 


1956 


Vice President-Treasurer, Anthony Var- 


Charter Presented to North Texas Section 


George P. Sutton of North American’s Rocketdyne Div., national director and member- 
ship chairman of ARS, presented the charter for the newly formed North Texas Section 
on Feb. 17. Accepting the charter on behalf of the Section was George H. Craig of Con- 


vair, Fort Worth, president. 


Shown at the ceremony are (left to right) Jim Nolan, Convair, program chairman; 
John Kerr, Convair, director-at-large; Jack Wurts, Chance Vought, secretary; Sutton; 


Harry 


Graham, Chance Vought, vice president; and Craig. Charles F. Crabtree, Convair, 


treasurer, was not present when picture was taken. - 


Cardullo is editing a newsletter for the 
group, entitled ‘“Rocket Study.” 

On February 24 a meeting was addressed 
by C. C. Miesse of Aerojet-General’s 
liquid engine division on combustion 
instability. 


Saint Louis University, (Parks College), © 


East Saint Louis, Ill. 


President, James M. Llewellyn, Jr. 

Vice President, Dale Anderson 

Secretary, Frank Brown 

Treasurer, Joe Kelly 

Technical Activities Chairman, 
Butkewicz 

Faculty Sponsor, Franz Hug 


Pete 


Donald J. Morin, President, Polytechnic 
Institute of Brooklyn Chapter 


Niagara Frontier 


Elects Smith of Bell 


ILLIAM M. SMITH, chief rocket 
engineer at Bell Aircraft Corp., 
Buffalo, has been elected president of the 
Niagara Frontier Section for 1956. He 
succeeds Tommie Zannes, also of Bell. 
Following are the other officers named: 
Stanley F. Lewinski, Bell, vice president; 


John H. Keefe, Becco Chemical Div., 
Food Machinery & Chem. Corp., trea- 
surer; Daniel Y. Sing, Bell, treasurer; 
Walter J. Zebrowski, Convair, recording 
secretary. 

Directors include, for two-year terms: 
Harry Ferullo, John van Lonkhuyzen, and 
Zannes, all of Bell; and D. E. Morrison, 
Orenda Engines, Ltd. Serving for one- 
year terms will be H. O. Kauffman, 
Becco; John Blanton, Bell; George 
Markstein, Cornell Aeronautical Labora- 
tory; and Delbert D. Thomas, Olin- 
Mathieson Chemical Corp. 


Section Doings 


Arizona. A March 7 meeting, at the 
Tucson residence of C. J. Green, president, 
featured James J. Harford, ARS Execu- 
tive Secretary, who discussed the Society’s 
history and future plans. The Section 
members also talked over plans to launch 
an educational program on rockets in local 
high schools, as well as a Section publica- 
tion. 

California Desert. A meeting at the 
Naval Ordnance Test Station, China 
Lake, was attended by some 66 people. 
Heinz Haber spoke on the subject of 
Space Medicine, followed by a tour of the 
station. Richard Schmidt, Assistant 
Chief, Rocket Engine Test Station, 
Edwards AFB, is in charge of organization 
of this Section, and prospective members 
in this area may obtain application forms” 
by writing him. 

Central Texas. James J. Harford, 
ARS executive secretary, visited the 
Section on March 6. He toured the 
plant of the Phillips Petroleum Co., 
Rocket Fuels Div., in McGregor, and had 
an informal dinner meeting with Section 
officers and directors Elton E. Rush, 
Fred Sauls, A. C. Keathley, and John 
Alden, all of Phillips. 

Chicago. Local Sections of the IAS 
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Even for the special *, 
“out of this world” ‘ 
wiring problem... i 
a down-to-earth ; 
solution = 
Continental Wire 


Whether it’s a rocket to the moon or a radio for a 
room—chances are you'll find Continental ready 
fo serve your wiring needs exactly—with both 
speed and efficiency. 


With many quality insulations of Asbestos... 
Glass . . . Nylon .. . Varnished Cambric . . . 
Polyethylene . . . Polyviny! . . . Teflon . . . Zytel, 
among others, Continental also offers a wide 
range of wire sizes—in stock and on special 
order. For instance, Continental's ELECTRONIC 
HOOK-UP WIRE. This nylon-insulated, hook-up 
wire saves TIME...LABOR...and GUESS- 
WORK in assembly. Resistance to abrasion, 
acids, alkalis and petroleum solvents—and tem- 
peratures ranging from —50° C to +125° C 
—assures dependability plus versatility. Avail- 
able in AWG SIZES 18 to 32. 


One source for your many wiring requirements 
—cContinental. Write today for Continental's 
complete catalog of heat-resistant, moisture- 
resistant wires, cables and cords. Serving 600- 
5000 volts. Sizes, 18 AWG—2,000,000 CM. 


Continental's industrial wire and cable special- 
ists are available to serve you at any time. 


Contact: Continental Sales, Box 363, 
Wallingford, Conn., Phone COlony 9-7718 


ontinental 
WIRE CORPORATION 


WALLINGFORD, CONNECTICUT © YORK, PENNSYLVANIA 


and ARS joined forces for the first time at 
the February meeting at Illinois Institute 
of Technology. A pre-meeting dinner was 
held at the IIT Student Union. Speaker 
was Emerson W. Smith of North Ameri. 
can Aviation, who discussed the practi- 
cal aspects of aircraft power plant installa- 
tion testing. Several charts and photos 
were used to illustrate his talk. 

Columbus. A_ brief discussion op 
proposed by-laws preceded the February 
meeting addressed by A. N. Tifford, 
Aeronautical Engineering Dept., Ohio 
State Univ. His lecture on escape from 
earth was followed by a motion picture 
on rockets. 

Indiana. S. V. Gunn of the Rocketdyne 
Div., North American Aviation, Ine, 
spoke to approximately 75 persons at the 
February meeting. Dr. Gunn’s subject 
included a discussion of problems associ- 
ated with the development of large thrust 
rocket engines. He also discussed the 
activities of the newly formed Rocketdyne 
Div., and briefly described the test facili- 
ties located at the Santa Susana Test 
Station. 

Maryland. Sixty members and guests 
attended the February meeting held at 
Johns Hopkins University. Guest speaker 
was Kurt Stehling whose subject was 
“Winged Rockets.’’ He showed slides 
and films of the Bell X-1A and X-2 
research aircraft. 

New Mexico-West Texas. Invitations 
in the form of perforated strips of ‘Space 
tour’ tickets were issued by the Section 
in connection with their ‘Annual Solar 
System Soirée” in March. The tickets 
were “good for passage for two” from the 
White Sands Proving Ground Navy 
Officers BOQ, stops at a space station and 
six planets, and return to Earth. Each 
portion of the ticket featured an appropri- 
ate symbol of the planet named. “In- 
flight’’ punch was served to the “space 
travelers.” 


Commander L. I. Stecher (center) who 
addressed the New Mexico-West Texas 
Section recently on ‘‘The Inter-Range 
Instrumentation Group.’’ Left, program 
chairman R. Gilbert Moore, and right, sec- 
tion president R. K. Sherburne. 


New York. At the February 10th 
meeting at IBM World Headquarters, 
Gilbert Parker of IBM discussed the field 
of scientific computing with specific 
emphasis upon large electronic computing 
equipment. A film was shown (‘Piercing 
the Unknown’’) in which the use of com 
puting equipment was demonstrated i 
modern aeronautical and rocket applica 
tions. Mr. Parker discussed general 
types of modern computers and how they 
can assist the engineer and scientist 1 
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Airborne Mule Skinner 
for today’s Defense 


Now to the service of supply comes 
a new idea in military logistics— 
air transports that rival the train in 
ton-mile cost. Newest and largest, now on 
the way, is the turbo-prop Douglas C-133. 

About to undergo flight tests, C-133 packs the 
load of five freight cars into its huge hull, loads 96% 
of military or construction equipment—fully as- 
sembled and ready for action—through an adjust- 
able platform ramp in its tail section. 


Speed and range are still secret but C-133’s 
ability to shuttle back and forth across oceans 
gives it the cargo potential of a 7000-ton ship. 
Cost drops drastically because C-133 gets material 
into action in hours, rather than weeks or months. 


Biggest cargo transport—the Douglas C-133A 


Development of the Douglas C-133 shrinks supply 
lines and bolsters America’s armed strength. But the core 
of that strength is the personnel manning posts and air- 
craft. Ask your local recruiting officer about the oppor- 
tunities in the U. S. Air Force. 


First in Aviation 
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gQineers 
With 16 years leadership in the vital field 
of missile research and development, Northrop Aircraft 
offers unusual opportunities for advancement in the 
categories listed below. Where better could you be, 


and grow, than with a pioneer? There’s an interesting 
position for you in one of the following groups: 


Guidance and Controls, encompassing research and 
development of advance automatic guidance and flight 
control systems for both missiles and piloted aircraft. 
Specific areas of development include: radio and radar 
systems, flight control systems, inertial guidancé 
systems, instrument servo systems, digital computer 
and magnetic tape recording systems, airborne analog 
computer systems, optical and mechanical! systems, 
and systems test and analyzer equipment. Within these 
areas activities include: original circuit development, 
electronic and electro-mechanical design, laboratory 
and field evaluation of systems under development, and 
reliability analysis both at a system and component level. 


Flight Test Engineering Section, which plans 

the missile test programs and establishes test data 
requirements in support of the programs. The data 
requirements are predicated on the test information 
required by the Engineering analytical and design groups 
to develop and demonstrate the final missile design, 

and are the basis from which the instrumentation 
requirements are formulated. 


= 
ECORDING 


TELEMETERING 


OPTICAL 


The analysis work performed consists of aerodynamic, 
missile systems, dynamics, flight control, propulsion and 
guidance evaluation. The Flight Test Engineering 

Section is also responsible for the field test program of 
the ground support equipment required for the missile. 


Flight Test Instrumentation Section, which includes a 
Systems Engineering Group responsible for the system 
design concept; a Development Laboratory where 
electronic and electro-mechanical systems and 
components are developed; an Instrumentation Design 
Group for the detail design of test instrumentation 
components and systems; a Mechanic Laboratory where 
the instrumentation hardware is fabricated; and a 
Calibration and Test Group where the various instru- 
mentation items and systems are calibrated and tested. 


There are now a number of openings available for 
engineers in each of these groups at all experience levels. 


If you qualify for any of these challenging opportunities, 
we invite you to contact Engineering Industrial 
Relations, Plant 2, Gate 3B, Broadway & Prairie, 
Northrop Aircraft, Inc., Hawthorne, California; or write 
Manager of Engineering Industrial Relations, Northrop 
Aircraft, Inc., 1021 East Broadway, Hawthorne, Calif. 


NORTHROP AIRCRAFT, INC. 
PIONEERS IN ALL WEATHER AND PILOTLESS FLIGHT 


Producers of Scorpion F-89 Long-Range Interceptors 
and Snark SM-62 intercontinental Missiles. 


calculating problems associated with rock. 
etry. The talk was followed by a demon. 
stration of computers operated at the 
data processing center. 

Southern California. Homer E. Newell, 
Jr., acting superintendent of the Atmos. 
phere and Astrophysics Division of the 
Naval Research Laboratory, was the 
speaker at the February meeting ip 
Pasadena. His subject was the use of 
rockets and satellites for geophysical re. 
search. Dr. Newell has specialized jn 
upper atmosphere research at NRL for 
nearly ten years. 

Wichita. About 100 members and 
guests attended a March 1 meeting at the 
Beech Employees Club which featured 
the appearance of James J. Harford, 
ARS executive secretary, as speaker, 
A film program followed. 


Five Companies Join ARS 


WO Michigan organizations, two from 

California, and one eastern company 
are recent applicants for Corporate 
Membership in ARS. The firms, their 
fields of interest, and their Corporate 
representatives are: 

Accessory Products Corporation, Whit- 
tier, California. Valves, regulators, fil- 
ters and allied equipment for control of 
fluids. Research on high flow fuel, high 
and low pressure pneumatics and environ- 
mental test equipment for aircraft, rockets 
and missiles. Robert G. Rogers, Presi- 
dent; Robert L. Maple, Executive Vice 


President; William M. Smith, Vice 
President, Chief Engineer; Orville A. 
Cordes, Sales Manager; William P. 


Taylor, Project Engineer. 

Chrysler Corporation— Missile Opera- 
tions, Detroit, Mich. Producer of U. 8. 
Army’s Redstone Missile. Aerodynamics, 
ballistics, guidance and control, handling, 
launching, propulsion of missiles. C. A. 
Brady, Operating Manager; J. P. Butter- 
field, Chief Engineer; B. J. Meldrum, 
Business Manager; C. W. Williams, 
Manager, Fabrication and Assembly; 
R. M. Graham, Manager, Quality Control. 

Consolidated Western Steel Division, 
United States Steel Corporation, Los 
Angeles, Calif. Maker of motor cases for 
various rocket and missile programs. 
Charles E. Fife, contracting Engineer; 
Charles R. Fenninger, Contracting Engi- 
neer; J. Greer Thompson, Assistant Vice 
President, Sales; George E. Boyd, Con- 
tracting Manager, General Sales; M. P. 
Klick, Chief Engineer. 

Electro-Mechanical Products Company, 
Garden City, Mich. Manufacturer of 
guided missile components, flight simu- 
lators, propulsion controls, optics, and 
varied electro-mechanical devices. Jack 
M. Beauchamp, President; George T. 
Spencer, Sales Manager; Chester C. 
Quantz, Chief Engineer; ‘Wallace H. 
Franklin, Staff Engineer; Frank W. 
Beauchamp, Executive Vice President. 

Norden-Ketay Corporation, New York, 
N. Y. Airborne automatic control sy® 
tems and components. G. D. Butler, 
Director of Sales; B. Levine, Director of 
Design and Application Engineering; 
C. F. Schaeffer, Technical Director, 
Norden Laboratories Division; R. Sey- 
mour, Chief Engineer, Instrument and 
Systems Division; R. Porter, Chief 
Engineer, Precision Components Division. 
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project a proud past 
into a new realm of 


ENGINEERING OPPORTUNITY 
FOR YOU! 


Drawing on a 50-year history of engineering accomplish- 
ment, Sperry now adds a Missiles Systems Division to its ever- 
expanding organization. 

“Firsts,” of course, are an old story with Sperry engi- 
neers. From the installation of the first gyro-compass aboard a 
Navy warship, back in 1911, to a myriad of electronic wonders 
today, Sperry has been busy marking milestones of progress. 
And Sperry Engineers are eminently qualified to embark on 
their newest project. Their vast experience with missiles and 
associated systems make them intimately acquainted in this 
field. As a matter of record, way back in 1918 Sperry engi- 
neers successfully developed the first radio-controlled ‘guided 


What all this means to engineers in search of a rewarding 
life work should be clear. In Sperry’s new Missiles Systems 
Division, major opportunities are unfolding. Not only can you 
now tap the tremendous potential in the field of missiles and 
pilotless air-borne devices, but you can do so from the well- 
established base of a stable organization. Over 1500 employ- 
ees have been employed by Sperry for more than 15 years. 
And, as Sperry grows, you will grow . . . in professional stature 
and in personal gain. 


Consider These Exceptional Openings For Engineers in 
the Following Fields: 

@ Feed Back Control Systems @ Telemetering 

@ Magnetic Amplifier Circuitry @ Aerodynamics 

@ Digital Computers @ Systems Analysis 

@ Radar Systems @ Environmental Test 

@ Stable Platforms @ Reliability 
RELOCATION ALLOWANCES + LIBERAL EMPLOYEE BENEFITS 
AMPLE HOUSING in Beautiful Suburban Country Type Area 

TUITION REFUND PROGRAM (9 graduate schools in area of plant) 
MODERN PLANT with Latest Technical Facilities 

ASSOCIATION WITH OUTSTANDING PROFESSIONAL PERSONNEL 


APPLY IN PERSON 
Daily (Including Sat.). Also Wed. Eves. 


OR SUBMIT RESUME 
to Mr. J. W. Dwyer 


Engineering Employment Supervisor 


OR PHONE FOR APPOINTMENT 
Fieldstone 7-3600, Ext. 2605 or 8238 


GYROSCOPE COMPANY 


Division of Sperry Rand Corp. 
GREAT NECK, LONG ISLAND, NEW YORK 
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unusual 
career 
opportunity 
at IBM's 


AIRBORNE COMPUTER LABORATORY 


At IBM’s Vestal*, New York, Laboratory 
is gathered a group of computer engineers _ 
whose entire efforts are devoted exclusively 
to the development and perfection of 
airborne computers. 


In support of these engineers are the 
finest IBM engineering facilities, experienced 
staffs of logicians, mathematicians and 
other specialists, and an engineering know- 
how accumulated during 42 years of building 
the world’s finest mechanical, electrical 
and electronic computing equipment. 

At the Vestal Laboratory, you’ll work in 
a climate of creativeness with stimulating 
associates on today’s most advanced projects. 

If you are an engineer with experience in any 
of the following fields, we’d like to tell you more 
about this unusual opportunity and outline other 
advantages of working with IBM: digital and 
analog computer circuitry and design— 
transistor circuitry—electronic display 
systems—microwave theory and wave 
guide design—component application and 
evaluation—electronic packaging—power | 
supply design—servo and servo-control 
systems—optics—systems analysis and 
operations research—instrumentation 
theory and design—field engineering. 

Write, outlining your interests and 
qualifications to: A. N. Hurst, Room 4504, 
International Business Machines Corp. 

Endicott, New York 
‘near Binghamton and Endicott, N.Y. 


Producer of electronic 

data processing machines, 
electric typewriters, 

and electronic time equipment. 


ARS SECTIONS AND PRESIDENTS 


Alabama 
Conrad Swanson, Redstone Arsenal 


Arizona 
Charles J. Green, Hughes Aircraft Co. 


Central Texas 
Elton Rush (acting), Phillips Petroleum Co, 


Chicago 
Gerald M. Platz, Armour Research Foun. 
dation 


Cleveland-Akron 


Walter T. Olson, Lewis Flight Propulsion 
Laboratory, NACA 


Connecticut Valley 
Charles H. King, Jr., United Aircraft Corp. 


Detroit 
Fred Klemach, Nat’! Des. & Res. Corp. 


Florida 
R. S. Mitchell, Pan American World Airways 


Fort Wayne 
Lloyd Wadekamper, Indiana Tech. College 


Holl Al gordo Group 
John P. Stapp, Holl Air Develop 
Center 


Indiana 
Philip M. Diamond, Purdue University 
Maryland 
William A. Webb, Aircraft Armaments, Inc. | 


National Capital 
J. D. Gilchrist, Aerojet-General Corp. 


New England 
Joseph Kelley, Allied Research Assoc., Inc. 


Niagara Frontier 
W. M. Smith, Bell Aircraft Corp. 


New Mexico-West Texas 


Russell K. Sherburne, New Mexico College 
of A.&M.A. 


New York 
Charles J. Marsel, New York University 


North Texas 
George Craig, Convair 


Northeastern New York 

A. H. Fox, Union College 
Northern California 

M. A, Pino, California Research Corp. 
Pacific Northwest 

Jim C. Drury, Boeing Airplane Co. 
Princeton Group 

J. Preston Layton, Princeton Univ. 
St. Louis 


J. J. Mazzoni (acting), McDonnell Aircraft | 
Corp. 


- Southern California } 


Richard D. Geckler, Aerojet-General Corp. 


Southern Ohio 
W. C. Cooley, General Electric (ANP Dept.) 


Twin Cities 
T. F. Irvine (acting), University of Minnesota | 


Wichita 
J. F. Reagan, Beech Aircraft Corp. 
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MISSILE 
TEST PROJECT 


is attracting Creative Engineers 
and Scientists who are bound 
toward success 


MTP. .. the Missile Test Project of the 
RCA Service Company, Inc... . at the Air Force 
Long Range Missile Test Center . . . is the largest 
missile-testing range and laboratory in the world! 
The responsibility for providing precision 
instrumentation for advanced stage missiles with a 
vast range of performance characteristics means 

é new engineering and planning organizations and— 
Specific fields real opportunity for rapid individual growth. 


include: Here you will encounter instrumentation problems 
that challenge the state of the art and stimulate 
TELEMETRY your creative ability in Aero, Ballistic and 


Space techniques. 


RADAR 


TIMING : Your rewards will include all the fascination, 
excitement and satisfaction of achievement 
CO M M U N | CATIO NS in new frontiers of scientific knowledge... 


plus top salary. 


DATA PROCESSING 


What’s more, you’ll enjoy the pleasant climate 


OPT ics and ideal living on Florida’s central east coast, 
where Patrick Air Force Base, site of the MTP, 
° is located. 
Request Your Interview NOW. a RCA offers many additional advantages: 
Complete facilities . . . Planned advancement 
will at the program ... Professional recognition . . . Liberal 
program of company-paid bene.its . . . Relocation 


time and place you find most 
convenient. Please send complete 
resume of education and 
experience to: 


Mr. H. N. Ashby, Technical Employment 
Missile Test Project, Dept. N-7D 

RCA Service Company, Inc., P.O. Box 1226 
Melbourne, Florida 


RADIO CORPORATION of AMERICA 


assistance. You should have a Bachelor’s or } 
advanced degree in EE, ME, Physics or 
Mathematics and two or more years’ experience. 
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PRESSURE 


ALTITUDE 
e AIR SPEED 


e ACCELERATION 


e PRESSURE 
SWITCHING 


e ALTITUDE 
CONTROL 


2 WEST 


Potentiometer type 
transducers for use 
at temperatures be- 
tween — 60° C. and 
+120° C. 


Designs available 
permit single or mul- 
tiple range linear 
output from 
non-linear input 
function, and 
operation over 
wide ranges of envir- 
onmental conditions. 


_ Designs available for 


operation to 160° C. 


Special instruments 
for use with corro- 
sive oxidants and 
fuels are available. 


TO THE FINE ENGINEERING MIND 
SEEKING THE CHALLENGING PROJECTS IN 


ROCKET PROPULSION ENGINEERING 


ROCKET PROPULSION ENGINEERS are offered unusual career opportunities now 
at Convair in beautiful, San Diego, California, including: Design Engineers for 
design and analysis of advanced high performance rocket engine systems and 
components including propellant systems, lubrication systems, control systems, 
mounting structure, and auxiliary power plants; Development Engineers for 
liaison with Engineering Test Laboratories and Test Stations in the planning,. 
analysis, and coordination of rocket engine system and component tests; Devel- 
opment Engineers for coordination with Rocket Engine Manufacturers in the 
installation design, performance analysis, and development tests in conjunction 
with Convair missile programs. Professional engineering experience in rocket 
missiles and aircraft propulsion system development will qualify you for an 
exceptional opportunity. 

CONVAIR offers you an imaginative, explorative, energetic engineering depart- 
ment... truly the “engineer's” engineering department to challenge your mind, 
your skills, your abilities in solving the complex problems of vital, new, long- 
range programs. You will find salaries, facilities, engineering policies, educa- 
tional opportunities and personal advantages excellent. 


Generous travel allowances to engineers who are accepted. Write at once 
enclosing full resume to: 


H. T. Brooks, Engineering Personnel, Dept. 1416 


CONVAIR 


A Division of General Dynamics Corporation 
3302 PACIFIC HIGHWAY . SAN DIEGO, CALIFORNIA 


SMOG-FREE SAN DIEGO, lovely, sunny city on the coast of Southern California, 
offers you and your family a wonderful, new way of life ...a way of life judged 
by most as the Nation’s finest for climate, natural beauty and easy (indoor- 
outdoor) living. Housing is plentiful and reasonable. 
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balanced or single- 
ended, with sensi- 
tivity of 1 mv to 
2 v/em (AC), 1 mv 
to 2 v/mm (DC). 


VERSATILITY 


A Sanborn "150 Series” System 
can be set up to record any of these 
inputs in any of the channels, 


Low Level Signals, 
with extreme sta- 
bility; high gain, 
and greater band- 
width than with 
150-1500 Low 
Level Preamplifier. 
STABILIZED DC Preamp 


Magnitude and Direction 
of Physical Variables, 


with variable re- 
sistance, differen- 
tial transformer or 
variable reluctance 
transducers. 


Average Value of AC Watts in a Circuit, 
in ranges from 25 
volts x 40 ma to 
250 volts x 2 amps. 
(with internal multi- 
pliers and shunts 
which can handle 


AC WATIMETER Preamp UP to 4 amps). 


in phase or 180° out 
of phase with a 
reference voliage 
(e.g., servo error 
signal). 


Higher Level Signals 
where maximum 
sensitivity of 
1 v/em, and input 
impedance of about 
200,000 ohms are 
adequate. 

INPUT COUPLING NETWORK 


DC Signals 

(push-pull, single- 
ended or difference 
between two). 
Basic sensitivity 
50 mv/em te 50 
v/em. 

DC COUPLING Preamp 


RMS Values of AC Voltages, Currents, 


50 ma — 1 amp. 


trom 25-250 volts, 


VOLT /AMMETER Preamp 


Voltage Levels Recorded Logarithmically 


Audio signals (20 

cycles to'20 KC) or 

DC voltages recore'- 

ed in logarithmic 

fashion on 50 

decibel chart. 
LOG-AUDIO Preamp 


Symmetric or Asymmetric Waveform Inputs, 


in 350-450 cycles 
(2 cycles/mm) and 
375-425 cycles (1 
cycle/mm) ranges. 


FREQUENCY DEVIATION Preamp 


Extremely Low Voltages and Currents, 
at sensitivities of 
100 pv and 1 ya per 
cm. (with external 
shunt of 100 
ohms), by means 
of DC chopper 
circuit. 
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BRASIC “150” design features include: inkless 
; ing in true rectangular coordinates, — 
improved overall linearity, numerous paper | 

travel speeds, and a choice of mobile-cabinet 
or portable-case packaging in 2-, 4-, 6-, and 
-8-channel systema. 

Sanborn Representatives will be glad to help S 

_ you select the equipment best suited to your _ 

needs, Complete catalog available. 


SANBORN COMPANY, Cambridge 39, Mass. 


worlds largest 
turbine-powered 
transport helicopter 


Statham 
Accelerometers 
to gather 
invaluable 
test data 


q 

‘J, Power Supply (A) for RE G 

AC or DC Signals, | 

AC Voltage Components 

t 
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THE Y-4 PERISCOPE BOMBSIGHT used in B-47 Stratojets has 
3,433 parts, nearly 2,000 of them in this head-end assembly alone. General 


Mills manufactures 
prime contract. 


No need to swamp your staff with piece part and assembly problems... 


Eliminate the thousands of detail problems in- 
volved in turning out complete electro-mechanical 
assemblies—and save money, too! 


The Mechanical Division of General Mills is 
ready to manufacture or purchase component 
parts, assemble to your requirements, and deliver 
assemblies performance-tested to rigid Government 
standards—on time. We have the experience and 


precision instrument in quantity under a USAF 


Let General Mills supply the whole package} 


equipment required to take over the complete job. 
The highest precision standards are maintained in 
engineering, manufacture, quality control, packag- 
ing and accounting. 


LET US BID on your requirements. Write, wire 
or phone: Dept. JP-2, Mechanical Division of 
General Mills, Inc., 1620 Central Ave., Minneapolis 
13, Minn. STerling 9-8811. 


Job opportunities available for creative engineers. Work closely with outstanding men on interesting projects. 


MECHANICAL Division oF General Mills 
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Equipment 
Blectrical, Electronic 


Missile Alternator. For short-time 
operation, weighs 15 lb; 400 cycle, 12,000 
_ Can be rated to 2000 va., 3-phase 
libv. Specialty Component Motor Div., 
General Electric Co., Schenectady 5, N. Y. 


line Pressure Transducer. Rated at 
2%) psig at differential pressures up to 
(0 psid. Statham Laboratories, 12401 
W. Olympic Blvd., Los Angeles 64, Calif. 


Missile Heating Element. Neoprene 
coated, can be used as rocket-tube heater 
where only moderate heat is required. 
Higher temperatures available. Safeway 
Heat Elements, Inc., Middletown, Conn. 


Direct-Reading Pre-Amplifier. Model 
#14 is a probe type for piezoelectric 
sickups. With VTVM can read accelera- 
fion, pressure, force. Endevco Corp., 180 
E. California St., Pasadena 1, Calif. 


Teledyne Strain Gage. New pressure 
transducer, 21/,-in. diam and in. long 
available in eight ranges up to 10,000 psi. 
Taber Instrument Corp., Section 36, 111 
Goundry St., N. Tonawanda, N. Y. 


Miniature Blower. Motor and blower 
in lin. cube moves 3 cu ft of air per min. 
Weight, 1 oz, power 4 watts. For air- 
qaft and missiles. Sanders Associates, 
Ine, Nashua, N. H. 


Subminiature Relay. Length 15/1, in., 
diam 1°/;. Wt. 0.225 lb, 2 amps at 26.5 
VDC or 115 VAC. 4 pdt. ill take 50 
G@ US. Relay, 1744 Albion St., Los 
Angeles, Calif. 


Silicon Diode. Model 1N23D is de- 
signed for X-band mixer at low noise 
(15 db). Normal and reversed polarity 
models. Bomac Laboratories, Inc., Be- 
verly, Mass. 


Klystron Radar Transmitter. 50 kw 
peak, 1 kw average power for MTI 
applications. Power required is 40 amps 
per phase for 208 v 50/60 cps, 3-phase 
supply. Levinthal Electronic Products, 
2910 Fair Oaks Ave., Redwood City,Calif. 


Signal Simulator. 400 cycle suppressed 
carrier modulated signal generator for 
checking aircraft or missile autopilots. 
Micro Gee Products Inc., 6100 W. 
Slauson Ave., Culver City, Calif. 


Sub-Miniature Rate Gyro. Weighing 
402, measures 2'/, in. long and 1 in. diam. 
Rate range of 10-500 deg/sec. Fre- 
—?- of 8-75 cps. Robey Rotor Div., 

JO Rea Co. 1723 Cloverfield Blvd, 
Santa Monica, Calif. 

High Temp Pressure Switch. Kilotrol 
switch allows operation to 800 F, burst 
pressure 10,000 psi; operating pressure of 
3100 psi at 800 F. SPST '/, amp at 28 
VDC contacts. Cook Electric Co., 2700 
Southport Ave., Chicago 14, IIl. 


Compact Power Rectifier. 30 W liquid- 
cooled rectifier occupies 
Volume of 190 cu. in. Current outputs of 
540-750 maps. Input voltages of 26 v 
to 66 v rms. International Rectifier 
ig 1521 E. Grand Ave., El Segundo, 


Materials 


Metals 


No Bead Tubing. Internal welding 
i eliminated from tubing by new 
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reverse process. Smooth and flawless 
interior. Crucible Steel Co. of America, 
Oliver Bldg., Pittsburgh 30, Pa. 


Small Self- Fasteners. Size 
ranges of 0-80 and 4-48. Standard 
Pressed Steel Co., Jenkintown, Pa. 


Ni-SpanC. Nickel - chromium - iron- 
titanium alloy retains uniform elasticity 
at elevated temperatures. H. A. Wilsom 
Co., Union, N.J. - 


Safety Well. Type 304 stainless steel 
protects thermocouple elements to 3000 
psi and-—300 F to +1850 F. Conax 
a 7811 Sheridan Drive, Buffalo 21, 


JAN Shield Insert. Corrugated cad- 
mium-plated brass insert conducts heat 
away from tubes. Inserts fit -all 
and 6'/, miniature envelopes. Birtcher 
ye 4371 Valley Blvd., Los Angeles 32, 


Thermocouple Connector. Kit features 
interchangeability with various materials. 
Merlin anufacturing Corp., 12410 
Triskett Road, Cleveland 11, Ohio. 


Plastics 


F Nitric Acid Filter. Woven 
Tefion housed in stainless steel. Pore 
sizes of 150 X 300 microns and 25 xX 40 
microns. Porous Plastic Filter Co., 30 
Sea Cliff Ave., Glen Cove, N. Y 


Acid Resistant Hose. Fluorinated 
inner tube permits acid resistant opera- 
tion from —100 to +500 F. Aeroquip 
Corp., 303 8. East St., Jackson, Mich. 


Resistant Tubing. Fiberglas-epoxy 
tubing in lengths up to 9 ft and diameters 
down to '/,, in. ID. Withstands temps 
to 350 F. Lamtex Industries, 51 State 
St., Westbury, L. I., N. Y. 


Stripper for Epoxies and Polyesters. 

Solvent Monastrip EP for reclaiming 
ttings and castings. Mona Industries, 
aterson 4, N. J. 


Vinyl Chloride. Pliovic AO in organo- 
sols. Chemical Div., Goodyear Tire & 
Rubber Co., 1144 E. Market St., Akron 16, 
Ohio. 

Plastic Patch Kit. Epoxy welds metals, 
nonmetals. Temperature and acid resist- 
ant. Amchem Corp., 823 Tuxedo, High- 
land Park 3, Mich. 


Formica Tooling. Sheet S-52 for stretch 
forming, metal spinning and die forming. 
Formica Co., 4614 Spring Grove Ave., 
Cincinnati 32, Ohio. 


Polyester Resin. Solid Atlac 382 resin 
for laminates, molding permits monomer 
to be added in desired amount. Atlas 
Powder Co., Wilmington 99, Del. 


Miscellaneous 


Silicone Polymers. New silicones are 
available with water solubility. For 
rubber, plastics, and lubricants. Linde 
Air Products Co., 30 E. 42nd St., New 
York 17, N. Y. 

Thixotroping Agent. Cabo-O-Sil, a 
finely divided silicon dioxide permits lay- 
sharp corners, etc., in polyester- 
reinforced plastics. Godfrey L. 
Cabot, 77 Franklin St., Boston 10, Mass, 


SYSTEMS ENGINEERS 


IMMEDIATE 
OPENINGS 


As designers of America’s 
finest aircraft instru- 
ments, Kollsman can pro- 
vide excellent opportunity 
for professional growth. 
And the congenial at- 
mosphere, the modern fa- 
cilities contribute to an 
environment in which a 
man can do his best work. 


Please submit resumes to 
T. A. De Luca. Confi- 


dential interviews will be 
arranged. 


For development, testing and field service of airborne 
navigational systems. Some experience with servos 
utilizing electronic and electro-mechanical components 
desirable. 


ELECTRICAL ENGINEERS 


Experience in small electro-mechanical devices, motors, 
synchros and induction pick-offs. 


FIELD SERVICE ENGINEERS 


For work associated with airborne navigational equip- 
ment. Should be versed in electronic circuitry and me- 
chanical repairs. Entails some travel. 


DESIGNERS 
Electronic, Electro-mechanical, Mechanical 
1. Electronic and electrical packaging. Knowledge of 


sub-miniaturized techniques. Work is associated with 
servos, amplifiers and components. 


2. Mechanical or electro-mechanical packaging of pre- 
cision flight instruments. 
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or 


ONE OF A SERIES — depicting Rocketry — “Yesterday, Today and Tomorrow” 


from the tartar wall to the moon 


The Battle of Kai-fung-fu in 1232 ... Ft. McHenry in 1814. . . Goddard's 
experiments in 1920... the supersonic Bell X-1 in 1947 . . . these were mile- 
stones in rocketry. Slowly evolving into today’s potent weapon, the rocket is 
still more than a weapon . .. it is man’s escape route into space. 


Before this frontier can be achieved, tremendously complex problems must 
be solved. Intriguing problems that challenge the ingenuity of even the most 
creative engineer. These problems are being solved every day by Bell Aircraft’s 
engineering teams. 

Bell's rocket engine division has designed and built numerous successful 
rocket engines for different Army, Navy and Air Force projects, including the 
Nike and Rascal guided missiles. In a continuing research program complete 
facilities are available, including 44 test cells capable of conducting over 400 
test runs a month. 


For the imaginative engineer with a sincere interest in his professional 
future and a wide choice in career opportunities, contact . . . 


MANAGER, ENGINEERING PERSONNEL 


e RESEARCH ENGINEER — heat transfer and 
fluid mechanics 


e GROUP ENGINEER — heat transfer and 
fluid mechanics 


@ AERONAUTICAL AND MECHANICAL 
ENGINEERS 
e ROCKET TEST ENGINEERS 


DEPT. 20N © P.O. BOX 1 e BUFFALO, N.Y. — 
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Sintered Metal Oxides. Alite 
minum oxide features high strength, hig 
temperature and chemical _resigtan. 
Stoneware Co., Box 350, 

io. 


Silicone Adhesives. C-269 and Gg, 
pressure sensitive adhesive can be gg, 
over temperature range of —80 F to +5 
F. Dow Corning Corp., Midland, Mig 

Casting and Laminating Resin. SCs; 
73 provides technical data on new gy. 
viscosity epoxy resin, Epon 815. She 
ag Corp., 380 Madison, New You 
17, N. Y. 


Plastics Equipment. Spiral flow jp. 
tensive mixers, presses, and calender 
Stewart Bolling & Co., 3190 E. 65th & 
Cleveland 27, Ohio. 


Care and Maintenance of Aerk 
Plastics in Aircraft. Folder no. PL 
Rohm & Haas Co., Washington Squar 
Philadelphia 5, Pa. 


Sub-Miniature Galvanometers. Seria 
M galvanometers. Fluid and _ eleetn. 
magnetic damping. Flat frequency pm 
sponses of 0-3000 and 0-240 cps. Bulletin 
301, Heiland Instruments, 5200 
Ave., Denver 22, Colo. 


Product List. Alloys and compounk 
of boron, columbium, hafnium, lithium 
selenium, tantalum, tungsten, titanium 
and zirconium. Kawecki Chemical Go, 
220 E. 42nd St., New York 17, N. Y. 


Engineering Hand-Book. Data m 
Bakelite laminated, vulcanized fibre, and 
plastics. Manne-Knowlton __ Insulation 
Co., 416 W. 13th St., New York 14, N.Y. 


Commercial Aluminum Alloys. Folder 
describes properties, fabricating applics- 
tions, and sizes. Peter A. Frasse & (Co, 
17 Grand St., New York 13, N. Y. 


Titanium. Bulletin lists chemistry and 
metallurgy, properties, heat-treat, forging 
and welding of MST-6A1-4V titanium 
alloy. Mallory-Sharon Titanium Corp, 
Niles, Ohio. 


Flexible Metal Hose. Reference Man- 
ual on bronze, monel, nickel, carbon, and 
stainless steel hose. Universal Metal 
— 2133 S. Kedzie Ave., Chicago 
23, Til. 


Arc-Cast Moly. 72-p: book con- 
piles technical data. Climax Moly) 
rE ae 500 Fifth Avenue, New York 
36, N. Y. 


Small Gas Turbines. GE Brochur 
GED-2600 covers history of gas turbine: 
for small aircraft. General Electric Co, 
Schenectady 5, N. Y. 


Shaped Tubing. Data Memo No. I’ 
illustrates various shapes available for 
aircraft and special work. Superior Tube 
1505 Germantown Ave., Norristow?, 

‘a. 


High Heat Treat Steel. Manufacturing 
and engineering forum report. Lockheed 
Aircraft Corp., 2555 N. Hollywood Way, 


Burbank, Calif. 


Aluminum Coated Steel. Folder tel 
how to protect steel from corrosion 
oxidation with Alumicoat. Arthur Tickle 
Engineering Works, 21 Delavan 5&t, 
Brooklyn 31, N.Y. ~ 


PLEASE NOTE: In the March issue, 
p. 205, middle column, the manufacturer! 
name for the Liquid Oxygen Storage mode! 
was regrettably omitted. It is: Ronan an 
Kunzl, Ine., Marshall, Mich. 
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The Chance Vought F8U-1 Crusader, newest of the Navy’s 
jet fighters, is a with a Marquardt ram-air emer- 
ckage. In case of engine failure the unit 

D-C, and hydraulic 

t control and communication 


gency power pa 
pops into the air stream and supplies A 
power sufficient to maintain fli 


through all speeds down to safe landing. 


Lifesaver... wHen POWER FAILS 


Marquardt Aircraft, the West’s largest jet 
engine research and development center, is an 
important contributor to the OMAR Joint 
Technical Committee. With Reaction Motors, 
first in the American rocket industry, and Olin 
Mathieson Chemical, a leading producer of 
chemicals, explosives, metals, Marquardt adds 


greatly to a continuous joint technical effort to 
achieve improved rocket and ramjet engines 
and special fuels. Exemplifying the weapons 
systems concept in action, the OMAR Com- 
mittee combines for the first time both chem- 
ical and mechanical experience applicable to 
high-energy power generation. 


Engineers, chemists, physicists; production and tool specialists...a wide 
variety of fascinating careers await you on this weapons systems team. For 
information write OMAR Employment Officer at the company nearest you. 


PROPELLANTS 


M dt Aircraft Company Mathieson Chemical 
"Ven Mabe, California New York, New = 


Apri. 1956 


MARQUARDT AIRCRAFT 


OLIN MATHIESON CHEMICAL Slike 


REACTION MOTORS 


Ali 
* 
/ 
i 
Kcheed 
Fey 
Way, 
AMJETS 
Reaction Motors, Inc. 
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New Patents 


Electric ignition control in starting de- 
vices for turboengines (2,730,862). Rene 
L. Lamy, Paris, France, assignor to So- 
ciété du Carburateur Zenith. 

Means responsive to the speed of aturbo- 
engine during the starting period whereby 
a switch in the ignition circuit is closed as 
the speed reaches a value enabling the en- 
gine to run by its own power. 


Gaseous fuel turbine power plant having 
parallel connected compressors (2,730,- 
863). Nathan C. Price, St. Helena, Calif., 
assignor to Lockheed Aircraft Corp. 

Plurality of parallel units each including 

a compressor, a turbine and a load-driving 
turbine. A combustion chamber is in re- 
ceiving relation to the outlets of the com- 
pressor, and delivers combustion gases to 
drive the load-driving turbines. 
Circuit for measuring the average intensity 
of disturbances (2,730,896). solde R. 
de Boisblanc, Bartlesville, Okla., assignor 
to Phillips Petroleum Co. 

System for determining combustion 
instability in a jet engine: An electro- 
magnetic detector picks up acoustical 
noise from the burner and measures the 
output of an integrating circuit. 


2,730,927 


Launching device for self-propelled mis- 
siles (2,730,927). Louis Marie E. Bour- 
guard, Courbevoie, France, assignor to 

iété d’Exploitation des Materials His- 
pano-Suiza. 

Support for holding a rocket missile in 
the direction of firing. When the support 
is retracted simultaneously from all points 
of contact, the missile drops with its axis 
parallel to direction, and the propelling 
charge is immediately ignited. 

Exhaust gas temperature responsive fuel 
system for gas turbine power plant 
(2,731,974). Bruce N. Torrell, Wethers- 
field, Conn., assignor to United Aircraft 


Fuel flow metered as a function of ex- 
haust gas temperature. System includes 
rotor driven by exhaust gases and system 


for maintaining a constant pressure ratio- 


across the rotor which actuates the fuel 
metering means. 

Acoustic pulses jet engine with acoustic 
air intake system (2,731,795). Albert G. 
Bodine, Jr., Van Nuys, Calif. 

Standing wave established within a 

resonant sonic column. An air-induction 
pipe has a length corresponding to one- 
quarter of the wave length of the wave 
pattern in the fluid column. The pi 
resonates with the fluid column and the 
end of the pipe furthest from the combus- 
tion zone represents a zone of minimum 
pressure variation. 
Gas turbine plant comprising a rotary 
regenerative heat exchanger (2,731,798). 
R. J. Welsh and R. D. van Millingen, 
Rugby, England, assignors to the English 
Electric Co., Ltd. 

Process of operating a turbine plant 
with a fuel mixture of air and a,me e 
concentration below 1 percent. Me- 
chanical energy extracted from heated and 
compressed elastic fluid is used partly for 
compressing the fluid and partly for 
driving an external load. 


2,731,885 


Rocket stowage and launching mech- 
anism (2,731,885). Auldin D. Nolan, 
Palos Verdes, Calif., assignor to Northrop 
Aircraft, Inc. 

Cylindrical streamlined shell contain- 

ing 19 rockets and mounted on a wing. 
Rockets stowed at the center of the shell 
are conveyed on pivotally mounted 
launching arms to external positions as the 
shell doors are swung open. 
Control valve apparatus for gs 
engine (2,731,983). Leighton Lee II, 
Rocky Hill, Conn., assignor to Pratt & 
Whitney, Inc. 

Spool valve with barometric control for 
increasing compressor pressure differential 
moves. A valve in the system is subject 
to tailpipe temperature. Lifting of a 
piston in the system permits increasing 
the torque. 


2,732,125 


Differential area compressor bleed con- 
trol (2,732,125). Lowell E. Ruby, Marl- 
borough, Conn., assignor to United Air- 
craft Corp. 

Multistage axial flow compressor in 
which the intermediate stage is bled to 


George F. McLaughlin, Contributo, 


move a piston as a function of compressor 
pressure ratio. Discharge pressure is ad. 
mitted to the flow regulating means to 
control compressor bleeding. 

Rocket sight for aircraft (2,733,006), 
Horace W. Babcock, Pasadena, Cali 
assignor to the U. 8. Navy. 

Fire control sight with motor coupled 
sighting member and computer. A po 
tentiometer with a sliding contact is me- 
chanically coupled so that incremental 
movements of the sight are followed by 
corresponding movements of the contact. 
Ejectable aircraft seat capsule (2,733,027), 
Doane R. Gero, Akron, Ohio, assignor to 
Goodyear Aircraft Corp. 

One piece shell, completely open at 
front, on launching rails. A lower and 
upper door completely closing the shell are 
operated by gas power located within the 
capsule. 


2.733,568 


Solid propellant jet reaction motor (2,- 
773,568). Lionel A. Dickinson, Ayles- 
bury, England, assignor to the British 
Government. 

Propellant grain within a casing pro 
viding an annular space between its outer 
surface and the inside surface of the cas- 
ing. The space communicates with the 
inside of the grain at both ends and an an- 
nular resilient sealing washer prevents the 
longitudinal flow of gas through the space 
in a rearward direction. 


2,733,569 


System for supplying liquid fuel to a com- 
bustion chamber (2,733,569). David R. 
Trowbridge, Hornchurch, England, a& 
signor to the Plessey Co., Ltd. 

Pressurized exhaust gases are initially 
supplied to the combustion chamber 
from a cartridge firing device. A regulat 
ing valve conducts the flow of low pressure 
fuel from the turbopump to the chamber. 
Jet engine power regulator (2,734,340). 
Cyrus F. Wood, Swarthmore, Pa., a& 
signor to Westinghouse Electric Co. 

A normally balanced electrical network 

operated in response to engine speed and 
temperature is used for varying the fuel 
input and the area of a jet nozzle. 
Gas turbine engine with means for reheat 
ing turbine exhaust gases (2,734,341). 
Alfred Cyril Lovesey, Derby, England, 
assignor to Rolls-Royce Ltd. 

Closed end casing structure in the e& 


haust duct with a restriction maintainiig i 


the combustion chamber pressure appro 
imately at that within the main combu® 
tion equipment. 


— 


Eprror’s Nore: The patents listed above were selected from recent issues of the Official Gazette of the U. 8. Patent Office a 


Printed copies of patents may be obtained at a cost of 25-cents each, from the Commissioner of Patents, Washington 25, D.C. 
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1S Me- 
nd” ELECTRICAL CONNECTOR 
3,027), 
nor to 
AHEAVY-DUTY WATERPROOF POWER QWL outstanding features: 4. The 
AND CONTROL CONNECTOR FOR USE tt combines the strength advantages of machined bar Sable ond effecve waterproofing for” 


2. A modified, double stub thread provides for speed 3+ The cable accessory is designed 
This new QWL Bendix* Electrical Con- and convenience in mating and disconnecting and the Kellems stainless steel wire strain vrelief g Grip for addtional 


special tapered cross-section thread design resists loosen- cable locking. 
nector was designed for and is being used ing under vibration. The threads can be easily hand 6. A left-hand thread is used on the cable accessory to 
rely on ground-launching equip- ed if d by a sub such as mud _ prevent inadvertent loosening. 
missil d ground radar equi padleae 7. High-grade copper alloy contacts are used which pro- 
ees tor — r 3. An Alumilite 225 hard anodic finish is used which gives vide for high current capacity and low voltage ue. 


ment. @ case tn to the aluminum surface. This finish | The famous Bendix closed-entry socket is used for contacts 
Obviously, for this important t of offers to jon and abrasion. _sizes 12 and 16. 3 
service a the highest standa of 


design and saseciele are acceptable. 
r (2, That’s Ke will pay you to specify SCINTILLA DIV ISION OF 


Shed the Bendix Q L Electrical Connector for AVIATION CORPORATION 


*TRADEMARK 


any job that requires exceptional per- Export Sales and Service: Bendix International Division, 205 East 42nd St., New York 17, N.Y. 
pro- formance over long periods of time. Canadian Representatives: Aviation Electric Ltd., 200 Laurentien Bivd., St. Laurent, Montreal 9, Quebec 
outer 
1 the 
O 
the 


space 


ION TESTING 


Realistic simulation of flight trajectories for 
supersonic ramjet engines is now possible 
with a ‘VARIABLE SUPERSONIC NOZZLE”’ 


operat ee for the U.S. Air Force’s 
Ab DEVELOPMENT 
CENTE of contributions to aero- 
ially pr progie ss by engineers staffing the ad- 
vanted research and development facilities of 
- this $250 million wind tunnel center. 
2 Current vacancies: TURBOJET and RAMJET 
PERFORMANCE ANALYSTS . . . EXPERI- 
MENTAL RESEARCH THERMODYNAMICISTS 

» SENIOR TEST PROGRAMMERS... 
SUPERSONIC AERODYNAMICISTS. 


wite.. Lee C. Kelley 


ARNOLD ENGINEERING DEVELOPMENT CENTER 
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REVOLUTIONARY, NEW HIGH 
PRESSURE REGULATOR 


Look at these outstanding new developments: 


© Widest adjustment range available for manual 
regulation due to greater sensitivity created 
by a new engineering innovation! 

e 4” tube straight thread gasket seal bosses per 


AND10050! A izati cializing in the 
anization §| 
Maximum inlet pressure of 4500 psi! scion tania and manufac: 


Internal relief valve and filter! ture of explosive ordnance, McCor- 


. mick Selph places first emphasis on 
Internal valve and filter are removable with- the 


out disturbing spring, diaphragm structure! located in @ 60,000-s9-ft plant on a 
200-acre site having perhaps the 
most complete facilities of its kind 
in existence. 


Adjustabic Range 
Pressure 


Your procurement and reliability 
problems in specialized explosives* 
can probably be solved at McCor- 
110100-2 8/150 mick Selph — either with standard 
items, tried and proved, or with units 
produced to meet your specific need. 
110100-6 30/750 Send for data 
110100-8 150/1500 or submit your * Ignition 

problems to: Actuation 
110100-10 400/3000 Ejection 
Fracturing 


Write today for new engineering data sheet and catalog covering 
our complete range of regulators. 


| ACCESSORY PRODUCTS CORP. 
Dept.-J-4 616 W. Whittier Blvd. 


j Whittier, Calif.» Phone: OXford 3-3747 
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AND. THE LATEST 
INTERCONTINENTAL MISSILES 
these 


great guided missiles carry vital 
electronic equipment protected by 


ROBINSON 


ALL-METAL SHOCK and VIBRATION MOUNTS 


Decisions by design groups to specify Robinson stand on exacting 
performance fests. These tests compared every availabie shock and 
vibration mount to prove which had superior characteristics and which 
could best stend the terrific punishment of guided missile operation. We 
invite you to submit your problem of moximum protection and optimum 
functioning of airborne electronic equipment,” 


ROBINSON AVIATION: INC. 


SEND FOR BULLETIN Neo, 
Robinson Vibration and Shock Moana for 


"Guided Missiles, Rockets, and 4 
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Book Reviews_ 


Mathematics of Engineering Systems, 
by Derek F. Lawden, M. A., Methuen, 


London; John Wiley, New York; 1954, 


viii + 380 pages. $5.75. 
Reviewed by Joun J. Martin 
Bendix Products Division 
Bendix Aviation Corporation 


“Mathematics of Engineering Systems” 
will serve as an excellent primer for the 
young engineer beginning in the dynamic 
analysis of physical systems, and will in 
addition be the ready reference, for the 
most part, for the practitioner of the art. 
The physical systems considered range 
through electronic amplifiers and oscilla- 
tors, electric circuits, servomechanisms 
and regulators, and other mechanical and 
electrical systems. The book will be a 
blow to the feelings of some undergrad- 
uates who hold that their courses in 
mathematics and mechanics are needless 
complications to receiving a baccalaureate 
degree. 

The book begins with a concise as well 
as complete introduction to fundamentals, 
including the definition of the function 
concept, limits and derivatives, integrals, 
and complex numbers. For the person 
without the mathematical background to 
take Chap. 1 in stride, further study would 
be indicated before the latter chapters are 
taken on. For the most part, Chap. 1 
will serve as a sufficient review for the 
above average graduate engineer. 

Chap. 2, without much ado, gets into 
the practice of solving linear differential 
equations with constant coefficients by 
classical methods, and applies these to the 
solution of practical problems. As the 
next obvious step, the chapter considers 
simultaneous differential equations of the 
type mentioned and brings this considera- 
tion to useful fruition in the study of the 
stability of linear systems in general, 
and of some common systems in particular. 

Chap. 3 extends the work of Chap. 2 
on linear differential equations with con- 
stant coefficients from the point of view 
of modern methods. For the most ‘part, 
the introduction and use of the opera- 
tional calculus are the substance of the 
chapter, and this is applied to linear sys- 
tems. The chapter is concluded by the 
generalization of the Heaviside operational 
calculus to the method of the Laplace 
transform. 

Chap. 4 deals with Fourier analysis and 
is a standard treatment beginning with the 
trigonometric series representation of 
periodic functions, continuing with a 
study of how the coefficients of this series 
are obtained, and giving some useful ap- 
plications. These practical applications 
include those relating to the band width 
of an R.F. amplifier and to the servo- 
mechanism. After this, some ofthe ad- 
vantages of the Fourier analysis are pre- 
sented, and the chapter is concluded with 
a consideration of the Fourier integral in 
the infinite interval. Hence the text 
treats the operational calculus from 
Heaviside to the generalizations named for 
Laplace and Fourier. 
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Ali Bulent Cambel, Northwestern University, Associate Edifes 


In the final chapter, Chap. 5, some 
modern work applying to nonlinear dif- 
ferential equations is given and these 
range over the same problems treated 
earlier in the text. If the subject of the 
last chapter becomes of considerable in- 
terest to the reader, other texts on non- 
linear mechanics may be of interest. Two 
good ones are “Non-Linear Mechanics” 
by Minorsky and ‘Introduction to Non- 
Linear Mechanics”’ by Kryloff and Bogoliu- 
boff. 

Running through the text are problems, 
the solutions for which are given in the 
back of the book. In adition, a short 
table of common Laplace transforms is in- 
cluded. In general, the book is highly 
recommended for those entering or en- 
gaged in systems analysis work. 


A History of Mechanical Inventions, by 
Abbott Payson Usher, Harvard Uni- 
versity Press, Rev. Ed., 1954, 450 
pp. $9. 

Reviewed by R. HarTENBERG 
Northwestern University 


Recent discoveries in physics, reduced to 
practice by advances in technology, have 
somewhat shaken the world. Just where 
and how these things will find their place 
in the future remains to be seen. Thenew 
horizons are only relative, however, for the 
past has had its technological and other 
impacts, too—a fascinating story un- 
known to many. An outline of this story 
is presented in this “History of Me- 
chanical Inventions.”’ 

History, the study of past events, is, 
contrary to a generally held belief, neither 
static nor well defined, for it is the study 
of man’s unrest. Knowledge derived 
from newly found archeological evidence, 
mislaid documents and the enlightened 
collation of old materials adds contin- 
uously to the story of how it was. This 
edition, a revision that covers a span of 
25 years, reflects the winnowings of those 
years. 

Abbott Payson Usher, Professor of 
Economics (Emeritus) at Harvard Uni- 
versity, has an interest and understanding 
that goes beyond mechanical devices. 
The interplay of the forces of nature and 
their economic implications has been the 
author’s concern in four other volumes: 
it is from this matrix of ideas that the 
present volume is justly derived. One 
may not agree with all points of Professor 
Usher’s thesis, but one will find each 
thought-provoking. 

Since mechanical inventions are many, 
and perhaps even impossible to list com- 
pletely, a consideration of the chapter 
headings is necessary to show the scope of 
the over 400 pages of text. They are: 
I, The Place of Technology in Economic 
History; II, Historical Analysis of Social 
Change; III, The Particular System of 
Events; IV, The Emergence of Novelty 
in Thought and Action; V, The Early 
History of the Pure Applied Mechanical 
Sciences; VI, The Mectianical Equip- 
ment of Pre-Christian Antiquity; VII, 


The Development of Water Wheels amg 


Windmills: 150 B.C.—A.D. 1500; Vit : 
Water Clocks and Mechanical Clogkum 


16 B.C.-A.D. 1500; 
Vinci: 


IX, Leonardo 
Engineer and Inventor; X, Tig 


Invention of Printing; XI, Machinerygim 


the Textile Industries: 


100-1800; Xi 


The Development of Clocks and Watelamm 


into Instruments of Precision: 1508 
1800; XIII, The Production and AR 
plication of Power: 


15C0-1830; 


Machine Tools and Quantity Production 


1450-1850; XV, The Production aa 


Distribution of Power Since 1832. 


The first four chapters are entircly new 
Here a theory of inven 
tion and a concept of social evolution aim 
presented. Gestalt psychology forms them 


to this revision. 


basis of the argument; according to thig 


the mark of man’s behavior is the manifeaaam™ 
tation or striving to achieve an end, whem 


the solution to the problem is reached 
an act-of insight, a grasp of the logicl 


significance of the whole situation, instegiy 


of by only the blundering tactics of tril 
and error. 


The remainder of the book is in the form 
of a eritical narrative, concerned in lam 
part with the production of power. Neceam™ 


sarily interdependence of industry ani 
powef, 
studied. , 

It is inevitable that a reader will fig 


that he has questions that are left ii 
answered. In justice to this fine singagy 
volume he must realize that there Way 


space for only selected topics, and thi 
even these could not always be discussed 
in ultimate detail. The Notes and 
Bibliography are extensive; they show 
not only the tremendous wealth of sount 
material but also the great diversificatiol 
of sources. 

This book isfve 
not often that one finds the story of how 
the particular is imbedded in a total pat 
tern, ie., the relation of a part to tht 
whole. 


HALEY ELECTED CHAIRMAN 
OF ARS BOARD 


At a meeting of the ARS Board 
of Directors in New York City on 
January 24, 1956, Mr. Andrew G. 
Haley, a Director of the Society, and 
Past President of the ARS, was 
elected Chairman of the Board of 
Directors, to serve in that capacity 
until the return to active duty of 
President Noah S. Davis, Jr. Mr. 
Davis is currently unable to serve as 
Chairman of the Board because of 
illness, although his return to this 
capacity is expected in the near 
future. This election was conducted 
in accordance with the provisions 
of Article VII, Section 4, of the ARS | 
by-laws. 

Mr. Haley was also named, at the 


ARS Space Flight Committee. 
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ranges, 0-10,000 P.$.1. 


and test personne! engaged 
B® projects on rockets, pulse jets and 
turbo jets can control higher pressures 
With the new TELEDYNE electrical 
Pressure transmitter that is both 


ACCURATE and DURABLE. These . 


features make the new Teledyne the 
Most. advanced piece of equipment of 
its kind available: 15 ranges (0- 
10,000 PSI) for wide application. 
Bonded strain gage construction (in- 

itive to vibration or shock). High 


7 
’ frequency response. Linear output 


over full pressure range. Easily dis- 
assembled for clean-out and repairs. 
Insensitive to accelerations in 2 hori- 
zontal planes and vertically less than 
-3% of full scale per “G’’, Tempera- 
ture compensated for zera shift and 


liquid and gas. 
integral control 
unit for 
proportioning 
type control. 


Write for Illustrated Literature 


TABER INSTRUMENT CORPORATION 
SECTION 211 
113 Goundry St., N. Tonawanda, N.Y. 


The SAFE-EFFICIENT 
EQUIPMENT for HANDLING 


Liquefied Orygen and Witrogen 


MODEL 
LOX 150 


150 GALLON 
CAPACITY 


VACUUM INSULATED 
STORAGE AND TRANSPORT CONTAINERS 


produced by 


CRYOGENIC DIVISION 


RONAN AND KUNZL, INC. 
MARSHALL, MICHIGAN 
ILLUSTRATED LITERATURE IS AVAILABLE. WE INVITE YOUR INQUIRY. 


COMBUSTION AND 
PROPULSION RESEARCH 


Experiment Incorporated, a leading contributor to the 
development of missile propulsion systems for over 11 
years, is again expanding its operations. El’s programs 
cover a broad range of activities from fundamental re- 
search to design and fabrication of engines. 


ENGINEERS and SCIENTISTS with experience are needed 


for challenging and responsible positions in: 


@ Air-Turborocket Development @ Rocket Design and Testing 
@ Engine Component Design @ Propellant Development 
@ Fundamental Combustion Research 


Located in suburban Richmond for over a decade, the 
company offers completely modern facilities, attractive 
working conditions and opportunity for individual responsi- 


bility. Living is pleasant in Richmond and the company 
maintains competitive salaries with liberal benefit programs. 
A relocation allowance is provided new employees. 


For descriptive brochure, address inquiries to Personnel Manager. 


EXPERIMENT INCORPORATED 
RICHMOND 2, VIRGINIA 
RESEARCH - DEVELOPMENT - ENGINEERING » PRODUCTION 
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RV-A-10 MISSILE 
powered by a 
‘THIOKOL’ 
solid propellant 


rocket motor 


The RV-A-10 missile, developed jointly by 
the Army Ordnance Corps, General Electric, 
and Thiokol Chemical Corporation, is an 
example of technical progress in solid pro- 
pellant rocketry. The rocket motor for this 
missile was designed and developed by the 
Thiokol Chemical Corporation. Its success- 
ful flight tests definitely established the feasi- 
bility of “Thiokol” solid propellant power 
plants for use in large, long range and high 
altitude missiles. 


Development of the power plant for the RV- 
4-10 stems largely from Thiokol’s “‘systems” 
concept of rocket motor development. Under 
this concept of rocket motor design, all of 
the requirements of the missile are taken into 
consideration before the propulsion unit is 
put on the drawing boards. The motor’s 
thrust, duration, size and shape are all gov- 
erned by the specific application for which 
the unit is intended. 


Close coordination is conducted with the mis- 
sile prime contractor to consider problems 
such as aerodynamic stability, center of grav- 
ity requirements, guidance system impact 
and vibration resistance, wing and fin attach- 
ments, and other over-all missile problems. 


Engineers and Chemists — become a mem- 
ber of Thiokol’s rocket development and manu- 
facturing team. We welcome inquiries from 
mechanical engineers, chemical engineers and 
chemists interested in the rocket field. 


THIOKOL CHEMICAL CORPORATION 


780 North Clinton Avenue, 
Trenton 7, N. J. 


Theoko€: 


CHEMICAL CORPORATION 


ROCKET DIVISIONS AT 
ELKTON, MARYLAND + HUNTSVILLE, ALABAMA « MARSHALL, TEXAS 


Aprit 1956 
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EMPLOYMENT 
OPPORTUNITIES 


for 


For Your 


Engineers . . Scientists 


Invitation 


To Join Experienced Technical 

Our Team 

Of Solid Specialists and Trainees 
Propellant 
ROCKET 
Specialists 

Call 

PYramid 4-6775 Work and Live in 


Ask for 
Mr. Graber Smog - Free Redlands 


to 


MISSILE 
POWER 
PLANTS 


SPECIAL - 
DEVICES — 


For Challenging Careers . . . P. O. Box 111, Redlands, California 
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Jet Propulsion Engines 


Noy. 1948, 26 pp. (Declassified from 
(Confidential, Dec. 2, 1955.) 

The Supercharged and Intercooled 

free-Piston and Turbine-Compound En- 
gine, a Cycle Analysis, by A. L. London, 
ASME Paper no. 55—A-147, Nov. 1955, 
pp. 
Some Design Aspects of the Free Piston 
Gas Generator Turbine Plant. Part I, 
Thermodynamic and Component Charac- 
teristics, by S. L. Soo and W. A. Morain. 
Part II, Control and Accessories, by W. A. 
Morain and 8. L. Soo, ASME Papers no. 
55-A-146 and A—55-155, Nov. 1955, 30 
and 20 pp. 

The Slotted Blade Axial Flow Blower, by 
H. E. Sheets, ASME Paper no. 55—A-156 
Nov. 1955, 11 pp. 

Closed Systems for Testing Compres- 
sors, by R. M. Johnson, ASME Paper no. 
53—A-157, Nov. 1955, 9 pp. 

Tests on a Typical Whittle Compressor, 
by A. Simons and C. K. Roberts, Gt. Brit. 
Aeron. Res. Counc. Rep. and Mem. no. 
2913, 1955, 10 pp. 

The Gyron—Most Powerful Jet Engine 
in the World; Announcement of the 
Gyron Junior, De Havilland Gazette, no. 
§9, Oct. 1955, pp. 132-133. 


Heat Transfer and 
Fuel Flow 


An Experimental Investigation of Two- 
Dimensional Turbulent Mixing of Parallel 
Compressible Jets at Large Temperature 
Differences and Approximately Equal 
Mach Numbers, by Ronald E. Walker, 
Bumblebee Rep. no. 232, May 1955, 78 pp. 
tab) Johns Univ., Appl. Phys. 

The Viscocity of Steam, Heavy Water 
Vapor and Argon at Atmospheric Tem- 
perature Up to High Temperatures, by 
C.F. Bonilla, 8. J. Wang, and H. Weiner, 
ASME Paper no. 55—A-6, Nov. 1955, 4 
pp. 
_ Momentum Diffusion from a Slow Jet 
into a Moving Secondary, by A. S. Wein- 
stein, J. F. Osterle, and W. Forstall, 
ASME Paper no. 55—A-60, Nov. 1955, 7 
Pp. 

_Liquid Droplet Heating and Evapora- 
tion in a High Temperature Gas Stream, 
by J. W. Rizika, ASME Paper no. 55— 
A-154, Nov. 1955, 25 pp. 

The Thermodynamics of Cooled Tur- 
bines, by W. R. Hawthorne. Part I, 
The Turbine Stage; Part II, The Multi- 
stage Turbine, ASME Papers no. 55— 
i and 55—A-191, Nov. 1955, 25 and 

Pp. 

Application of Internal Liquid Cooling to 
Gas Turbine Rotors, by Sumner Alpert, 
Ralph 1. Grey, and Delson D. Drake, 
acl Paper no. 55—A-202, Nov. 1955, 

Pp. 


An Experimental Investigation of Flow 
Separation in Nozzle at High Mach Num- 
bers and Low Density, by Norman B. 
Tucker, Toronto Univ. Inst. of Aerophys. 
Rep. no. 33, Aug. 1955, 33 pp. 


Aprit, 1956 


lechnical Literature Digest 


M. H. Smith, Associate Editor, and M. H. Fisher, Contributor 
The James Forrestal Research Center, Princeton University 


Tables for Quasi-Exponential Decay of 
Explosive Shock in Air, by G. F. Kinney, 
Naval Postgrad. School, Res. Rep. no. 8, 
Nov. 1955, 10 pp. 

Engineering Relations for Heat Transfer 
and Friction in High Velocity Laminar and 
Turbulent Boundary Layer Flow Over Sur- 
faces with Constant Pressure and Tem- 
tosgp, by E. R. G. Eckert, ASME 

‘aper no. 55—A-31, Nov. 1955, 11 pp. 

Calculation of a Cascade for a Distribu- 
tion of Subsonic Mach Numbers Given as 
a Potential Function, by Robert Legendre, 
Recherche Aéron. no. 47, Sept.-Oct. 1955, 
pp. 3-9 (in French). 


Combustion 


Sampling Studies of Cool Flames, b 
K. G. Williams, J. E. Johnson, and H. W. 
Carhard, Indust. Engng. Chem., vol. 47, 
Dec. 1955, pp. 2528-2532. 

Pulsating, Pressure Generating Com- 
bustion Systems for Gas Turbines, by 
F. H. Reynst, ASME Paper no. 55—A-56, 
Nov. 1955, 12 pp. 

The Temperature Stability of the Lami- 
nar Combustion Wave, by J. F. Wehner, 
Princeton Univ. Chem. Kinetics Project, 
TN 24 (Off. Sci. Res. TN-55-426), July 
1955, 13 pp. 

The Explosive Reaction of Carbon Mon- 
oxide and Oxygen at the Second Explo- 
sion Limit in Quartz Vessels, by Alvin S. 
Gordon and R. H. Knipe, J. Phys. Chem., 
vol. 59, Nov. 1955, pp. 1160-1165. 

Spectrum of Carbon-Monoxide Oxygen 
Explosion. I. The Visible Spectrum, by 
R. H. Knipe, and Alvin 8. Gordon, J. 
ore Phys., vol. 23, Nov. 1955, pp. 2097- 

101. 

On the Origin of the Electronically 
Excited C: Radical in Hydrocarbon 
Flames, by R. E. Ferguson, J. Chem. 
Phys., vol. 23, Nov. 1955, pp. 2085-2089. 

Experimenting with Rocket System 
Stabilizers, by Y. C. Lee, A. M. Pickles, 
and C. C. Miesse, Aviation Age, vol. 25, 
Jan. 1956, pp. 110-115. 

Burning Velocities of Hydrogen-Air 
Flames, by Herman Burwasser and Robert 
N. Pease, J. Amer. Chem. Soc., vol. 77, 
Nov. 20, 1955, pp. 5806-5808. 

Integration of a Simplified Kinetic 
Model for a Hydrogen Bromine Flame, by 
Edwin 8S. James, Wisconsin Univ. Naval 
Res. Lab. CM-849, July 1955, 82 pp. 

An Investigation of High-Frequency 
Combustion Oscillations in Liquid-Pro- 

ellant Rocket Engines, by Adelbert O. 

ischler, Rudolph V. Massa, and Ray- 
mond L. Mantler, NACA RM E53B27, 
June 1953, 37 pp. (Declassified from 
Confidential, Dec. 2, 1955.) 

The Effect of Structure on Combustion 
Stability of Liquid Propellant Rockets 
During Flight and During Static Test 
Stand Firings, by Robert S. Wick, Calif. 
Inst. of Tech. Jet Propulsion Lab. Prog. 
Rep. 20-274, July 1955, 36 pp. 

Combustion Instability in Liquid Pro- 
pellant Rocket Motors, Princeton Univ. 
Dept. Aeron. Engng., Rep. 216-M (13th 
Quart. Prog. Rep.), May 1-July 31, 1955, 
23 pp. 

Spectrographic Study of a Combustion 
Wave in the Laminar Regime, by S. 
Barrére, Recherche Aéron., no. 46, July- 
Aug. 1955, pp. 15-24 (in French). 


Burning Rate Studies. II. Variation of 
Temperature Distribution with Consump- 
tion Rate for Burning Liquid Systems, by 
D. L. Hildenbrand and A. Greenville 
Whittaker, J. Phys. Chem., vol. 59, Oct. 
1955, pp. 1024-1028. 


Fuels, Propellants, 
and Materials 


Vapor- and Liquid-Phase Reactions 
Between Nitrogen Dioxide and Water, by 
M. 8. Peters and J. L. Holman, /ndust. 
Engng. Chem., vol. 47, Dec. 1955, pp. 
2536-2539. 

Nitric Acid—Nitrogen Dioxide-Water 
System, by Anderson B. McKeown and 
Frank E. Belles, Indust. Engng. Chem., vol. 
47, Dec. 1955, pp. 2540-2543. 

Stability of Dilute Alkaline Solutions of 
Hydrogen Peroxide, by W. D. Nicoll and 
A. F. Smith, Indust. Engng. Chem., vol. 47, 
Dec. 1955, pp. 2548-2554. 

Corrosion of Metal Rods by the Prod- 
ucts of Combustion of Perchlorate Propel- 
lants, by Robert W. Geene, Aberdeen 
Prov. Ground, Ballistic Res. Labs. Mem. 
Rep. no. 928, Sept. 1955, 10 pp. 

The Development of Cermets as Struc- 
tural Materials, by J. T. Norton, ASME 
Paper no. 55—A-196, Nov. 1955, 9 pp. 

A Report on Jet Fuel Stability Studies 

y C. M. Barringer and others, SAE J., 
vol. 63, Dec. 1955, pp. 39-44. 

Oxidant Resistant Coatings for Molyb- 
denum, by J. R. Blanchard, Wright Air 
Dev. Center, TR no. 55-205, June 1955, 36 
pp. 
Effect of Dissolved Oxygen on the 
Filterability of Jet Fuels for Temperatures 
Between 300° and 400° F., by Anderson B. 
McKeown and Robert R. Hibbard, NACA 
RM no. E55118, Dec. 1955, 22 pp. 

Research on the Safety Characteristics 
of Normal Propyl Nitrate, Bur. Mines, Div. 
Explosives Tech. Prog. Rep. no. 1, Sept. 
14-Oct. 31, 1955, 10 pp. 

Studies on Boron Hydrides, Ninth 
Annual Report of Investigations on Water 
Reactive Chemical Compounds, by Anton 
B. Burg, Univ. Southern Calif. Dept. Chem., 
Nov. 1955, 28 pp. 

Which Alloys for Jet Hot Spots, by S. G. 
Demirijian, Materials and Methods, vol. 
42, Oct. 1955, pp. 116-118. 

Jet and Rocket Applications of the 
‘‘New’’ Metals, by Gerald D. Johnson, J. 
Space Flight, vol. 7, Nov. 1955, pp. 1-4. 

Physical Properties of Fourteen API Re- 
search Hydrocarbons, Cy to C;;, by David 
L. Camin and Frederick D. Rossini, J. 
Phys. Chem., vol. 59, Nov. 1955, pp. 1173- 
1179. 


Instrumentation and 
Experimental Techniques 


A Method of Telemetering Temperature 
Data, by R. S. Weil, J. Instrum. Soc. 
Amer., vol. 2, Nov. 1955, pp. 502-504. 

Note on a Convenient Recording Tech- 
nique for Rapid Changes (in Tempera- 
ture), by A. D Misner, Can. J. Tech., vol. 
34, Jan. 1955, pp. 41-43. 

Guide to Instrumentation Literature, by 
W. G. Brombacher, Julian F. Smith, and 
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Effects of Variation in Fuel Pressure on 
Combustion Performance of Rectangular 
Ramjet, by Wesley E. Messing and — 
Dugald O. Black, NACA RM _ E8T28, 
| 
| | 
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From advanced research into the fundamental 
forces of the universe—gravity, nucleonics, astro- 
physics—to the launching of man’s first stepping 
stones into space itself, Martin engineering activities 
are among the most exciting in the aircraft industry 
today. 

The sky is our world, and outer space is the next 
frontier! 
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the most advanced projects now in the research and 
development stage, contact J. M. Hollyday, Dept. JP- 
04, The Martin Company, Baltimore 3, Maryland. 
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1955, 157 pp. 

A Bibliography of Available Digit 
Computer Codes for Nuclear React, 
Problems, by A. Radkowsky and R 
Brodsky, Atomic Energy Comm. ABCY. 
3078, Oct. 14, 1955, 102 pp. 

Metallurgical Problems in Design 
Nuclear Power Reactors, by Vincent P 
Calkins, Metal Prog., vol. 68, Dec, 19% 
pp. 73-76. 

Preparation of Metal Powders fy, 
Nuclear Reactors, by P. Chiotti and 
Wilhelm, Metal Prog., vol. 68, Dec. 1955 
pp. 77-80. 

An Introduction to Reactor Physics, }y 
George Leppert, J. Amer. Soc. Naw 
Engrs., vol. 67, Nov. 1955, pp. 855-874, 

Powder Metallurgy—Its Role in th 
Design of Nuclear-Power Reactors, by HI 
H. Hausner and M. C. Kells. J. Ame 
Soc. Naval Engrs., vol. 67, Nov. 1955, pp 
875-882. 

Thermonuclear Power Reactors—Ar 
They Feasible? by Hans Thirring, Ny. 
cleonics, vol. 13, Nov. 1955, pp. 62-65., 

Highlights of Reactor Technology, 194?- 
1955, by M. C. Leverett, Gen. Elec. Rep 
vol. 58, Nov. 1955, pp. 23-26. 

Materials for Atomic Plants, by BlairR 
Elder, Gen. Elec. Rev., vol. 58, Nov. 1955, 
pp. 38-41. 
~ “How to Control a Boiling Reactor, by 
John Macphee, Nucleonics, vol. 13, Dec 
1955, pp. 42-45.4@ « - = 

Structural Materials for Heat Transfer 
in Nuclear Reactors, Part I—Selecting 
Materials for Liquid Sodium Systems, hj 
R. F. Koenig and E.G. Brush, Materialsani 
Methods, vol. 42, Dec. 1955, pp. 110-114 


Astrophysics, Aerophysies, 
Space Flight 


Astronautics, by Eugene Sanger, Ae 
Digest, vol. 71, Dec. 1955, p. 21. 

The Satellite Vehicle for Communite 
tions and Navigation, by Hayward = 
Canney, Jr., Aero Digest, vol. 71, Det 
1955, pp. 40-46. 

Interplanetary Magnetic Fields ai 
Cosmic Rays, by Levett Davis, & 
Phys. Rev., vol. 100, Dec. 1, 1955, pp 
1440-1444. 

Radio Observations of the Ionosphere 
at Oblique Incidence, by J. H. Chapman 
Kenneth Davies, and C. A. Littlewood, 
a J. Phys., vol. 33, Dec. 1955, pp. 7 


Cosmic-Ray Experiments with a Prot § 


Velocity Selector, by K. W. Ogilvie, Gas 
J. Phys., vol. 33, Dec. 1955, pp. 746-198 

Studies of Small Scale Turbulent Dit 
fusion in the Atmosphere, by F. 4 
Frenkiel and I. Katz, The Johns Hopkins 
Univ. Appl. Phys. Lab. CM-852, Allg 
1955, 30 pp. 

Circulation in the Atmosphere, 
Engng. Sci. Rep. no. 1, May 1955, 43 pp 
(AP. Cambr. Res. Center, TN no. 55-651), 


A Method of Sealing and Gas-Cont® 


tioning an Orbital Duty Survival Suit, 
M. C. Paul, J. Brit. Interplan. Soc., 
14, Nov.-Dec. 1955, pp. 300-306. 

The Radio Exploration of Space, ly 
R. C. Jennison, J. Brit. Interplan. 
vol. 14, Nov.-Dec. 1955, pp. 307-314. 


The Volcanic Theory of Martian Gre@® 


Areas, by A. E. Slatter, J. Brit. Interplat 
Soc., vol. 14, Nov.-Dec. 1955, pp. 319-38. 

Orbit Lifetimes of the U.S. Artificid 
Satellites, by Harold Ketchum, J. Spat 
Flight, vol. 7, Oct. 1955, pp. 1-5. 
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AIR-TO-GROUND 
TY SYSTEM 


Transmits Combat 
Pictures on FM 


Airborne military television crams a self-contained 
transmitting station into a small reconnaissance 
plane, then flies this ever-moving station over 
unpredictable terrain. Taking these adverse con- 
ditions into account, Admiral developed an ex- 
tremely compact television system which uses FM 
transmission for the picture. It is now in produc- 
tion for the U.S. Army Signal Corps. Even under 
difficult conditions, this equipment provides ex- 
cellent definition. 

The transmitting plane, flying at approximately 
1,000 feet, would have a line-of-sight range of 25 
or 30 miles. This would enable a battle commander 
aided by a panel of TV screens, each screen show- 
ing a different sector, to coordinate military 
operations over a wide area. 

In addition, a mobile ground-to-ground TV 
system is under development. Inquire about 
Admiral’s exceptional capabilities in the field of 
military electronics. Address inquiries to: 


CORPORATION 
Government Laboratories Division, Chicago 47 


took To Admiral ror 
RESEARCH DEVELOPMENT PRODUCTION 
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Facilities Brochure describing 

Admiral plants, equipment and ex- 

perience sent on request. ‘ ‘ 
ENGINEERS: The wide scope of work in progress at 
Admiral creates challenging opportunities in the field of 
your choice. Write Director of Engineering and Research, 
Admiral Corporation, Chicago 47, Illinois. 


groups. Address Public Relations Director, Admiral Corporation, Chicago 47. 


Note: NEW COLOR SOUND FILM on Admiral A 


1956 


317 


GASDYNAMICS 
RESEARCH 
ENGINEER 


needed to supervise and for- 
mulate research programs re- 
lated to high energy gas- 
dynamic equipment such as 
shock tube and rotating wave 
machines. Duties will also in- 
clude the design and construc- 
tion of test a and the 
preparation of reports. 


cludes a hy nic 
tunnel; propulsion test — 
shock 
supersonic wind tunnels 
associated instrumentation. 


An especially attractive fea- 
ture of employment at C.A.L. 
is the system of “‘internal re- 
search” wherein promising re- 
search ideas of engineers and 
scientists are spon and 
supported by the Laboratory. 


CORNELL 


ny. Tue. 


7 
ENGINEER « 


Design 
Engine Units 
for Nuclear Powered Flight 


The next big advance in aviation is 
nuclear powered flight, a field which 
offers exceptional possibilities for 
professional achievement. General 
Electric now has a position open in 
this new field. 


Qualificati lude 2 to 6 years’ 
experience in thermodynamics and 
aerodynamic development. Work 
involves thermodynamics and fluid 
flow aspects of turbine type aircraft 
engines and components. 


Openings in Cincinnati, Ohio 
and Idaho Falls, Idaho. 
Address replies to location 
you prefer 


Aircraft Nuclear Propulsion Dept. 
GENERAL ELECTRIC 


Mr. oe J. Kelly © Att: Mr. L. A. Munther 

P 0. Box 132 P. O. Box 535 

Cincinnati, Idaho Falls, Idaho 
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TORQUE WRENCH’ 
MANUAL 


Applications 

Engineering Data 

Screw Torque Data | 
Adapter Problems 
General Principles 


wrenches used in industry 


Got problems in design 
miniaturization that remind 
you of this? Send for the 
inside info on how over 500 
types and sizes of M.P. B.'s * 


OGEG>- 


BALL BEARINGS ACTUAL SIZE 


such as these will help to 
make things easier 
and easier and easier. || ten 


*k 
MINIATURE PRECISION BEARINGS, INC. 


20 Precision Park, Keene, N. H. 


ORTHAM 


Miniature Magnetic 
Tape Recorder 


Height: 5’—Diameter: 4’ 


Weight: 41 Ibs. 


Eight channels of information can be applied 
Simultaneously to a %-inch tape in the 
Model MR-1 Recorder (illustrated). This 
self-contained unit includes transistor tim- 
ing oscillator and battery power source and 
operates under axial accelerations up to 
500 g. A pre-record carrier system permits 
direct coupling with D.C. devices such as 
thermocouples and strain gauge transducers. 
MODEL MR-1 SPECIFICATIONS: 
Sensitivity :....28 millivolts D.C. full scale 
across 40 ohms impedance 


Frequency Bandwidth:.......... 0-250 cps 
Timing: ...... 1000 cps transistor oscillator 
Applications: ....Running time, ber of 


channels and recording characteristics are 

adaptable to meet requirements of appli- 

cations—including rocket and aircraft 

flight testing, oil well logging and 

atmospheric studies. 

WRITE FOR BULLETIN NO. ...JP-104 
NORTHAM PRODUCTS INCLUDE... 


s for p , acceleration and 


electronic equip t for systems. 


NORTH AMERICAN INSTRUMENTS, INC. 


2420 North Lake Avenue © Altadena, Colif 


PRESSURE 
ACTUATED 
‘SWITCH 


Accurate operation under vibration. Withstands rocket 
fuels. Small: 4.750” L. x 1.520” Dia. overall — Light: 6 0z. 


DIAPHRAGM MATERIAL Stainless steel. 


DIA Standard aireraft fuels, oil, water-alcohol 
ee solutions, nitric acids, hydrogen peroxide, 
and all other commonly used rocket fuels or 

products of combustion. 


OPERATING RANGES Adjustable from 30-200, 200-500 psig. 
OPERATING TOLERANCES Approximately 5% of control point setting. 
PROOF PRESSURE 2000 psi. 
BURST PRESSURE 4000 psi. 
AMBIENT TEMPERATURE —65° F. to +250° F. 
VIBRATION Up to 2000 cps at 20 g. 


TS SPDT—SPST_ (normal-open or normal- 
closed) —D PST. 


N 7/16-20 per AND 10050, 5/16-24 per AND 


MOUNTING Bulkhead 


R_ Hermetically-sealed pe receptacle (mates 


The Manning, Maxwell & Moore Type 6867 switch is designed for use in aircraft 
systems to control electrical circuits by pressure of a liquid or gas. Reliable, ac- 
curate operation under vibration makes the unit particularly attractive for rocket 

licati Our lete line of airborne pressure switches covers operation from 
inches of water through 4000 psi. We believe our long experience and extensive 
facilities for developi facturing and testing pressure switches for modern 
aireraft can be helpful to you. Our engineering counsel is at your service. Write 
for details about our pressure switches, including Type 6867. 


ASK FOR DATA SHEETS 


= MANNING, MAXWELL & MOORE, INC. 
j M ad AIRCRAFT PRODUCTS DIVISION, Danbury, Conn 
Aircraft Products include: Turbojet Engine Tem 
ifiers Electronic Amplifiers © 


2 ches for Rockets, Jet Engine and Airframe Applications * Pressure Gauges 
Thermocouples © Hydraulic Valves * Jet Engine Afterburner Control Systems 


Aprin 1956 


missile 
Sub-systems 


ENGINEER 


For the Systems Engineer bent on utilizing past experience 

to ‘‘go places,’’ GE's new Aircraft Products Department offers 

a fine opportunity to guide and integrate all sub-contracting 

of the missile for bomber defense missile system .. . regulating 
and coordinating schedules, costs, specifications. He will 
also provide technical coordination between all project 
sub-systems. 


7 to 9 years’ experience in missile field including structures, 
guidance, power plant, and accessory sub-systems. 

B.S., E.E., M.E. or A.E. (minimum). 

This 3-year old GE Department has expanded at an 
unusual rate. There is a fine career potential here. New 
York ‘‘Vacationland”’ location means satisfied living, too! 


Write in confidence to: 


Mr. C. E. Irwin 
Engineering Admin. 


Aircraft Products Dept. 


600 Main Street, Johnson City, N. Y. 
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Here’s important news! Du Pont research has developed 
a new oscillographic paper designed for critical record- 
ing at intermediate frequencies. It’s Lino-Writ 2 and it 
gives you clear, high contrast records every time because 
it is a fast, fog-free, non-staining paper with exceptional 
latitude. It’s made in two types: extra-thin Type W and 
standard-weight Type B. Both types take ink or pencil 
marks readily. 
NEW SPEED 

New Du Pont Lino-Writ 2 has plenty of extra speed to 
handle critical recording of intermediate frequencies at 
high paper travel rates. Thanks to this extra speed, you 
can use lower lamp voltages with Lino-Writ 2 and mini- 
mize lamp failures during a critical test. 


WON'T CRACK OR TEAR 
You can forget about cracks or tears caused by folding 
and handling when you use Lino-Writ 2. That’s espe- 
cially true with Type W—the strongest, toughest, extra- 
thin photorecording paper ever produced for commercial 
oscillographic use. In test after test, Type W outlasted 
every other paper on the market...regardless of whether 
the paper was handled and folded wet or dry, processed 
or unprocessed. And you can handle it without a worry 


DU PONT 


OSCILLOGRAPHIC PRODUCTS 


REG. U.S. PAT.OFF. 


BETTER THINGS FOR BETTER LIVING . . . THROUGH CHEMISTRY 
320 


ANNOUNCES 
LINO-WRIT 


The NEW oscillographic paper for critical 
recording at intermediate frequencies 


right after stabilization processing. : 
Type W withstands abnormal handling and abuse ham 
cause it’s the only photorecording paper that’s maim 
from 100% rag stock. : 
SPLICE-FREE ROLLS 

To decrease the chances of paper failure at the spiitam 
during an expensive test, Lino-Writ 2 is supplied ii 
splice-free rolls up to 450’ in length at no extra charge 
FITS ALL OSCILLOGRAPHS 

New Lino-Writ 2 is available in popular widths ail 
specification numbers to meet the demands of practical 
all types of oscillographs, including Century, 
dated, General Electric, Hathaway, Heiland, Midweseagl 
and Miller. 


SEND FOR FREE BOOKLET 


We’ll welcome-the opportunity to show you how this iff 4 
portant new photorecording paper can save you tif 4 
and money and improve your records of expensive testi 3 
Simply fill out the coupon and we’ll immediately semi 
you a free booklet that gives a detailed description iim 
Lino-Writ 2, as well as other Lino-Writ photorecordilil 
papers and chemicals. It’s worth doing right now! 


E. I. du Pont de Nemours & Co. (Inc.) 
Photo Products Department 

2420-17 Nemours Building, Wilmington 98, Delaware 
Please send me information about new Lino-Writ 2, other 
photorecording papers, and Du Pont Lino-Writ processing 
chemicals. 


Name 


Firm 


Street 


City State 


(244226 
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